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ABSTRACT 

A rigorous solution is presented for the case of stiffened anisotropic cylindrical shells with general 
imperfections under combined loading, where the edge supports are provided by symmetrical or 
unsymmetrical elastic rings. The circumferential dependence is eliminated by a truncated Fourier series. 
The resulting nonlinear 2-point boundary value problem is solved numerically via the "Parallel Shooting 
Method". The changing deformation patterns resulting from the different degrees of interaction between the 
given initial imperfections and the specified end rings are displayed. Recommendations are made as to the 
minimum ring stiffnesses required for optimal load carrying configurations. 
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1. INTRODUCTION 

All modern aerospace structures are subject to weight restrictions and must satisfy severe reliability 
criteria. Thus the aerospace community is always looking for ways to produce lighter and more reliable 
structures. Because thin walled stiffened and unstiffened shells exhibit very favorable strength over weight 
ratios they are often used in aerospace applications. Initially they were made mostly of light-weight metal 
alloys but recently bonded sandwich constructions with honeycomb core and advanced composite designs 
have also been proposed. Unfortunately, thin walled shells are prone to buckling instabilities. 

Numerous investigations in the past have shown that the load carrying capacity of orthotropic and of 
anisotropic shells can be significantly affected by a variety of factors such as the presence of initial 
geometric imperfections^’^ and of different in-plane and out-of-plane boundary conditions'^, by the 
nonlinear prebuckling deformations caused by the edge constraints^®’®], by eccentricities of the load 
application point^'®] and by inelastic effects 

Measuring and recording initial geometric imperfections has become a standard experimental 
procedure both for laboratory scale^®’^] and full scale structures^]. Recently Singer and his 
coworkers^®’^] have developed an experimental technique which makes it possible to estimate the 
degree of elastic support present in a particular test set up. 

When investigating the effect of initial geometric imperfections, in the beginning it has been a standard 
practice to neglect the effects caused by the edge restraint and the different boundary conditions, in order 
to simplify the mathematical problem that had to be solved. Later Arbocz and Sechler^®’^®] and Simitses 
and his coworkers^ ^®] have presented solutions for isotropic, orthotropic and anisotropic shells with 
both axisymmetric and asymmetric imperfections under different external loads, where a rigorous 
satisfaction of the specified boundary conditions was included. 

In this paper a rigorous solution is presented for the case of stiffened layered composite shells with 
general axisymmetric and asymmetric imperfections under combined axial compression, internal or 
external pressure and torsion, where the edge supports are provided by symmetrical or unsymmetrical 
elastic rings. The analysis is based on a combination of the nonlinear Donnell type anisotropic imperfect 
shell equations^ ®’^®] with the ring equations of Cohen^]. 

2. THEORETICAL ANALYSIS 

In an effort to gain insight into the possible nonlinear interaction between elastic boundary conditions 
and the initial imperfections the following analytical investigation is carried out, whereby the elastic 
boundary conditions are modeled by attaching rings of general cross-sectional shape eccentrically at the 
shell edges. The sign convention used for shell and ring analysis is shown in Fig. 1. Whenever necessary 
the corresponding variables will be distinguished by superscripts ( ) s for shell variables and by superscripts 
( ) r for rings variables. 
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2.1 Anisotropic Shell Equations 

Using the sign convention defined in Fig. 1 and introducing an Airy stress function F such that 
^x = F >yy ’ Ny= F, xx ; N X y= -F >X y (i) 

then the Donnell type nonlinear governing equations for imperfect, anisotropic shells^ 1 9 ’ 2 °] can be written 


L A +(F) " L B +(W) = " i W,xx ' J L NL(W,W + 2W) (2 ) 

L B * (F) + *"D = ^ F ’ xx+ L NL( f - w + W) + p e (3) 


The linear operators are 


*~A *0 ~ A 22^ )-xxxx 2A 2 g( ), X xxy + ( 2A i 2 + A 66^ )’xxyy 

* , . _ * 


2A 16^ )’xyyy + A n( )>yyyy 


Lb*( ) = B 2 i< )-xxxx + ( 2B 26 - B 61^ )-xxxy + ( B n + B 22“ 2B 66^ )>xxyy 

+ ( 2B 16 - B 62^ )’xyyy + B i2^ )’yyyy 


L D*( ) “ B n( )>xxxx + 4B 16^ )'xxxy + 2 ( B i2 + 2B 66^ )>xxyy 


£ + 

+ ^ B 26^ )’xyyy + B 22^ )’yyyy 


and the nonlinear operator is 


( 4 ) 

( 5 ) 

(6) 


L NL(S.T) “ S, xx T,yy 2S, x yT, x y+ S,yyT, xx ( 7) 

Commas in the subscripts denote repeated partial differentiation with respect to the independent variables 
following the comma. The stiffness parameters A^, B^, D^, A^, ... etc. including smeared stiffeners 
are defined in Appendix A. W is the component of the initial stress free imperfection and W is the 
component of the displacement normal to the shell midsurface, here both positive inward. 

These equations, together with the appropriate boundary conditions provided by the elastic end-rings, 
govern the behavior of circular cylindrical shells 

- in the prebuckling stress- and deformation state; 

- at the limit or bifurcation point (if there is one); 

- in the postbuckling stress and deformation state. 
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2.2 Reduction to an Equivalent Set of Ordinary Differential Equations 

If the initial imperfections are represented by 

W(x,0) = tA 0 (x) + tA-|(x) cos n0 + tA 2 (x) sin n0 (8) 

where A 0 (x), A-| (x) and A 2 (x) are known functions of x = x/R, then Eqs. (2) and (3) admit separable 
solutions of the following form 

W(x,0) = t(W v + Wp e + W t ) + tw 0 (x) + tw-|(x) cos n0 + tw 2 (x) sin n0 (9) 


— P Rt 1 O 1 _ — o 

F(x,0) = — _ — { - L - L p e x - t x0 + f 0 (x) + f-|(x) cos n0 + f 2 (x) cos 2n0 

+ f 2 (x) sin n0 + f^x) sin 2n0 } 

where x = x/R and 0 = y/R. 

Assuming the axial dependence of the response to be an unknown function of x will reduce the 
stability problem to the solution of a set nonlinear ordinary differential equations. This will allow the 
rigorous enforcing of the boundary conditions. The values of the Poisson’s expansions 


W v 


'12 


a. ; 


W Pe= 


22 - . 
— p e , 



(ID 


are obtained by enforcing the circumferential periodicity condition (see Appendix B). 

Substituting the expressions for W, W and F into the compatibility Eq. (2), using some trigonometric 
identities, and finally equating coefficients of like terms, results in the following system of five nonlinear 
ordinary differential equations 


A 2*2 f o V " j -pj- B 2*1 w o V + cw o “ j p2 { W 1< W 1 + 2A l) + 2w i( w -| + 2A.,) + w 1 (w 1 + 2A.,) 

// / / / ft n 

+ w 2 (w 2 + 2A 2 ) + 2w 2 (w 2 +2A 2 ) + w 2 (w 2 + 2A 2 ) } = 0 


( 12 ) 


A 22 f 1 " p2 ( 2A i2 + A 66^1 + n 4a i 1 f 1 2nA 26 f 3 + 2n3A 16 f 3 \ ^ B 21 w 1 

- n (B 11+ B 22 - 2B 66 )w 1+ n^B 12 wi+ n(2B 26 - B 61 )w 2 - n J (2B 16 - B 62 )w 2 } 
+ cw f |^ n2 ( w 0 (w 1+ 2A.|) + w 1 (w Q + 2A q ) } = 0 


(13) 


— * IV 9 — * — * a — * — * /// q — * ' p t 9 " 

^22^2 ” (2^1 2 + Agg)f2+ 16 n A-| -jf2“ ^ n ^ 26*4 + **® n "4 "r n ^ w -|( w 1 + 2A«|) 

ft/ tt tr tt /tt /t tr 

- 2w^(w^ + 2 A.j) + w-|(w^+ 2 A^) - [ w 2 (w 2 + 2A 2 ) - 2w 2 (w 2 + 2 A 2 ) + w 2 (w 2 + 2A 2 ) ] } = 0 


( 14 ) 



J TT { B 21 w 2 V " n 2 ( B iV B 22- 2B 66 > w 2 + n 4 B* 2 w 2 - n(2B 2 * 6 -Bg 1 )w 1 + n 3 (2B*g- Bg 2 )w' } 

+ CW2 ~ II" 2 * w o (w 2 + 2A 2> + w 2 (w q + 2A q ) } = 0 


* X IV 


A 22 f 4 " 4n ^( 2A i 2 + A 66 )f 4 + 16 n 4 A 1 ’' 1 f 4 + 4nA 2 gf 2 - 16n 3 A*gf 2 - £ i. n 2 { w 2 (w 1 + 2A^ 

- 2w 2 (w 1 + 2A^ + w 2 (w 1 + 2A.,) + w 1 (w 2 + 2A 2 ) - 2w 1 (w 2 + 2A g ) + wi(w g + 2A 2 ) } = 0 

Substituting in turn the expressions assumed for W, W and F into the equilibrium Eq. (3) and applying 
Galerkin’s procedure yields the following system of three nonlinear ordinary differential equations 


,47* 


(16) 


5* ,iv 1 ti7* iv 
B 0 if„ + ^ ^ D...W 


21 f o + 2 -R D 11 W o" 2c - f o + 2c ^o + V +cn ( fi(wi + A!) + 2f 1 (w 1+ A^ + f 1 (w 1 + A^ 


// // 


+ f 3 ( w 2 + a 2> + 2f 3 ( w 2 + A 2> + f 3( w 2 + a 2 ) > - 0 


(17) 


B 21 f 1 ~ n2 ( B n + B 22 2B 66^1 + n ^ B 1 2 f 1 + n ( 2B 26~ B 61^ f 3 ~ n3 ( 2B lV B 6*2^2 


An* 


(18) 


1 t r ^ * ...iv 


+ ~2 R" ^ D 11 W 1 2r|4 ( D i2 + 2C W W 1 + n4 °22 w 1 + 4nD 16 w 2 " 4n3 °26 w 2 } ' 2c y f i + 2c7(w 1 + A^) 


3pk * 


R 


rr // 


- 2cn P e (w-j + A-j) - 4cnt(w2+ A 2 ) + 2cn^ { f Q (w-j+ A«|) + fi(w Q + A Q ) } + cn^ { f 2 (w-|+ A-j) 
+ 4f 2 (w 1 + A^) + 4f 2 (w 1 + A^) + f 4 (w 2 + A 2 ) + 4f 4 (w 2 + A 2 ) + 4f 4 (w 2 + A g ) } = 0 


B 21 f 3 V n2 ( B n + B 22~ 2B 66 )f 3 + n4B 12 f 3" n ( 2B 26 - Bg 1 )f 1 + n 3 (2B 1 g- Bg 2 )f 1 


(19) 


1 t r f\ * iV 2/r> * * " A — * — * ,,f o — jk ' D " n n 

+ ~2 "R ^ °11 W 2 ” 2n (°12 + 2C W w 2 + n D 22 w 2“ 4nD 16 w 1 + 4n °26 w 1^ “ 2c y f 3 + 2cA.(w 2 + A 2 ) 
- 2cn 2 p e (w 2 + A 2 ) + 4cnx(w 1 + A.,) + 2cn 2 { f Q (w 2+ A 2 ) + f 3 (w Q + A Q ) } + cn 2 { f 4 (w-, + A^ 

+ 4f 4 (w 1 + A.j) + 4f 4 (w^ + A^) - [ f 2 (w 2 + A 2 ) + 4f 2 (w 2 + A 2 ) + 4f 2 (w 2 + A 2 ) ] } = 0 


The nondimensional stiffness coefficients Ay*, By* and Dy* are listed in Appendix A and ( )' = d/dx. 
Notice that Eq. (12) can be integrated twice yielding 

A 22 f o" ~2 r- B 21 w o + cw o" j -p p2 { w l( w 1 + 2Ai) + w 2 (w 2 + 2A 2 ) } = 0 


( 20 ) 
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The constants of integration are set equal to zero in order to satisfy the circumferential periodicity condition 
(see Appendix B). in order to be able to use the "Parallel Shooting Method" for the numerical solution of 
the above system of nonlinear ordinary differential equations it is necessary to reformulate Eqs. (13) - (19) 
such that a single 4th order derivative appears on the left-hand side. This can be accomplished as follows: 

- with the help of Eq. (18) eliminate the term w^ v from Eq. (13); 

- with the help of Eq. (19) eliminate the term w 2 v from Eq. (15); 

iv ,f 

- with the help of Eqs. (12) and (20) eliminate the terms f Q and f Q from Eq. (17); 

iv 7/ 

- with the help of Eqs. (13) and (20) eliminate the terms f 1 and f from Eq. (18); 

iv ,r 

- with the help of Eqs. (15) and (20) eliminate the terms fg and f from Eq. (19). 

Carrying out the details makes it possible to write the resulting equations as 


jy " rr rr rt /// / rrr / 

f 1 = C 1 f 1 ” C 2 f 1 “ C 3 w i ” C 4^( w i + A.,) + C 5 w-| + C 2 oif 3 + C 2 02 f 3 + C 203 w 2 + C 204 w 2 (21) 

/ r rr rr rr rr 

+ C 2 i 1 p e (w-|+ A-j ) + C 212 t(w 2 + A 2 )+ C 6 [ w Q (w-| + 2A^+ w-,(w 0 + 2A 0 ) ] - C 7 w Q (w 1 + A-|) 

rr rr 

+ CqW 0 (w-|+ A-|) - 2Ciofi(w Q + A 0 ) - Cg [ w-)(w-j+ 2A-|) + w 2 (w 2 + 2A 2 ) ] (w-| + A-)) 

“ c 10 ( f 2^ w 1 + A l) +4f 2 ( w i + A -|) + 4f 2 (w 1 + A.,) + f 4 (w 2 + A 2 ) + 4f 4 (w 2 + A 2 ) + 4f 4 (w 2 + A 2 } 


jv 


rr rr 


V = C 11*2 _ ^12 f 2 + ^205 f 4 + ^206 f 4 + ^13 ( w -|( w 1 + 2A-j) - 2w^(w 1 + 2A^ + w-|(w.| + 2A^) ( 22 ^ 

rr r r r rr rr 

~ [ w 2 (w 2 + 2A 2 ) - 2w 2 (w 2 + 2A 2 ) + w 2 (w 2 + 2A 2 ) ] } 


jv 


rr rr 


f 3 “ C 1 f 3 _ C 2 f 3 - c 3 w 2 _ c 4^( w 2 + Ag) + C 3 w 2 - C 201 f 1 - C 202 f.| - C203 w i - C 2 q 4W 1 (23) 

' ,f rr rr rr 

+ ^21lPe( w 2 + a 2) ~ ^212 x ( w i + A -|) + ^6 f w 0 ( w 2 + 2A 2) + w 2( w o + 2A o^ J _ C 7 w 0 (w 2 +A 2 ) 

rr rr 

+ CqW 0 (w 2 + A 2 ) - 2 C-|o^ 3 (w 0 + A q ) - Cg [ w-|(wi+ 2 A-j ) + w 2 (w 2 + 2A 2 ) ] (w 2 + A 2 ) 

rr r r r rr rr rr r r r rr rr 

~ C 10 { f 4 (w-i+ A-j) + 4f 4 (w 1 + A^ + 4f 4 (w 1 + A.,) - [ f 2 (w 2 + A 2 ) + 4f 2 (w 2 + A g ) + 4f 2 (w 2 + Ag) ] } 


jy ft rrr r rf r r r rr rr 

f 4 = c 11 f 4" C 12 f 4“ C 205 f 2 “ C 206 f 2 + C 13 { w g (wi+ 2A-,) - 2w 2 (w 1 + 2A.,) + w 2 (w 1 + 2A^ 

rr rrr rr rr 

+ w.|(w 2 + 2A 2 ) - 2w 1 (w 2 + 2Ag) + w-|(w 2 + 2A 2 ) } 


(24) 
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IV „ " „ 

~ ^14 w o - G 15 w o~ G 19^-( w o + A 0 > + G 17 t w l( w l + 2A-j) + W 2 (W 2 + 2 A 2 ) ] - C 10 { w^(w-)+ 2A-|) 

' ' ' n n n , , , „ „ 

+ 2w 1 (w 1 + 2A 1 ) + w 1 (w 1 + 2A 1 ) + w 2 (w 2 + 2A 2 ) + 2w 2 (w 2 + 2A 2 ) + w 2 (w 2+ 2A 2 ) } 

- c 18 { f-(( w 1 + A l) + 2 V W 1 + A l) + f l( w l + A i) + f 3 ( w 2 + A 2) + 2fg(w 2 + A 2 ) + f 3 (w 2+ A 2 ) } < 25 ) 


IV _ " " " " r nr , 

W 1 = C 20 w 1 - C 21 W- 1 + C 2 2 f -j ~ c 24 f 1 - C 23 ?l ( w l + A -|) + c 209 w 2 + c 210 w 2 + C 207 f 3 + c 208 f 3 

+ c 213Pe( w 1 + A l) + c 214 t ( w 2 + A 2^ ' C 25 I w 0 ( w 1 + 2A l) + w l( w D + 2A Q ) ] - C 2 6W 0 (w-,+ A-|) 


+ C 27 w 0 (w-|+ A-i) - C 2 e [ wi(w i+ 2A-|) + w 2 (w 2 + 2A 2 ) ] (w-| + A-)) - 2C 2 gf-|(w 0 + A Q ) 



jV _ " " " ” / w , 

w 2 = c 20 w 2 _ C 21 w 2 + C 22 f 3~ C 24 f 3~ c 23^-( w 2 + A 2^ “ C 209 w i - c 210 w i" c 207 f i ~ c 208 f i (2 
+ c 213Pe( w 2 + a 2) “ c 214 t ( w i + a -|) " c 25 [ w 0 ( w 2 + 2A 2> + w 2( w 0 + 2A 0 ) 1 “ c 26 w 0 ( w 2 + A 2 ) 

+ c 27 w o( w 2 + a 2) " c 28 I w l( w i+ 2A-j) + w 2 (w 2 + 2A 2 ) ] (w 2 + A 2 ) - 2C 2 9f3(w 0 + A Q ) 

' c 29 { f 4 (w 1 + A 1> + 4f 4 (w 1 + A 1> + 4f 4( w i + A 1 ) - [ f 2 (w 2 + A 2 ) + 4f 2 (w 2 + A g ) + 4f 2 (Wg + Ag) ] } 


The constants C-j, C2 C214 are *' s * ec * * n Appendix C. Besides the traditional simply supported and 

clamped boundary conditions (following Refs. [3] and [13] usually designated as SS-1 through SS-4 and C- 
1 through C-4, see Appendix D for details), in this paper boundary conditions will also be modeled by 
attaching elastic rings of general cross-sectional shape eccentrically at the shell edges. 


2.3 Equilibrium Equations of the End Rings 

Cohen’s ring equations! 2 ^ are based on moderate rotations and are derived using the principle of 
virtual work (see Appendix E for details). They can be written as 


ei z u >0990~ GJu -ee + El xz( v -eee + w -eeee) ~ a ( El z + GJ )Py,ee = a3 ( aF x r ) (28a) 

El xz u >eee + ( El x + a2EA ) v .ee + ( El x w >eee~ a 2 EAw, e ) - aEi xz p y e = - a 3 (aF y r ) (28b) 

El xz u -eeee + ( El x v >eee~ a 2 EAv,e) + (El x w, 0000+ a 2 EAw) - aEl xz Py00= a 3 (aF^) (28c) 

(El z + GJ)u,00+ El xz (v,0+ w,0 6 ) + a(GJpy 00- El z p y ) = - a 3 M t r (28d) 


where the definition of the variables and the sign convention used is displayed in Fig. 2. 

In order to be able to satisfy the displacement compatibility conditions between the end-rings and the 
edges of the shell, one must express the ring displacements in the same form as the one assumed for the 
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shell displacements. Further the expansion assumed for the load terms must not only be consistent with 
the terms assumed for the displacements but they must also form a self-equilibrating force system. Thus 
the Fourier decomposition of the ring equations will be based on the following expressions 


u r = Uq+ u| cos n0 + u 2 r cos 2n0 + Ug sin n0 + u 4 r sin 2n0 
v r= V V 1 s ' n n6 + v 2 sin 2n0 + v 3 cos n6 + v 4 cos 2ne 

w r = wj+ w| cos n0 + w 2 r sin n0 

p y = p y 0 + p yi cos n0 + p y2 sin ne 


and 


fJ = F v r cos n0 + F v r cos 2n0 + F v r sin n0 + F v r sin 2n0 
a x-j X2 Xg X4 

F w = F „ r sin n0 + F «, r sin 2n0 + F., r cos n0 + F w r cos 2n0 

y yi y2 y3 y4 

fJ = F r + F r cos n0 + F_ r sin n0 
z z o Z 1 z 2 

r r r r 

M t = M, + cos n0 + M. sin n0 
t t 0 t-| t 2 


(29) 


(30) 


Substituting these expressions into Cohen’s ring equations and equating coefficients of like terms results in 
the following 


2.3.1 Separated Set of Cohen’s Ring Equations: 
For n = 0 


S r U r 
s»1 u 1 



t - _ n 


r 1 


_ r 

a 2 EA 0 

i 

O - 

III 

aF z 

z o 

1 

N 

LU 

CM 

(13 

O 


K„) 


aM* r 

l 0 


For n> 2 


S 2 

0 

0 

0 




F 2 

0 

S 3 

0 

0 

< 

«4 


F 3 

0 

0 


0 


•4 



0 

0 

0 

1 

^ LO 

</> 


U 5 


r.. 


(31) 


(32) 
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where the s[ are matrices and the u[ , f[ are vectors defined as follows 


2 ^ 


n 2 (n 2 El z+ GJ) -n 3 El xz 


n^EI 


xz 


an 2 (El z + GJ) 


-n 3 EI 


xz n‘(El x +a*EA) -n(n 2 El x+ a*EA) -anEI 


n 4 EI 


xz 


k 2/tri . n 2cA\ «/«2n . 

-n(n 2 El x + a 2 EA) n 4 El x + a 2 EA 


XZ 


an^EI 


an^(EI 2 + GJ) 


-anEI 


xz 


an El 


xz 


xz 


a 2 (El z + n 2 GJ) 


( 33 ) 


3 ~ 3 
a d 


4n 2 (4n 2 El z + GJ) 


-8n d EI 


xz 


-8n d EI 


xz 


4n 2 (El x + a 2 EA) 


(34) 


S 4 = — 

4 a 3 


n 2 (n 2 El z + GJ) 


n J EI 


xz 


n^EI 


xz 


n 3 EI 
n 4 EI 


xz 


xz 


n 2 (El x + a 2 EA) n(n 2 El x + a 2 EA) 
n(n 2 El x + a 2 EA) n 4 El x + EA 


an d (El z + GJ) 


anEI 


an^(El z + GJ) 
anEI-. 


xz 


an El 


xz 


'XZ 


an 2 EI 


xz 


a 2 (El z + n 2 GJ) 


(35) 




4n 2 (4n 2 El z + GJ) 8n 3 El xz 
8n 3 El xz 4n 2 (El x + a 2 EA) 


(36) 


Further 
■>\T = 


(U 2 ) T = [ u{ v/ w j ] 
(°i) T = l u 3 r v 3 r w 2 r pj, 2 ] 


(U^) T = [u 2 r v 2 r ] 
( U 5) T - [U 4 r v 4 r ] 


(37) 


finally 


r r \ T 
2> 


. r 
X 1 

< 

aFz i 

a^] : 

CO 

IL 

ll 

&T 

** 

aF Y2 : 

. r 
x 3 

aF r 
V3 

aF 4 

aM £ ; 

: < F 5> T 

= f aF x 4 

aF V4' 


(38) 


2.3.2 Forces and Moments Acting at the Ring Centroid 

Next one must express the line loads (F^, Fy, F^) and the torsional line moment (M^ r ) acting at the 
ring centroid in terms of the applied external loads and the stress and moment resultant of the shell wall 
attached to it. Considering the free-body diagrams of the shell edges from Fig. 3 one obtains the following 
relationships at x = 0 (at the lower edge) 
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aF x r = (R - q)N 0 + RN® 

aF y r = - (R - q t )s 0 + RN x s y 

aF z r = RH s 

aM t r = - e z RN x S + (q - e z )(R - q)N 0 + e x RH s + RM® 

Whereas considering Fig. 4 at the upper edge, at x = L, the following expressions hold 

aF x = - (R - q)N 0 - RN® 
aF y= (R - q t )S 0 - RN xy 
aF z r = - RH s 

aM t r = e z RN x S - (q - e z )(R - q)N 0 - e x RH s - RM® 


(39) 


(40) 


Notice that at the upper edge (at x = L) one obtains the same expressions as at the lower edge (at x = 0) 
except for the minus sign in front of the terms on the right hand side. Further N Q is the external 
compressive line load applied at a distance q from the shell midsurface, whereas S Q is the external 
torsional line load applied at a distance q t from the shell midsurface. 

2.3.3 Displacement and Rotation Compatibility Conditions 

Considering now the compatibility of the ring and shell displacements and rotations, with the help of 
Fig. 5a one obtains the following expressions at x = 0 (at the lower edge) 


r s s 
u = u - e z w, x 

r a s s s 
v'= _ v b - e x u, y - e z w 


7 


r q 

w = w + 


e x w, 


s 

X 



Notice that the displacement v r consists of three parts, namely 
- a rotation around the centerline of the shell 


v 


r 


a 

R 


v 


s 


- a rotation about the x-axis 


(41) 


(42a) 


r s 

v = " e z w, y 


(42b) 
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- a rotation about the z-axis 

v r= _ e x u -y (42c) 

Further, by looking at Fig. 5b, it can easily be seen that at the upper edge (at x = L) one obtains identical 
compatibility conditions. 

2.4 Fourier Decomposition in the Circumferential Direction 

Using the previously shown Fourier decomposition of the ring variables (Eqs. (29)-(30)) and the 
following Fourier decomposition of the shell variables 


u s = t(u^ + u^cos n6 + u 2 S cos 2n0 + u 3 s sin nG + u 4 s sin 2n8) 
v s = t(v Q s + v® sin nG + v 2 s sin 2nG + v 3 s cos nG + v 4 s cos 2n6) 
w s = t(W v + W p + W t + w w^cos nG + w 2 s sin nG) 
w,®= t(w,Q x + w^ s x cos nG + w 2 s x sin nG) 


and 


s Et 


N x = (N Xq + N Xi cos nG + N x ^cos 2n0 + N x ^sin nG + N x ^sin 2n0) 

2 

N ^ = (N xy 0 + Nxy^in n6 + N xy 2 sin 2ne + N xy 3 C0S ne + N xy 4 C0S 2ne ) 


H s = 


t 


— L. D (HJ+ H^cos nG + H 2 s sin nG) 


R 

st s s s 

M x = — _ D (M Xq + M^cos nG + M^sin nG) 

R 


where 


(44) 


_Et 

4c 


— ; c 2 = 3(1 - v 2 ) and G = X. 

2 R 


(45) 


one obtains upon substitution into Eqs. (39) and equating coefficients of like terms, the following 


2.4.1 Separated Form of the Forces and Moments Acting at the Ring Centroid at the Lower Edge 
(at x = 0) 

For n = 0 



™ N 0 (1 - £ ) = - Mi 
Et 2 R 



( 46 ) 
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Et 2 


xy 0 ' 


— ) = t (1 - f!) 

R ' ' R ' 


~ T 


and 


c r d c s c e 

f i = b i f i + f i 


r s 6 

where the B-| matrix and the , F 1 , F 1 column vectors are defined as 


(47) 


(48) 


Bi = 


1 0 
e x R 


(49) 


<F'j) T - [aF z r o aM t r o ] 


(50) 


(F^) T -[H 0 S m x s 0 1 


(51) 


-e,T t 


EV 


(Fp'= ‘ D[0 RXq(1-^)] where D = 

R 2 H 4c 2 


; q = 4c — q 
t 2 


(52) 


B 2 0 0 0 


0 Bg 0 0 


0 0 b 4 0 


0 0 0 B 


5 J 


(53) 


where the Bi are matrices and F[, F® are column vectors defined as 


i ’ i 


B 


2 = 


Et‘ 


0 

0 


1 

0 

0 


0 

JL(1) 2 

4c R 
1 , t,2. 


0 
0 

1 / t %2r 


_!_(_.) e x _ (__) R 
4c R x 4c R 


(54) 
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Bq = 


Et‘ 


1 0 
0 1 


(55) 


B/i 


Et" 

c 


1 

0 

0 


0 

1 

0 


0 

0 

— (— ) 2 
4c R 


0 

0 

0 


V 1 2 e 1 ( 1 > 2 r 

4c n 6x 4F ( R } R 


(56) 


B 5 = 


Et‘ 


1 0 
0 1 


(55) 


Further 


< F 2> T - [<, aF y' ^1 

( p i) T ’ [<3 aF yg 


and 


< F 3> T - I aF 4 aF y^ 

< F 5> T ’ t aF x4 <1 


(58) 


(Fo) T = [N 


2 
.s,t 


I s 

*1 

N s 
xy-| 

H, S 

M xV i 

: < f §) T 

- [N x s 2 

N S 
xy 2 

I s 

*3 

N S 
xy 3 

h 2 s 

m x S 2 ] : 

: «1> T 

- < 

N S 
xy 4 


(59) 


At the upper edge (at x = L) one obtains the same expressions except for a minus sign in front of the 
terms on the right hand side. Thus 

For n = 0 


F ( --BlF? - F i 


(60) 


For > 0 


F 2 


b 2 f® 




B3F3 


-B4F4 


(61) 




- Bcf; 


5 r 5 
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Substituting next the corresponding Fourier expansions into Eq. (41) and equating coefficients of like terms 
one obtains the following 

2.4.2 Separated Form of the Displacement and Rotation Compatibility Conditions at the Lower 
Edge (at x = 0) 

For n = 0 


U r 1 -E^® + U® 


(62) 


where the E, matrix and the u!J, U^, U® column vectors are defined as 


R e x 
0 -1 


(63) 


<u r i> T -[»d O 

< U ?) T -[»o S 


i e \T 


(Up 1 = [t(W v + w pe + W t ) 0] 


(64) 

(65) 

(66) 


and ( ),- = _( ), x 


For n > 2 




E 2 

0 0 

0 


'"S' 

“3 


0 

O 

CO 

Hi 

0 


<4 


► = 






“4 


0 

o 

m 

0 



“5 


0 

0 0 

E 5 . 


U S 

re the Ej 

are matrices and U 




R 

0 

0 

_e z 


t 

ne. 

a 

ne z 

0 

E 2 “ 







"R 

0 

0 

R 

e x 



o 

0 

0 

-1 


(67) 


( 68 ) 
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E-, = ‘ 


R 0 
2ne x a 


(69) 


E 4 = — 
4 R 


-ne> 

0 

0 


0 

a 

0 

0 


-e- 


-ne z 0 


R 

0 


-1 


(70) 


Ek = ‘ 


R 0 
-2ne x a 


(71) 


Further 


’2> T 

.M 


J , , 

Vi 1 ’ 

: (U') T - [u 2 r 

v 2 r ] 

C\J 

<U5> T * I u 4 

V 4J 


(72) 


and 


AT 


AT 


<U-)' - [u 3 - v 3 " w 2 - *£-] 


(U 2 )' -tu," v," w" w“-] ; (U|) T - [u 2 “ v 2 5 ] 


/i \T r S S- 
(U 5 ) =[u 4 v 4 ] 


(73) 


Recall that, as stated earlier (see Figures 5a and 5b), at the upper edge (at x = L) the compatibility 
conditions are identical. Thus Eqs. (62) - (73) are valid both at x = 0 and x = L. 

2.5 General Elastic Boundary Conditions 

Considering the separated set of Cohen’s ring equations (Eqs. (31 -(32)), the separated set of 
equations of the forces and moments acting at the ring centroid (Eqs. (46)-(61)) and the separated set of 
equations for the displacement and rotation compatibility conditions (Eqs. (62)-(73)), then one has in matrix 
notation 

at the lower edge (at x = 0) 
for n = 0 


(31) 
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^ + F® (48) 

U r 1 -E 1 U* + U® (62) 


for n > 2 


S r U r = F f 

i i i 


F[. Bi Ff 


i = 2,3, 4,5 


Combining these equations one obtains the general elastic boundary conditions 

at the lower edge (at x = 0) 
for n = 0 


S\ (E 1 U® + U®) = B^ + F® 


for n > 2 


s !( E i u ?) = B i F ! 


i = 2, 3, 4, 5 


at the upper edge (at x = L) 
for n = 0 


(32) 

(53) 

(67) 


(74) 


(75) 


(E 1 U® + U®) = -B^-F® (76) 

for n > 2 

s[(EjU®) = - BjF® i = 2, 3, 4, 5 (77) 

Multiplying Eqs. (74) and (75) out yields the general elastic boundary conditions at the lower edge, at x 
= 0. These can be written in two different forms, namely the 

2.5.1 General Elastic Boundary Conditions Valid in the Limit as E r 0 (free edges) 

For n = 0 


a 1>o* W v * W Pe * W t ) * a,°w 0 S j - Hq 
a 2 >o S * W v . W Pe . W t ) * a 22 w oJ ’ M xV * 


( 78 ) 
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For n > 2 


a 11 u 1 + 


a 21 u 1 + 


a 31 u 1 + 
1 s 
a 41 u 1 + 


1 s 1 s 1 s 

a 12 v 1 + a 13 w 1 + a 14 w 1,x 

Is Is Is 

a 22 v 1 + a 23 w 1 + a 24 w 1,x 
Is Is Is 

a 32 V 1 + a 33 w 1 + a 34 w 1,x 
Is Is Is 

a 42 v 1 + a 43 w 1 + a 44 w 1,x 


a 11 u 2 + 


a 21 u 2 + 


a 12 v 2 = N x 2 
a 22 v 2 = ^xy 2 


a 11 u 3 + 


a 21 u 3 + 


a 31 u 3 + 


a 41 u 3 + 


a 1 3 2 v 3 S+ a 1 3 3 w 2 S+ a 1 3 4 w 2 S x 

a 22 v 3 S+ a 23 w 2 S+ a 2 3 4 w 2 S x 
3 s 3 s 3 s 

a 32 V 3 + a 33 w 2 + a 34 w 2,x 
3 s 3 s 3 s 

a 42 v 3 + a 43 w 2 + a 44 w 2,x 


a 11 u 4 + 


a 21 u 4 + 


4 s .. s 
a 12 v 4 " N x 4 

a 22 v 4 ~ N xy 4 


2.5.2 General Elastic Boundary Conditions Valid in the Limit as E r °° (fully clamped edges) 


For n = 0 


b 11 H o + b 12< M x 0 + ^ w o + w v + w p e + W t 

b 21 H o + b 22< M x n + W = w ox 


b 11 N x-| + b l2 N xy-| + b l 3 H 1 + b 14 M x-| = u i S 
b 21 N x-| + b 22 N xy-| + b 23 H 1 + b 24 M x-( = V -| S 
b 31 N x-( + b 32 N xy-) + b 33 H l + b 34 M x-| = w -f 
b 41 N x-j + b 42 N xy-( + b 43 H 1 S+ b 44 M x-| = w -| S ,x 
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.2 .. S 2 .. S s 

b 1 lV b 12 N xy 2 = u 2 

b 2 2 1 N x 2 + b l2< 2 = V 2 S 


K 3 M S ^ K 3 M S i 3 I. s . 3 s S 
b 11 N x 3 b 12 N xy 3 + b l3 H 2 + b 14 M x 2 u 3 

l 3 M s 3 .. s ,3 u s .3..s s 
b 21 N x 3 + b 22 N xy 3 + b 23 H 2 + b 24 M x 2 " v 3 

.3 m s u 3 k, s ,3 u s ,3,.s s 
b 31 N x 3 + b 32 N xy 3 + b 33 H 2 + b 34 M x 2 _ w 2 

b 41 N x 3 + b 42 N xy 3 + b 43 H 2 + b 44 M x 2 = w 2,x 


u4 M s 4 s s 
b 11 N x 4 + b 12 N xy 4 = u 4 

b 21 N x 4 + b 22 N xy 4 = v 4 


Notice that in turn the general elastic boundary conditions at the upper edge, at x = L, are obtained by 
multiplying Eqs. (76) - (77) out. These operations yield the same type of expressions as Eqs. (78) - (87), 
with the exception that all right hand terms are preceded by a negative sign. It turns out that it is 

k 

convenient to incorporate these negative signs in the definition of the ring stiffness coefficients a. - and the 

k 

ring flexibility coefficients b» , respectively. The components of the boundary stiffness and of the boundary 
flexibility matrices are listed in Appendix F. 

2.6 Reduced General Elastic Boundary Conditions 

Next the general elastic boundary conditions must be expressed in terms of the variables used in the 
anisotropic shell analysis. Recalling from Eq. (44) that 

s Et 3 s s s s s 

N x = __ (N Xq + N Xi cos n0 + N^cos 2n9 + N Xg sin n0 + N^sin 2n0) = F,yy (88) 

then substituting for F,yy from Eq. (10) and equating coefficients of like terms one gets 


" n2f 1 


N xV " n2f 3 


- 4n 2 f 2 


- -"A 
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Recalling further from Eq. (44) that 


N xy = "Hr {N xy 0 + N xy 1 sin ne + N xy 2 sin 2ne + N xy 3 C0S n0 + N xy 4 C0S 2n0 ) = 
then substituting for F >X y from Eq. (10) and equating coefficients of like terms yields 


N S = T 
'Vo T 


N xyi = nf 1 
N xy 2 = 2nf 2 


N xy 3 ” nf 3 
N xy 4 = ~ 2nf 4 


where ( ) = ( )_ and x = 21 . 

x R 

Using the expression derived in References [6] and [20] for the transverse shear force H s one can write 

H S= R^ D(H ° S+ HlSc0S 06 + "*** n0) = Mx ^ + (M) ^ + (W ’ x+ (w 7 + W.y )N x S y (92) 

From the semi-inverted constitutive equation (see Appendix A for details) 


M x,x < C 11 N x + *^12^y + ^16 N xy + ^ll K x + D i2 K y + ^16 K xy^’X 

= C 11 F -xyy + c i2 F ’Xxx' c i6 F -xxy- D n w >xxx- D i*2 W, xyy~ 2D i6 W *xxy 


(M xy + M yx ), y" 2(C 61 N x + C 62 N y + C 66 N xy + D 16 K x + D 26 K y + D 6*6 K xy)>y g4 

= 2C 61 F ’Wy + 2C 62 F, xxy- 2C 66 F ’xyy- 2D i*6 W -xxy~ 2D 26 W -yyy' 4D 66 W ’xyy 

Substituting for W, W and F from Eqs. (8)-(10), regrouping and equating coefficients of like terms one 
obtains 


"o* 2 t c lVo - D 1*1 w o" 4c t Mw o + A o> 

-2c .5. n 2 [ (w^ + A.|)fi+ (w 2 + A 2 )f 3 + (w-|+ A-))f.| + (w 2 + A 2 )fg ] 

H 1 ' 2 — f ^12 f 1 + ^62" ^i6) nf 3~ ^11” 2 ^66^ p2f 1" 2 ^61 n3f 3 1 " 4c y Mw 1 + A.,) 

- [ D^w^ + 4D^gnw 2 ~ (D^ 2 + 4Dgg)n 2 w^- 2D 2 gn 2 w 2 ] + 4c -!1 nt(w 2 + A 2 ) 

- 4c _ n 2 [ (w Q + A 0 )f-|4 2(w^ + A 1 )f 2 + 2(w 2 + A 2 )f 4 + (w-|+ A-|)f 2 + (w 2 + A 2 )f 4 ] 


( 96 ) 
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S R * j|c ~~~ 3^ ^ ^ ^ ^ Q ^ jJ, Q p% / f 

H£ = 2 _ [ “ (2002 ~ ^16^1" ^11" ^^66^ n *3 + 2 ^61 n *1 ] “ 4c — A,(w 2 + ^) 

— *'"—* " — * — * 9 ' — * O R — 

- [ D 11 W 2 ■ 4D i 6 nw 1 ~ ( D 12 + 4D 66 )n w 2 + 2D 26 n W 1 1 “ 4c — nt ( w 1 + A l) 

- 4c Ji n 2 [ (w 0 + A 0 )f 3 + 2(w 1 + A 1 )f 4 - 2(w 2 + A 2 )f 2 + (w 1 + A-,)f 4 - (w 2 + A 2 )f g ] 


(97) 


But now from Eq. (20) 


irr 



1 

2 


lliiw 





0 


- — w 0 + _ — [ W1 (W 1+ 2A-)) + w 2 (w 2 + 2A 2 ) ] 

A * 4 H A * 

m 22 a 22 


Thus Eq. (95) becomes upon substitution and regrouping 


H„ S - 


- (d,*, - 2c ^ fn w ;- 4c « m»;* 


22 


22 


C * 

,2 °12 


+ cn _ — (w-|W.|+ A-|W^+ w-|A^+ w 2 w 2 + A 2 w 2 + w 2 A 2 ) 


'22 
R „2 


(98) 


2c _ n [ (w 1 + A 1 )f 1 + (w 2 + A 2 )f 3 + (w 1 + A-,)^ + (w 2 + A 2 )fg ] 


Finally, recalling that Cj* = - Bj* T one can write Eq. (98) and Eqs. (96)-(97) as 


H, 


o *2 R * 

— ^ Dpj fff p Daj f p / / 

( D 11 + -=^- ) w Q + 2c li. _ll w Q - 4c « X (w Q+ A q ) 
A 22 A 22 

.2 ®21 


cn _ — (w-|W^ + A-|W^ + w-|A.| + w 2 w 2 + A 2 w 2 + w 2 A 2 ) 


'22 
R „2 


2c — n [ (w^+ A^)f-|+ (w 2 + A 2 )f 3 + (w-| + A-|)f^+ (w 2 + A 2 )f 3 ] 


(99) 


Hi 8 " 


R r m * 


2 ~ t ® 21^1 + ^ 2 ^26 ® 61^3 (® 11 “ 2 ® 66^ n ^1 2 ^ 16 n ^3 1 _ 4c -p M w -j + A -|) 

» 2„, on* -> 3 ,., i , > 4 ^ R 


R 


[ D 11 w 1 + 4D 16 nw 2 - (D 12 + 4D 66^ n w 1” 2D 26 n w 2 ] + 4c — nx(W2+ A2) 

4c .5. n 2 [ (w Q + A Q )f-|+ 2(w 1 + A 1 )f 2 + 2(w 2 + A 2 )f 4 + (w-,+ A 1 )f 2 + (w 2 + A 2 )f 4 ] 


(100) 



Next combining the last Eq. (44) with the semi-inverted constitutive equation (see Appendix A for details) 
one can write 


5 t S S S 

M x = — - D(M Xq+ M x ^cos n0 + M x ^sin nG) 


R 

C 11 N x S+ c l2 N y + c l*6 N xy + D 1*l K x + D iVy + D i6 K xy 


•C..F, 


* * 

/+ C HO F, XX - C^gF )X y- 


ii r, yy °i 2 ’xx" 


4 c 

Dh.W, 


i i vw -xx _ L 'i 2 vv -yy 


DhoW iV 


2d; 6 w, 


xy 


( 102 ) 


Substituting for W and F from Eqs. (9)-(10), regrouping and equating coefficients of like terms one gets 




R r ^ * 


o" 2 — I c n x + c i2^Pe _ V " c i6 T J " D ii w c 


(103) 


R * 


M xi 2 — ( C l2 f 1 c il n * f 1~ C l6 n y " ( D 11 w 1~ Di 2 n 2w i+ 2Di* 6 nw 2 ) 


(104) 


M x 2 2 — ( C l2 f 3 C ll n2f 3 + C l6 nf l) ” ( D 11 w 2“ D i*2 n2w 2" 2D l*6 nw 1^ 


(105) 


But now from Eq. (20) 


f”= 1 1 B21 


o = 2 R W o" -=V W ° + 7 TT -=~T [ W 1 (W 1 + 2A 1> + w 2(w 2 + 2A 2 ) ] 

A A n A 

A 22 


A 22 A 22 


thus Eq. (103) becomes upon substituting and regrouping 


M V S = 


- (D* r 2^31 ) w”- 2c fli w„- " R 


o 2 — (Cfi^- + C-^Pe - ^ 16 x ^ 


22 


1 V 

'22 


-* (106) 
C ? i2 

+ n — — [ Wi(w-|+ 2A-|) + W 2 (w 2 + 2A 2 ) ] 

d A * 
a 22 


Finally, recalling that C^= - one can write Eq. (106) and Eqs. (104)-(105) as 
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M v 


g * 2 g * 

" ( D n + -=j- ) w 0 + 2c ~r -zrp w o + 2 -r ( B ii^ + B 2 i p e" B 6i T ) 
A 22 A 22 


c ? B 21 

— n _ — [ w-|(w-(+ 2A-|) + W2(W2+ 2A2) ] 
2 * 

a 22 


M 


x-j" 2 — ( B 21 f 1~ B 1 1 n “ B 61 n V ( D 11 w 1~ B 12 n ^ w 1 + 2D 16 nw 2) 


* " — * 2 


M x S 2 = ‘ 2 y ( B 21 f 3 _ B 1*1 n 2f 3 + B 61 nf l) - ( D 1*1 w 2 _ D 1 + 2 n2w 2“ 2D l6 nw 1 ) 

Next one must express the in-plane shell displacements u s and v s in terms of the shell 
W, W and F. Eliminating initially v from the strain-displacement relations 

Ey= V,y- 1. W + ' (W,y+ 2W,y)W,y 


Yxy= u, y+ v, x+ (W, x + W, x )W, y+ W >x W,y 


yields 


u> yy = ’Yxy,y~ e y,x - -p W >x~ ( w -x + W, x )W,yy- W, x W,yy 
Introducing the semi-inverted form of the constitutive equations (see also Appendix A) one gets 


u >yy= ( A i6 F, yy + A 26 F, xx A 66 F '*y B 61 W, xx B 62 W, yy” 2B 66 W ’xy)’y 
_ ( A 1 2 F, yy + A 22 F ’XX“ A l6 F, xy“ B 21 W, xx~ B 22 W, yy -2B 26 W, xy)’X 


1 


^W, x - (W, x+ W, x )W,yy- W )X W,yy 


Substituting for W, W and F from Eqs. (8)-(10) one obtains after some regrouping 

u,yy= _Lu,00= _L{U-|(x)cosne + U2(x)cos2n9 + U3(x)sinn0 + U4(x)sin2n0} 
R 2 R 2 


where 


U^x) = 1{-[A 2 Vi - ( A iV A 66 )r,2f 1 ” 2A 26 nf 3 + A 1*6 n 3f 3] +c -^ n 2 [ < w 0 + A o> W 1 +w o A 1 ^ 
1 t ro * ... /n * on * \ „ 2... . /OD * _Q * \ _,,, . Q * „ 3„ 


+ -2r- [B 21 w 1 '( B 22 _2B 66 )n W 1 + ( 2B 26~ B 61> nw 2 +B 62 n w 2l- cv V 


(107) 

(108) 

(109) 
variables 

(110) 
(111) 

( 112 ) 

(113) 

(114) 


( 115 ) 
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~ _ -J “ )(t w •" j|t *”+ o' . // , 

^2(*) ~— { ~ [A pp fp ~( A i p + Agg)4n f2 -AA 26 n U +8A^gn (4] 

+ -|-ir |2 [( w 1 +A 1 )w 1 +w 1 A 1 -(w 2 +A 2 )w 2 -w 2 A 2 ]} 


(115b) 


^ 1 


— * — 


U 3 ( x ) » -i{ [A 22 f 3 - (A 12 + A 66 )n^f 3 + 2 A 26 nf 1 - A 1 gn 3 f 1 ]+c_Ln 2 [(w 0 +A 0 )w 2 +w 0 A 2 ] 


1 trn* 


* _OD * 2,., /on * n * 


+ 2-R [B 21 W 2 - ( B 22' 2B 66) n w 2^ 2B 26- B 6l) nw 1- B 62 n w lJ' cw 2} 


(115c) 


^ ^ J — p ^ /// _ ^ ^ ^ ^ // _ 

U4XX) =_.{ “^22^4 “(^12 + ^66^ n U + ^26 n ^2 6 n ^2^ 

+ |-^ n2 [( w 2 +A 2) Wl +w 2 A 1 + ( w 1 +a 1 )w 2 +w 1 A 2 ]} 


(11 5d) 


and where Eq. (20) has been used to obtain the following expression for f 


r-i±f? iw"- c 

0 2 R T» 0 TT 


.w 


c t n 


0 + 2B~ [(W 1 +A 1 )W 1 +W 1 A 1 + ( w 2 +A 2) w 2 +w 2 A 2] 


22 


'22 


'22 


(116) 


Integrating Eq. (114) twice with respect to 6 yields 


u =t ( U-| (x) cos n0 - _L U 2 cos 2n0 - _L U 3 (x) sin n6 - __L_ 0 4 (x) sin 2n0 + U 5 (x) -0 + U 6 (x)} 

4n 2 n 2 4n 2 

Notice that because of the periodicity condition 

2n 

J u, e d0 =0 -> U 5 (x) = 0 

o 


(117) 


(118) 


In order to determine the form of the remaining function U 6 (x) one must consider an alternate derivation 
of u. 

Using the strain-displacement equation 


£x =u, x+ ^(W, x+ 2W lX )W, x 


(119) 


and the appropriate semi-inverted constitutive equation one gets 


U .x=ex-^(W, X +2W, X )W, X 


(120a) 


U ’X =A 1* 1 F -yy +A 12 F -xx - a i* s F ’Xy W >xx -B 1 2 W, yy -2B* 6 W lXy -i.(W, x + 2W, X ) W,, 


(120b) 
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Substituting again for W, W and F from Eqs. (8) - (10) one obtains after some regrouping 

-i / + A' A' A' __ A' A' 

u, x = _u = _{U 0 (x) + U 1 (x)cosn0 + U 2 (x)cos2n0 + U 3 (x)sinn0 + U 4 (x)sin2n0} (121) 


where 


A ' m 1, /7 * i . 7 * r 7 * Tn . 1 t ,5 * A 12 7 * 


12. 


U 0 W = - { “ (A 11 X+A 12Pe- A l6 x ) + 2-R (B 21 — - B 11 )w o- c — w o 


'22 


'22 


(122a) 


c t A 12 _ 2 


c t 


t r f ft t ft 


+ '4'R^r n [(Wl +2A 1> W 1 + (w2 + 2A 2 )w 2 ]-__[2(w 0 +2A 0 )w 0 +(w., +2A i) w i +(w 2 +2A 2 )w 2 ]} 


'22 


U 1 ( x ) _ 7T^ A 12 f 1 A 1 1 n 2f 1 A 16 n V 


-it - * " — * o — * 'ct ' ' ' ' ' ' (122b) 

'■^■j=f (B 11 w 1 _B i 2 n *2B 16 nw 2 )-^[(w 1 *2A 1 )w 0 *(w 0 .2A 0 )w,]) 


t r r t r 


U 2 (x) =— { (A^f^ -4A-I*! n 2 f 2 -2A* 6 nf 4 ) -°JL[(w' +2A 1 )w 1 -(w 2 +2A 2 )w 2 ]} 

C 4 n 


(122c) 


U 3 (x) =1{ (A 1 2 f 3 - A^ n \ + A 1 * 6 nf 1 ) 


t t r f / 


'^ (B 1*1 W 2 _B 1*2 n 2w 2 - 2 B 1 6 nw 1 ) '|^ [(W 2 + 2 A 2> W 0 + < w o +2 V w 2l > 


(122d) 


t t f ft 


U4(x)=l{(A* 2 f 4 -4A* 1 n 2 f 4+ 2A 1 gnf 2 ) -£^[(w 2+ 2A 2 ) Wl +( Wl +2A 1 )w 2 ]} 


(1 22e) 


and where Eq. (20) has been used to eliminate f Q . 
Integrating Eq. (120) once with respect to x yields 


A A _ A A_ A A 

u =t {U 0 (x) + U -) (x) cos n0 + U 2 (x) cos 2n0 + U3(x) sin n0 + U 4 (x) sin 2n0 + Uy(0) } 


(123) 


Comparing the 2 expressions obtained for u (Eqs. (117) and (123)) one sees that they will be identical if 
and only if, 


U o (x)=0 6 (x) ; U 1 (x) = -_Lu 1 (x) ; U 2 (x) =-_Lu 2 (x) ; U 3 (x) = - _L0 3 (x) 

n ^ /I n ^ n ^ 


U 4 (x)=- 


4n‘ 


U 4 (x) 

(124a) 


U 7 (0) =0 


(124b) 
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Returning to Eq. (117) and substituting for u=u s from Eq. (43) yields 


s s 


t(u Q +u 1 cosne+u 2 cos2n9 + u 3 sinn0 + u 4 s sin 2 n 0 ) 


=t {^(x) _L U-| (x) cos n0 Lu 2 (x) C qs 2n0 -_L U 3 (x) sin n0 1 U 4 (x) sins2n0 } 


( 125 ) 


4n‘ 


4n‘ 


Equating coefficients of like terms one obtains 


u 0 S =U 6 (x)=U 0 (x) 


(126a) 


U 1 " -^ t J 1 (x)-_J_{A 22 f 1 -( A i*2 +A 6*6 )n2f 1 " 2A 26 nf 3 +A 1*6 n3f 3 +cw i 
cn 


n 

1 trr * 


(126b) 


-£" j = j - [ B 21 w 1 _ ( B 22 _ 2B 6 * 6 )n2w 1 + ( 2B 2 * 6 _ B 6l ) nw 2 + B 6 * 2 n3w 2l- c -^ n2 [( w 0 +A 0 ) w 1 + w o A 1 ^ 


u 2 s - ' o 2 ®. » 

4n 4cr\*- 

, ~ ” — * q c t ? ' ' ' ' ' ” (126c) 

~ 4A 26 nf 4 +8A 16 n f 4'|R n K w i +A 1 )w 1 +W 1 A 1 -(w 2+ A 2 )w 2 -w 2 A 2 ]} 


U 3 ~2 ^3( x ) = ~2 { A 22 f 3 “ (A 1 + 2 +A 66 > n % +2A 26 nf 1 ~ A 1*6 n3f 1 +cw 2 


cn 


(126d) 


"2 'r^ B 21 w 2 ~^ B 22 2B 66) n2w 2 ( 2B 26 _B 61^ nw 1 _B 6*2 n3w ll“ c 4- n2 f( w 0 +A o^ w 2 +w o A 2J} 


u 4 S = --U°4(x) = 
4n^ 


1 r7 * 


* - * 


4cn‘ 


( A 22 f 4 ( A 12 +A 66^ 4n f 4 


+4A 26 nf 2 8A 16 n ° f 2"-|-^- n2 [( w 2 +A 2) w 1 +W 2 A 1 + ( w -j +A -|) W 2 +w i A 2 ]} 


(126e) 


Thus once the solution of W and F has converged at a specified load level, one can then evaluate the 
above expressions for u. The only exception is u^ , the evaluation of which requires the solution of the 
following ordinary differential equation 


A A * 

d S 1/7** 7*— T*“ It— *^19— * "A io 

-= u 0 -~-5-( A 11 X+ A 12Pe- A i6 T ) + 2^-p( B 21 — _B 1l ) w 0 -— w o 


22 


'22 


(127a) 


1 A 

4^^T n2f(w 1 +2A l) w 1 + ( w 2 +2A 2) w 2]"^[ 2 ( w o +2A 0 ) w o +( w 1 +2 A i) w 1 + (w 2 +2A 2 )w 2 ] 


'22 


u o ( x = 0)=0 


- to suppress rigid body motion 


(127b) 
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To obtain a similar expression for v s one begins with the following strain-displacement relation 


e y =v >y -_Lw + l(W, y +2W, y )W, y 


and introduces the appropriate semi-inverted form of the constitutive equation to obtain 
v, y =e y+ .lw - i.(W, y +2W, y )W, y 

v.y = A 12 F >yy +A 22^’XX “^e^’Xy _B 21^’XX -B 22^’yy _2B 26^’xy + -j^-W--l(W, y +2 ^.y) W, y 
Substituting for W, W and F from Eqs. (8) - (10) one obtains after some regrouping 

v, y =_L v,0 =_L {Vi (x) cos n0 + V 2 (x) cos 2n0 + V3(x) sin n0 + V 4 (x) sin 2n0} 

* R R 


where 


Vi (x)=-l{(A 2 * 2 f 1 -A* 2 n 2 f 1 -A 26 nfg)-i.l.(B 21 w 1 -B 22 n 2 w 1 +2B 26 nw 2 )+cw 1 } 

V 2 (x)=_L{(A 22 f 2 -4A 12 n 2 f 2 -2A 26 nf 4 ) +^.l.n 2 [(w 1 +2A 1 )w 1 -(w 2 +2A 2 )w 2 ]} 

V 3 (x) =-l{(A 22 f 3 -A 12 n 2 f 3 +A 26 nf 1 ) -ll.(B 21 w 2 -B 22 n 2 w 2 -2B 26 nw.,)+cw 2 } 

V 4 (x)=l{(A 22 f 4 -4A 1 * 2 n 2 f 4 +2A 26 nf 2 ) +-£JLn 2 [(w 2 +2A 2 )w-| +(w-| +2A 1 )w 2 ]} 

C 4 H 


(130a) 


(130b) 


(130c) 


(130d) 


and where Eq. (20) has been used to eliminate f . 

Integrating Eq. (129) once with respect to 0 yields 

v=t{^Vi(x)sinn0+_Lv 2 (x)sin2n0-J-\/3(x)cosn0-_Lv 4 (x)cos2n0+V5(x)} (131) 

In order to determine the form of the unknown function V 5 (x) one must consider an alternate derivation of 
v. Solving for v, x from the strain-displacement relation for y xy (see Eq. (Ill)) yields 

v>x = Yxy _u, y - (V/>x + W, X ) W, y - W, x W, y (1 32) 

Introducing for y xy the appropriate semi-inverted form of the constitutive equation one gets 


v, x _A l6 F, yy +A 26 F, xx ^66^’xy B 61 ^’XX B 62 W, yy 2B 66 W ’xy -u -y _ ( w ,x+W, x )W,y-W I)< W, y (133) 
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Substituting for W, W and F from Eqs. (8) - (10) and for u from Eq. (123) one obtains after some 
re grouping 


A' 


1 t A ' - A/ — , - 

v> x = -p ^ = -p { v 0 ( x ) + V 1 (x) sin n0 + V 2 sin 2n6 + V g (x) cos n0 + V 4 (x) cos 2n0 } 


(134) 


where 


V 0 (x) -1{-(A* 6 X + A 26 p e -Agg-r) +-l-i(B 21 ^ -Bg^w' -c^w 0 


22 


22 


+ T-R-^? n2[(w 1 +2A l) w 1 + (w 2 +2A 2 )w 2 ]-|i.n[(w 1 +A 1 )w 2 +w 1 A 2 -(w 2 +A 2 )w 1 -w^]} 
A 22 


(135a) 


Vi (x) =-!{ (A 26 f3- A 1 gn 2 f 3 + A 6 * 6 nf ') 
1 t /5 * ' - 


' 2 'r (B 61 w 2 B 62 n2w 2 _2B 66 nw l) _ ^1 +c -^ n t( w 0 +A 0 ) W 1 +w 0 A lH 


c r, . _ t 


f / 


(135b) 


V 2 (x) =1{( A 26 V 4A 1 6 p2f 4 +2A 66 nf 2 ) 


+A 1 )W 1 +W 1 A 1 ~( w 2 +a 2 ) w 2- w 2 a 2]} 


(135c) 


V 3 (x) {(A 2 gf 1 A 1 g n 2 f-| - Ag 6 nf 3 ) 

4^ (§ 61 w 1 - § 62 n2w 1 +2 B 6 6 nw 2) + ^°3-cln[(w 0+ A 0 )w 2+ w 0 A 2 ]} 


(135d) 


V 4 (x) -1 { (A 26 f 2 -4A 1 * 6 n 2 f 2 -2Ag 6 nf 4 ) 


+ -^^4"|-^ n [( w 2 +A 2 ) w 1 +w 2 A 1 + ( w l +a 1 )w 2 +w'a 2 ]} 


and where Eq. (20) has been used to eliminate f 
Integrating Eq. (134) once with respect to x yields 

A — A — . A — A a _ A 

v = t(V 0 (x) + Vi (x) sin n0 + V 2 (x) sin2n0 + V 3 (x)cos n0 + V 4 (x)cos2n0 + V 6 (0)} 


(135e) 


(136) 


Comparing the 2 expressions obtained for v (Eqs. (131) and (136)) one observes that they will be identical 
if and only if, 


V 0 (x)-V 5 (x) ; V 1 (x)=lv 1 (x) ; V 2 (x) =_L.V 2 (x) ; V 3 (x) = -lv 3 (x) ; V 4 (x) =-_Lv 4 (x) (137a) 
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and 

V 6 (0)=O 

Returning to Eq. (131) and substituting for vsv s from Eq. (43) 

t (v® + v^ s sin n0 + v 2 s sin 2n0 + cos n0 + v 4 cos 2n0) 

=t { V5(x) + V-| (x) sinsn0 + _L V2M 2n0 - — ^(x) cos n0 -_L V 4 (x) cos 2n0 } 

Equating coefficients of like terms one obtains 
v q =V 5 (x) =V 0 (x) 


(137b) 


(138) 


(139a) 


V-j - — V-j(x) = — ^22^1 ”^12 n ^26 n ^3 “o‘p‘^21^1 ^22^* w-j +2B2gnw2) + cw^ } (139b) 

n cn 2 n 


V 2 ^ 2 (x) = 23iT {A 22 f 2 _4A l*2 n 2f 2 " 2A 26 nf 4 + |^ n2 [< w 1 +2A -|) W 1 -(w 2 +2A 2 ) w 2 ]} 


(139c) 


v 3 S = -J r V3( x )=-_L { X 2 V 3 / 


A^n 2 f3+A * 


' 12 ' 


26 nf 1 ^^ B 21 w 2 B 22 n2w 


l 2 -2B 26 nw 1 )+cw 2 } 


(139d) 


v 4 S = -^V 4 (x) = -_l_{A2 2 f 4 -4A 1 + 2 n 2 f 4 +2A 2 gnf 2 +£i.n 2 [(w 2 +2A 2 )w 1 + (w-| +2A-,)w 2 ]} (139e) 

Thus once the solution of W and F has converged at a specified load level, one can then evaluate the 
above expressions for v. The only exception is v®, the evaluation of which requires the solution of the 
following ordinary differential equation 


7 * - 7 * T\ . 1 t ,o * A 26 ~ * 


— V o“ -( A 16 X+A 26 p e A 66 x) + 2 ^-r (B 21^T' B 61 )w o “ 


26. 


dx 


-W r 


22 


'22 


(140a) 


+ ^^^ n2 [( Wl +2A-,)w-| +(w 2 +2A 2 )w 2 ]-ll.n[(w 1 +A 1 )w 2 +w 1 A 2 -(w 2 +A 2 )w 1 -w 2 A-,] 
A 22 


v o S (x=0)=0 


- to suppress rigid body motion 


(140b) 
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Turning now to the general elastic boundary conditions, in the limiting case as E r -> 0 one can 
regroup Eqs. (78)-(82) as follows 

h 1= a 1°l(w 0 s + W v+ W pe+ W t ) + a 12 w 0 ®- - H q s = 0 


o /... s 


h 2= a 2 V w o + W v + w p e + W t) + a 22 w o & ,x “ ( M Xo + = 0 

h 3 = a 11 u 1 S+ a 12 v 1 S+ a 13 w 1 S+ a 14 w 1^x " N xV 0 

. 1 S is 1 s 1 s s 

h 4 - a 21 u 1 + a 22 v 1 + a 23 w 1 + a 24 w 1x - N xyi = 0 


*31 U 1 + a 32 v 1 + a 33 w 1 + a 34 w Tx " H 1 S 


h 6 = a 4 1 u 1 s + a^v®* aj 3 w® + aj 4 w®- - M x * = 0 

h 7 = + a^Vg- N x ® 2 = 0 

h 8~ a 21 u 2 + a 22 v 2~ N xy 2 = 0 

u 3 s 3s 3 s 3 s ai s 
h 9 = a 1lU g - a 12 Vg + a 13 w 2 + a 14 w 2 - - N Xj} = 0 

h 10= 4l U 3 + a 22 v 3 + a 23 w 2 S+ a 24 w 2^ " N xy 3 = 0 


,3s 3s 3 s 3 s ..s 
5 31 u 3 + a 32 v 3 + a 33 W 2 + a 34 w 2,x " H 2 


. 3 s 3s 3s 3s s 

h 12= a 41 u 3 + a 42 v 3 + a 43 W 2 + a 44 w 2,x " M x 2 = 0 

h 13= a 1 4 1 u 4 + a 1 4 2 v 4 S - N x 4 = 0 

h 14 = a 21 u 4 + a 22 v 4 S ~ N xy 4 = 0 

Notice that upon replacing the shell variables involved u s , v s , w s , w, s - and N ®, N ® H s , M® by their 
equivalent expressions in terms of the variables used in the anisotropic shell analysis W, W and F 
derived earlier, one obtains a set of highly nonlinear boundary conditions. In short these equations shall be 
denoted as 


h(x, Y(x) ; l, p e , x) = 0 


An alternate set of nonlinear boundary conditions can be obtained by rewriting the general elastic 
boundary conditions in the limiting case as E r — > oo given by Eqs. (83)-(87) as follows 
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91- f,°,H 0 s . b° 2 (M x S o * Xq) - (w 0 s . W v * Wp e - W,) - 0 

92- » 2 >o * b 2>xV - "oV 0 

93= b n N Xl + b i2 N xyi + b 13 H 1 + b 14 M xi" U 1 ' 0 

94- b 2lN x S , * ♦ b^H, 8 . b^ 4 M 8 - V 8 - 0 

95- bJ,N x 8 * bggNj^ + bg 3 H 8 + b 34 M 8 - w, 8 - 0 

96- b 41 N x S , * b « N x?, * b 43 H 1 S * bX - », S X - 0 

97- b l 2 l N xV b 1 2 2 N xy 2 - u 2’ 0 

(143) 

98- »2 2 iN X 2 » b 2 2 2 N X y 2 - v 2 8 - 0 

g 9 = b 11 N X3 + b 12 N X y 3 + b 13 H 2 + b 14 M X2 u 3 = 0 

910= ^21^X3 + ^22^xy0 + b 23^2 + b 24^X2 v 3 = ® 

.3 .,s u s K1 s , 3 u s .3 , . s s _ 

9l1 = b 3i N X3 + b 32 N xy3 + b 33 H 2 + b 34 M X2 w 2 = 0 

912 = b 41 N X3 + b 42 N xy3 + b 43 H 2 + b 44^X2 w 2,x = 0 

.4 ..S . 4 .. S s . 

913= b n N x 4 + b l2 N xy 4 - u 4 =0 

914= b 21^x 4 + b 22 ^*xy 4 ~ v 4 = ® 

Notice that upon replacing the shell variables involved u s , v s , w s , w, s - and N s , N s , H s , M ^ by their 

X x xy x 

equivalent expressions in terms of the variables used in the anisotropic shell analysis W, W and F 
derived earlier, one obtains indeed an alternate set of highly nonlinear boundary conditions. In short the 
equations shall be denoted as 


g(x, Y(x) ; X, p e , t) = 0 


3. NUMERICAL ANALYSIS 

All the known numerical techniques for the solution of nonlinear equations involve iterative 
improvements of initial guesses of the solution. Working with different shell configurations^ 1 5,1 it was 
found that it is more efficient to adjust the initial guesses of the solution at a limited number of "matching 
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points rather than trying to adjust the solution at all points along the interval simultaneously, as required 
by the standard finite difference schemes. Thus, in the following, the problem will be cast into a form that 
is suited for the so-called "Parallel Shooting Method"^^ that will be used for obtaining numerical solutions. 

3.1 The Governing Equations in Terms of the Vector Variable Y 

Introducing now the 28-dimensional vector variable Y defined as follows 


Yl=fl 

y 8 = f-| 

YiS- f i 

II 

CM 


n 

Y 9= f 2 

Y 16= f 2 

/// 

Y 23” f 2 


Y 3 =f 3 

Y 10= f 3 

Y 17= f 3 

Y 24= f 3 


ii 

5? 

Yi 1 = f 4 

Yl8= U 

/// 

Y 25= f 4 

(145) 

Y 5= w 0 

r 

Y 12= w 0 

// 

y 19= w q 

/// 

Y 26 = w q 


ii 

J? 

/ 

y i 3 = Wl 

ff 

Y 20 = w 1 

y 27= w 1 


Y7= W2 

f 

y 14~ w 2 

// 

Y 2 i= w 2 

rrr 

Y 28 = w 2 


where ( )' = 

d/dx ( ), then 

the system of 

governing equations (Eqs. (21)-(27)) 

and the general elastic 


boundary conditions (Eqs. (142) and (144)) can be reduced to the following nonlinear 2-point boundary 
value problem 


4y=1(x, Y ; X, p e , T) 
dx 

for 0 < x < Jl 
R 

(146) 

g(x = 0, Y(o) ; X, p e , x) = 0 

at x = 0 

(147) 

h(x = 1 , Y (ii) ; X, p e , x) = 0 

at x = Jl 
R 

(148) 


where the general nonlinear boundary conditions are specified by the 14-dimensional vectors g and h . 
The solution of this nonlinear 2-point boundary value problem will then locate the limit point of the 
prebuckling state, whereby one of the 3 possible load parameters (X, p e , x) is chosen as the variable load 
A. The remaining two load parameters are assigned fixed values. By definition, the value of the variable 
load parameter A corresponding to the limit point will be the theoretical buckling load A s . 

As can be seen from Fig. 6, using load increments AA the solution fails to converge close to and 
beyond the limit point. This situation is somewhat unsatisfactory, especially since without the appropriate 
starting values the nonlinear iteration scheme will fail to converge also at load levels less than the 
theoretical buckling load. A closer look at the solution curve presented in Fig. 6 reveals, however, that one 
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should be able to extend the response curve around the limit point by using increments of the appropriate 
deformation parameter instead of increments in loading. 

Using the concept of generalized load-displacement relationships defined initially by Budiansky^ 2 ^, if 
the variable load A is axial compression, that is A - X, then the appropriate generalized displacement is 
the unit end-shortening defined as 

2kR L 

E = 2^RL / / {U,X " qW,xx) dXdy (149) 

o o 


Notice that the significance of e is, that the product X • e represents the decrease in potential energy of 
the applied load. 

Substituting for u, x from Eq. (120b), for W, xx from Eq. (9) and carrying out the y-integration one 
obtains after some regrouping 


-* A, 


L/R 


L/R 


8 - ■ *i> * A 1*2 p e- a,V * JZ ( ®i*i ' 5 2 *i * tit 5) f J w °' 


dx 


where 


P _ 1 t 
eor c R 


A 22 


A 2 2 o 


L/R 


- — _L _li n 2 J { w-|(w-)+ 2A^) + W 2 (W 2 + 2 A 2 ) } dx 

A22 0 

L/R 

+ I r / { (W ° + 2A o )w o + J t ( w 1 + 2A l) w l + ( w 2 + 2A 2) w 2 1 } dx 


(150) 


(151) 


Notice that the terms involving the integrals represent the nonlinear part of the end-shortening. Once the 
solution of the boundary value problem (146) has converged, it is advantageous to evaluate 5 nC by solving 
the associated initial value problems, rather than using numerical integration schemes. Here one must 
solve 


iL 8 nC = C 5 ow 0 + C 5 -|W 0 - C 52 [ wi(w i+ 2A-|) + w 2 (w 2 + 2A 2 ) ] 
dx 

/ / / j t / t / / / 

+ c 53 { w 0 ( w 0 + 2A o^ + \ [ w i( w i + 2A i) + w 2 (w 2 + 2A 2 ) 1 ) 

8 n{ (o) = 0 


(152) 

(153a) 


5 nc 



= 5 0 


(153b) 


where 8 0 is the prescribed “end-shortening”. 
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If, however, the variable load A is external pressure, that is A = p e , then following Hutchinson^ 24 } the 
appropriate generalized displacement is the average normal displacement defined as 


W 


ave - 


1 

2ttRL 


2rcR L 

J J W dxdy 

o o 


(154) 


Substituting for W from Eq. (9) and carrying out the y-integration one obtains after some regrouping 


Wave - 


Wave 

"WcT 


L/R 

= A* 2 A + A 2 * 2 p e - A 26 x + c ^ J w 0 dx 

o 


(155) 


where 

W c{ 


J. 

c 


(156) 


The term involving the integral represents the nonlinear part of the average normal displacement. Notice 
that the significance of W ave is, that the product p e • W ave represents the decrease in potential energy 


of the applied load. When the solution of the boundary value problem 
advantageous to evaluate W ave by solving the associated initial value 
numerical integration schemes. Here one must solve 

(146) has converged, it is 
problems, rather than using 

d W - cR w 
- 7 = vv ave - -j— w o 

dx L 


(157) 

w ave(°) = 6 


(158a) 

Wave(^) - W<& 


(158b) 

where is the prescribed "average normal displacement". 



If, on the other hand, the variable load A is torsion, that is A = t, then the appropriate generalized 
displacement is the apparent shearing strain^ 23 } defined as 


2xR L 

f J V ’X dxdy (159) 

0 0 

Substituting for v, x from Eq. (134) and carrying out the y-integration one obtains after some regrouping 


Y = Y - 


A* L/R 

* ~ . 7 * Tx . 1 t ~ * a 26 z * 


■,dx 


Yen ( " Al6 ^ " A2ePe+ AeeT) + 2 L (B 21 ■ = T" " B 61 } f w o dx ~ ^ -Jr f w o' 

^22 0 A 22 o 

A* L7R 

C t 26 2 C 

+ 4 -j- — n J ( w^w-jh- 2 A-|) + w 2 (w 2 + 2 A 2 ) } dx ( 160 ) 

Ag2 0 

L/R 

c t r ' ' 9 f 

~ ~2 "l n J ^ w -|( w 2 + A 2 ) + A^W 2 ~ [ w 2 (w-( + A-|) + A 2 w-( ] } 


dx 
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where 

Wt? (161) 

The terms involving the integrals represent the nonlinear part of the apparent shearing strain. Notice that 
the significance of y is, that the product x • y represents the decrease in potential energy of the applied 
load. Once the solution of the boundary value problem (146) has converged, it is advantageous to 
evaluate y by solving the associated initial value problems, rather then using numerical integration 
schemes. Here one must solve 


rl n 7/ 

— T n " c 550 w o “ C 551 w 0 + C 552 

[ w-)(w-|+ 2 A-|) + W 2 (W 2 + 2A2) ] 


dx 

f t t t 

(162) 

~ c 553 

{ w^W 2 + A2) + A^2- [ w 2 (w-|+ A-j) + A 2 wi ] } 


y ne (o) = 0 


(163a) 

— nC/ Lv — o 
y (-p-) = y 


(163b) 


where y ° is the prescribed "apparent shearing strain". 

One of these initial value problems, expressed in terms of the unified vector variable Y, represents the 
additional equation needed when solving the 2-point boundary value problem (Eq. (146)) using increments 
of the appropriate deformation parameter instead of increments in applied loading! 1 6 ] 

3.2 Numerical Analysis using Load Increments 

Due to the highly nonlinear nature of the problem anything but a numerical solution is out of question. 
Trying to make use of the readily available coded subroutines for solving nonlinear initial value problems 
lead to the decision to employ the "shooting method" to solve the nonlinear 2 -point boundary value 
problem represented by Eq. (146). Though, due to the numerical instabilities, it may become necessary to 
employ "parallel shooting" over (say) 8 intervals to carry out the integration over shell lengths greater than 
L/R = 1 . 0 , for the purpose of describing the method let us consider just "parallel shooting" over 4 intervals. 
Let us associate the following 4 initial value problems with the above 2-point boundary value problem 

-i 

— — Ug= f(x, Ug ; A., p 0 , x) 0 < x < x-| 

dx 

Ug(0) = 8 

Vp= f(x, Vp ; X, p e , x) Xi < x < X 2 Backward Integration 

dx 


Forward Integration 


(164) 


V p (x 2 ) = P 


(165) 
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■^= Up-J(x, Up ; X, p e , x) x 2 £ x < x 3 

U p (x 2 ) = P 


Forward Integration 


( 166 ) 


4 V t = f(x, V, ; X, p e , x) 
dx 


Xo < X < Jl 

d R 


Backward Integration 


(167) 


Notice that the unknown boundary conditions are specified by the 28-dimensional vectors s, p and t In 
addition, at the shell edges (at x=0 and x=L/R) the unknown vectors s and t must satisfy the general 
elastic boundary conditions, specified by the 14-dimensional vector equations g or h, derived earlier. 

Under appropriate smoothness conditions on the nonlinear vector function f(x,Y ; X, p e , ?) one is assured 
of the existence of unique solutions of these initial value problems. As can be seen from Fig. 7 these 
solutions must satisfy the following matching conditions 

U s (x = s ; X, p e , x) = V p (x = *i, p ; X, p e , x) 

- - _ ~ _ (168) 
U p (x = x 3 , p ; X, p e , x) = V t (x = x 3 , t ; X, p e , x) 


Introducing the new 56-dimensional vector function 9 these conditions can be written as 


2(§) - 


where 


U s (x = x 1( s ; X, p e , x) - V p (x = x 1( p ; X, p e , x) 
U p(x = x 3 , p ; X, p e , x) - V t (x = x 3 ,J ; X, p e , x) 


(169) 


s 


S - 



is an 84-dimensional vector 


(170) 


Finally, by introducing the 84-dimensional vector function 9 one can combine the given nonlinear 
boundary vectors g and h and the above matching conditions into a single vector equation. 


g(x = 0, Y(0) = s ; X, p e , x) 
9(S) 

= -R- - Y (i) = t ; X, p 6 , t) 


<KS) = 


(171) 
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Thus the solution of the nonlinear 2-point boundary value problem (Eqs. (146)-(148)) has been 
transformed to the solution of the 4 associated initial value problems (Eqs. (164)-(167)) and to the finding 
of the roots S of the system of simultaneous equations defined by Eq. (171). 

Using Newton’s method for finding the roots of <|>(S) = 0 one has the following iteration scheme: 

S v+1 = S v + AS V v =1 ,2 (172) 

where AS V is the solution of the 84th order linear algebraic system 
3d) 

gg- (S v ) • AS V = - $(S V ) (173) 

To apply Newton’s method one must be able to find the Jacobian matrix J(S V ) 



9gi 

9s 2 8 

^914 

9s 2 8 

9<Pl 9<Pi 

3s 2 s 9pi 


3(p 2 8 9(p 2 8 

3s2q 9pi 

3cp 2 9 

9Pl 


9<P56 

9P1 


0 


9<P1 

9P28 


5<P28 

9p 2 8 

3<p 29 

9P28 


9<P56 9<p 5 0 

9p 2 8 9*1 

ah-, 

ITT 

gh 14 

9*1 
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Notice that those components of this Jacobian, which involve derivatives of the components of the 
specified boundary vectors g and h, can be calculated analytically. However, the components involving 
derivatives of the matching conditions must be obtained by solving the appropriate variational equations. In 
order to solve for these components let us introduce the following new 28-dimensional vectors 


Wr-^Us for i = 1,2 28 


Wj= Up for i = 29,30 56 

JPj ~ H 

j = i-28 



for 


i = 29,30,...,56 
j = i-28 


Zj=^_V t for i = 57,58,.. .,84 

dtj 

j = 1-56 


(175) 


These vectors are then found as the solutions of variational equations obtained by implicit differentiation of 
the corresponding associated initial value problems. Thus the following variational equations must be 
integrated forward 


~?= W i = -?■ . 1 ( x > U s > Pe' T ) * W i 

dx oU s 

0 < x < x-( 

Wj(0) =Jj 

fori = 1,2, ...,28 

h " lu p -I (x - : p e- T> 

x 2 < x < x 3 

Wj(x 2 ) -J| 

for i = 29,30,... ,56 
j = i-28 

1 the following variational equations must be integrated backward 

•^= - __ f (x, V p ; X, p e , x) • Zj 

X-| < x < x 2 

Zj(x 2 ) =Jj 

for i = 29,30,. ..,56 
j = i-28 

(x, v t;X, p e . x) • Zj 
dx oV { ~ 

Xq < x < — 
d R 

?i (^) -Ji 

for i = 57, 58,. ...84 
j = i-56 


(176) 


(177) 


(178) 


( 179 ) 
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Here I; = {0,...,0,1,0,...0}'*" is the i® 1 generalized unit vector in the n-space. Notice that the components of 
the Jacobian matrix J' 


J'= 




at-i at! 3t-j 

au 1 au 2 9u 28 

at 2 

aUj" 


af 28 at 28 

aCTT " 


(180) 


can be calculated analytically. For details see Appendix G. 

Since the Jacobian J' is a function of U (or V), the variational equations (176)-(177) depend step-by- 
step on the results of the associated initial value problems (164) and (166), whereas the variational 
equations (178)-(179) depend step-by-step on the results of the associated initial value problems (165) 
and (167). Thus the variational equations depend on the initial guess S v . Hence it is advantageous to 
integrate per interval the 28 variational equations simultaneously with the corresponding associated initial 
value problem. This results for parallel shooting over 4 intervals in four 812-dimensional, 1st order, 
nonlinear ordinary differential equations. 

One of the greatest difficulties in "shooting" consists of obtaining a starting estimate of the initial data 
which is sufficiently close to the exact initial data, so that the iteration scheme used to find a solution of the 
nonlinear problem will converge. Fortunally, in the cases of imperfect cylindrical shells loaded by axial 
compression, external pressure or torsion the nonlinear solution approaches the linearized solution 
asymptotically as the variable load parameter A -» 0. Thus for sufficiently low values of the variable load 
parameter A one can use the linearized solutions as starting values for the nonlinear iteration scheme. 
Solutions of the linearized problem are also obtained by the "shooting method". It is well known that for the 
linearized 2-point boundary value problem Newton's method yields the correct initial value S directly 
without the need of iterations^!. The solution of the associated initial value problems and the variational 
equations is done by the library subroutine DEQ from Caltech’s Willis Booth Computing Center. DEQ uses 
the method of Runge-Kutta-Gill to compute starting values for an Adam-Moulton corrector-predictor 
scheme. The program includes an option with variable interval size and uses automatic truncation error 
control. For proper convergence (5 digits accuracy) at low variable load levels 2 iterations are sufficient; 
however, at variable load levels close to the limit point 6-12 iterations may be needed to obtain the same 
level of convergence. Thus, upon reaching a variable load level corresponding to point A in Fig. 6, instead 
of further increments in the variable load AA increments in the nonlinear part of the appropriate 
deformation parameter are used to continue the integration. As can be seen from Figure 6 this switch in 
increments makes it possible to integrate around the limit point at B and get converged solutions on the 
decreasing part of the solution curve. 
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3.3 Numerical Analysis using Increments of Deformation 

Here one must find a solution for the nonlinear 2-point boundary value problem defined by Eqs. (146)- 

(148) under the restriction that the solution must satisfy also a constraint condition specified by either Eq. 

(149) if the variable load is axial compression (A = X), or by Eq. (154) if the variable load is external 
pressure (A = p e ), or by Eq. (159) if the variable load is torsion (A =?). Thus if one increases the 
dimension of the unified vector variable Y by one (-» Y 29 = A, f 29 = F A ) then one must solve the following 
nonlinear 2-point boundary value problem 


H A A _ _ 

- 5 = Y = f (x, Y; X, p e , X) 
dx 

for 0 < x < ^: 
R 

(181) 

g(x = 0, Y(o); X, p e , x) = 0 

at x = 0 

(182) 

h( x = ^ - Y (1); X, p e , x) = o 

at x = Ji 
R 

(183) 

A = A 0 

prescribed deformation 

(184) 


Notice that the length of the capped vectors is incremented by one with respect to their uncapped 
equivalents in Eqs. (146)-(148). Thus the unified variable Y and the nonlinear vector function t are 29- 
dimensional, whereas the general nonlinear boundary conditions are specified by the 14-dimensional 
vectors g and h. 

Proceeding as before let us associate the following 4 initial value problems with the above 2-point 
boundary value problem 


d , A , A — A _ _ 

y s !(** y s . k Pe» 


U s (0) 


s 

0 


0 ^ ^ *1 Forward integration (185) 


HA A — A _ _ 

Yp = K X - Yp ^ P e - *) 



x i ^ x < x 2 Backward integration (186) 


d A A — A __ _ 

Up = f(x, Up ; X, p e , t) 


x 2 < x < x 3 


Forward integration 


y p (x 2 ) = 


P 

0 


(187) 
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H A A — A - - 

V t = f(x, V t ; X, p e , x) 
dx ~ 


A L 


t 

0 


x 3 < x < — 
J R 


Backward integration (188) 


As mentioned earlier, under appropriate smoothness conditions on the nonlinear vector function 
f(x, ^ ; X, p e , x) one is assured of the existence of unique solutions of the above initial value problems, 
which besides the usual matching conditions 


U s (x = x 1t s ; X, p e , x) = V p (x = x 1; p ; X, p e , x) 
U p (x = x 3 , p ; X, p e , x) = V t (x = x 3 , t ; X, p e , x) 


(189) 


must also satisfy the following constraint condition 
A s (xi) - A p (x-() + A p (x 3 ) - A t (x 3 ) = A 0 


(190) 


Introducing the new 57-dimensional vector function tp these conditions can be written as 


$(§) = 


U s (x = x 1( s ; X, p e , T) - Vp(x = x-i , p ; X, p e , x) 
U p (x = x 3 , p ; X, p e , x) - V t (x = x 3 , t ; X, p e , x) 
A s (xi) - Ap(x-j) + A p (x 3 ) - A t (x 3 ) - Aq 


= 0 


where 


^ = 



is an 85-dimensional vector 


(191) 


(192) 


Finally, by introducing the 85-dimensional vector function $ one can combine the given nonlinear boundary 
vectors g and h and the above matching and constraint conditions into a single equation 



g(x = 0, Y(o) = s ; X, p e , x) 

<p(f) 

h(x = ± , Y(i.) - 1 ; X, p e , x) 


(193) 
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Thus the solution of the nonlinear 2-point boundary value problem (Eqs. (1 81 )-(1 84)) has been 
transformed to the solution of the 4 associated initial value problems (Eqs. (185)-(188)) and to the finding 
of the roots § of the system of simultaneous equations defined by Eq. (193). 


Using Newton’s method for finding the roots of * 

§ v+1 =£ v + A§ v 


0 one has the following iteration scheme 


(194) 


where AS V is the solution of the 85th order linear algebraic system 
\ $ (§) ' A§ v = - $(§ v ) 

as ~ ~ 


(195) 


To apply Newton’s method one must be able to find the Jacobian matrix J(S V ) 



a§ 


4> (S v ) = 


.A 

9<j>1 

A 

4i 

4 1 

as-, 

as£ 

as 8 4 

a7T 

A 

a<t> 2 


A 

302 

A 

302 

as7 


Ss 84 

37T 


A 

3084 

C D 
CD 

A 

3084 

as 1 * 

as 8 4 

aA 

A 

3085 

4s5 

00 

CJl 

as-, 

as 84 

SA 


(196) 


Looking at Eq. (1 93) it is obvious that there are components of this Jacobian which involve derivatives of 
the specified boundary vectors g and h, which can be calculated analytically. However, the components 
involving derivatives of the matching and constraint conditions must be obtained by solving the appropriate 
variational equations. In order to solve for these components let us introduce the following new 29- 
dimensional vectors 


A 

Wj= 


U s for i = 1 ,2,. ..,28 



for 


i = 29,30,...,56 
j = i-28 


A 3 A 

"r^Up for 


i = 29,30 56 

j = i-28 



i = 57,58, 
j = i-56 


84 


( 197 ) 
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and 



II 

Q_ 

<N i 


II 

<N l 


r) A 

S~ v r 


3 A 

— v t 

3 A ~ l 


(198) 


Notice that all capped vectors have been incremented with respect to the corresponding uncapped vectors 
of Eq. (175) by one component, namely by Y 2 g= A, such that 



W T 3A 
’ asT 


(199) 


These vectors are then found as the solutions of variational equations obtained by implicit differentiation of 
the corresponding associated initial value problems. Thus the following variational equations must be 
integrated forward 


HA A — A — — A 

iL Wj = f (x, U s ; X, p e , t) • Wj 0 < x < X! 
dx 3U S ~ 

Wj(0) =fj for i = 1,2,.. .,28 


( 200 ) 


HA a A — A _ _ A ___ 

^ Wj = A_f (x, Up ; X, p e , x) -Wj x 2 < x < x 3 

dx ~ au p 

Wj(x 2 ) =fj for i = 29,30,... ,56 

j = i-28 

and the following variational equations must be integrated backward 

HA 2 a — A _ — . A _ _ _ 

— Zj = — — f (x, Vp ; X, p e , x) • Zj x-| < x < x 2 

dx ~ ay p ~ 

Z j(x 2 ) = t j for i = 29,30,... ,56 

j = i-28 


(201) 


( 202 ) 


HA 3 A — A A 

JLZj = -Vf (x, V t ; X, p e , x) • Zj 
dx 3V t 


Xo < x < — 
d R 

for i = 57,58 84 

j = i-56 


(203) 
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where Jj * {0, 0, 1, 0, ... 0} T is the i th generalized unit vector in the n-space. 

In addition, the following inhomogeneous variational equations must be integrated forward 


d A 

iw s = 

dx 


A 

3U< 


A — A _ _ 

1 (x, u s ; X, p e , t) 


A 

yy s + 


® t 

3A ~ 


W s (0)=l A 


(204) 


iw.-l 

<Jx ■ St), 


a A 


~ ( X ’ y P ’ Pe’ yy p+ f 


A _ 


Wp(x 2 ) = 0 


(205) 


and the following inhomogeneous variational equations must be integrated backward 


4 


dx 


d ?p 


Pe’ x ) ^p + 1 


Z p (x 2 ) = 0 


(206) 


? (x, Pe> *) V 4 t- f 

jw - dA ~ 


dx 


3 <>t 


A I A 

?4>-ia 


(207) 


The generalized i th unit vectors t j in the 29-dimensional space are obtained from the corresponding ones 
in the 28-dimensional space by adding a zero in the 29th row. The forcing function 3?/3A and the unit 

A 

vectors \ A for the inhomogeneous variational equations are summarized in Appendix H. Notice that the 

A' 

components of the Jacobian matrix J 


3f-| 3f-, 

3U-| 31)2 

3f 2 

307 


A' 


J = 


JL- f (x, 0 ; X, p e , t t ) = 

3U 


3f 28 

3U7 

3f A 

■507 


3f-| 3f-| 

^^28 ^ 


3f 2 e 3f 28 
3 U 2 q 3A 
3f A 3f A 
3U 2 8 3A 


(208) 
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can be calculated analytically. For details see Appendix G. 

A' A A 

Notice that since the Jacobian J is a function of U (or V) , therefore the variational equations depend 

step-by-step on the results of the associated initial value problems. Thus the variational equations depend 

on the initial guess S v . Hence the variational equations must be integrated together with the 

corresponding associated initial value problems. For parallel shooting over 4 intervals this actually involves 

the numerical integration of four 870-dimensional, 1st order, nonlinear ordinary differential equations. The 

final values of these integrations are then used to assemble the vector function <j>(S v ) and the Jacobian 

matrix — $ (S v ) given by Eqs. (193) and (196), respectively. For details please consult Appendix I. 

~ 

4. DISCUSSION 

To test the accuracy and reliability of the collapse load calculations, initially using Khot’s 
unsymmetrically laid up glass-epoxy shell^ 26 ! close approximations of the bifurcation buckling loads for 
various boundary conditions are computed. This is done by assuming a very small asymmetric 
imperfection 

W = t^sinmniicosn-^. ; % 2 ~ O'* •‘I0 -5 (209) 

L R 

where m.n are integers denoting the number of axial half-waves and the number of full waves in the 
circumferential direction, respectively. By choosing n such that it corresponds to the critical buckling mode 
of the perfect shell, the limit point of the prebuckling state for a vanishingly small initial asymmetric 
imperfection (£^ = 0-1 -ICT^) approximates closely the bifurcation point of the perfect shell calculated by 
ANILISAt 29 l 

As can be seen from the results listed in Table 1 in all cases but one (the SS-2 boundary conditions) 
the limit loads calculated by the present program COLLAPSE are, as expected, slightly lower than the 
bifurcation buckling loads of the corresponding perfect shells calculated by ANILISA. The explanation for 
the sole exception, the case of SS-2 boundary conditions, can be found by looking at the load vs 
deformation plots displayed in Figs. 10-18. 

It is known from Koiter’s Imperfection Sensitivity Theory^ ^ that if the buckling load of a structure is 
sensitive to small asymmetric imperfections, then buckling occurs at a limit point. The existence of a limit 
point implies that for the corresponding perfect structure the initial slope of the postbuckling path at the 
bifurcation point is negative (downward). As can be seen from Fig. 11b the load vs deformation (end 
shortening) curve does not have a limit point. That is, after buckling the structure continues to develop 
considerable postbuckling strength. Thus the loss of stability of the primary path does not result in 
structural collapse. The same type of behavior is displayed in Fig. 11a, where the maximum value of the 
asymmetric mode w-) max is plotted as a function of the axial load parameter p s =X s /^ c . Notice that up to 
close to the bifurcation point the asymmetric mode w-| has vanishingly small amplitudes (of the order of 
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10 7 ). After passing the bifurcation point the asymmetric mode w-| grows very fast, but as can be seen 
distinctly, with increasing axial load level. The load vs deformation path has in this case a point of 
inflection in place of a limit point. Collapse occurs here at a higher load level not shown in fig. 11. 

To illustrate this point further the results of Koiter’s Asymptotic Theory computed with ANILISA^ 2 ^ are 
shown in Table 2. Notice that for SS-2 boundary conditions both the second postbuckling coefficient b and 
®c’ initial slope of the postbuckling path are positive. This implies that the buckling load of the perfect 
structure is not sensitive to small asymmetric imperfections, a prediction that is fully confirmed by the 
present nonlinear analysis. 

Further, as can be seen in Fig. 13, for the SS-4 boundary conditions the present nonlinear 
computaional module COLLAPSE has encountered convergence difficulties shortly after passing the limit 
point. To ascertain that the solution curve has not turned back on itself, the results close to the limit point 
have been replotted using a blown-up scale. As can be seen in Fig. 14, the limit point is very sharp but 
indeed one has obtained two distinct solution branches. 

The fact that in Fig. 15b the postbuckling path turns back on itself can be attributed to the fact that in 
this case in the deep postbuckling region the response cannot be described accurately with a single 
asymmetric mode with n full waves in the circumferential direction only, as is done in COLLAPSE. Due to 
nonlinear coupling effects additional asymmetric modes with 2n circumferential full waves are triggered. 
The suppression of these additional response modes is the cause of the turn-back phenomena. 

The response mode shapes at the limit points calculated for the different boundary conditions are 
displayed in Fig. 19. Recalling Eq. (9), for each boundary condition three curves are shown, namely the 

- axisymmetric part W v + w 0 (x) 

- asymmetric part w-|(x)cosn0 + W 2 (x)sinn 0 

- compound shape at 0=0 ~ ^v + w o( x ) + w-j(x)cosn0 + W 2 (x)sinn 0 

Comparing the response modes at the limit points calculated with COLLAPSE with the corresponding 
bifurcation buckling modes calculated with ANILISA^ 2 ^ shown in Fig. 20, one sees an excellent agreement 
as far as the mode shapes go. The amplitudes differ, since the bifurcation buckling modes in Fig. 20 are 
normalized to unity, whereas the corresponding collapse modes in Fig. 19 are not. 

Finally some remarks about the normalization and re-normalization used are necessary. All external 
load parameters, X for axial load (=(cR/Et 2 )N 0 ), p e for external lateral pressure (=(cR 2 /Et 2 )p) and? 

n 

for torsion (=(cR/Et )N X y) are normalized according to the conventions introduced by the Harvard School 
under Budiansky and Hutchinsont 2 ’ 23 ' 24 l Notice that with this notation the critical buckling load of an 
axially compressed isotropic shell using membrane prebuckling (the so-called clasical buckling load) is 

1-0. 

In order to be able to make direct comparisons with the many published results and to assess the 
possible imperfection sensitivity of the buckling load of the perfect shell for both, the load scale and the 
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deformation scale a second normalization has been introduced. For a particular boundary condition the 
total end-shortening given by Eq. (150) is normalized by the value of the linear part of the end-shortening 

8 l,m =Ai> + A-|* 2 P e - A^gt (210) 

evaluated at the limit point. The re-normalized axial load p s is obtained by using the value of the 
nondimensional perfect shell buckling load using nonlinear prebuckling and calculated by ANILISA^ 2 ^ 
as a normalization factor. 

5. CONCLUSIONS 

In order to explore the capabilities of COLLAPSE, a program which can solve the nonlinear Donnell 
type imperfect anisotropic shell equations under combined axial load, internal or external prressure and 
torsion, a series of test cases are investigated. Notice that besides the standard simply supported and 
clamped boundary conditions listed in Tables 1 and 2, also elastically supported edge conditions 
represented by elastic rings of various shapes are considered. Specifically the following topics are 
examined: 

1. Imperfection Sensitivity Predictions for various boundary conditions 

- asymptotic vs nonlinear solutions 

2. Effect of elastic boundary conditions 

- symmetric cases 

- eccentric cases 

3. Interaction curves for axial compression and torsion 

- modal imperfections 

- affine imperfections 

The results and their interpretaion for the stringer stiffened shell AS-2^ 0 ^ and for Khot’s glass-epoxy 
(-40°,40 o ,0°) shell^ are presented in Part II of this report. 
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APPENDIX A: The Laminate Constitutive Equations including Smeared Stiffeners 

In deriving the constitutive equations for the shell wall made out of filamentary laminae, both 
unidirectional or woven, it is assumed that the laminae are homogeneous orthotropic layers in a plane 
stress state. The constitutive equations of the individual lamina in turn are used to derive the stress-strain 
relations for the laminate. 

Normally, as can be seen from Fig. 8, the lamina principal axes (1,2) do not coincide with the 
reference axes of the shell wall. When this occurs, the constitutive relations for each individual lamina 
must be transformed to the shell wall reference axes in order to determine the shell wall constitutive 
relations. 

Following Ashton et al^ one obtains the stress-strain relationship for the k th -lamina in terms of the 
midsurface strain- and curvature tensors of the shell wall as 


[o] k = [QJk [e] + z[Q] k [k] 


where 


[Qlk= 


Q 11 q 12 q 16 
Q 12 q 22 Q 26 
Q 16 q 26 ^66 


Jk 


t CT lk = l^x 

ay 

T xy] k 

[e] = [e x 

£y 

Yxy] T 

[k] - [k x 

Ky 

K xyl"*" 

Oi 

ii 

O 

c 4 + 

2(Q 12 + 2Q 66 )c 2 s 2 + Q 22 s 4 


q 12 = ( q 11 + q 22" 4Qg6) c 2 s 2+ Q-|2(c 4+ s 4 ) 

Q22= QhS 4 + 2(Q 12 + 2Q 6 6)c 2 s 2 + Q 2 2C 4 
^66 = ( Q 11 + q 22“ 2Q 12“ 2Q 66) c2s2+ q 66( c4+ s4 ) 
Ql6 = (Qll- q 12~ 2 Q66) c3s+ ( q 12 - q 22 + 2 Q66) cs3 
Q 26 = ( q 11 _ q 12“ 2Q 66) cs3+ ( q 12“ q 22 + 2Q 66)° 3s 


(A.1) 


(A.2) 


(A.3) 


(A.4) 



53 


and 


Q 11= y - 11 ! q 22= 

1_v 12 v 21 


e 22 

1 " v 12 v 21 


; c = cos 0|<; ; s = sin 0« 


V 21 E 1 1 v 12 e 22 


(A.5) 


n -n - 1 2 l- i'2 . n r 

° 12 - q 21 = -= = -5 . Q66 = g 12 

1 _v i2 v 21 1 _v 12 v 21 


Notice that the factor (1/2) which appears in front of the shearing strain y X y in the strain tensor has been 
incorporated into the Q matrix. 

The laminate constitutive equations are obtained by integrating the constitutive equations for the laminae, 
(Eq. A.1), over the thickness t of the laminate (see Fig. 8) yielding 


[N] = [A] [e] + [B] [k] 
[M] = [B] [e] + [D] [k] 

where 


N 

Ajj= E (Qij)k (^k _ ^k-1 

k=1 

N 1 - 2 
Bjj= E i (Qij)k (hk- h 

k=1 ^ 


N 

E 

k=1 


l) 

a 1 A- 

• A r Ef 

k-1> 

’ B U~ — p B| J 


Et^ 

•k-1> 

- 4c 2 

: D r ^ °ij 


Et J 


for i,j = 1 ,2,6 


and 


E = reference modulus 
v = reference Poisson's ratio 

N 

t = laminate thickness = E (^k - ^k-l) 

k=1 

c = y 3(1 -v 2 ) 

Inverting the first constitutive equation one gets 


[e] = [A *] [N] + [B *] [k] 


(A. 6) 


(A.7) 


(A.8) 


where 
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[A *] = [A] -1 ; [B *] = - [A]" 1 [B] 


(A. 9) 


Substituting the above equation into the second constitutive equation (A.6) yields 


[M] = [C *] [N] + [D *] [k] 


(A. 10) 


where 


[C *] = [B] [A]" 1 = - [B *] T ; [D *] = [D] - [B] [A]" 1 


[B] 


(All) 


Finally 


A r EtA ij* ; B ij = ( 2 c/t) B* ; Dj% ( 4 c 2 /Et 3 ) Dj* 


(A. 12) 


The layered composite cylindrical shell body can further be reinforced by using closely spaced rings placed 
parallel to the y-axis and closely spaced stringers placed parallel to the x-axis (see also Fig. 9). Let 


A s’ A r 
d s- d r 
E s* E r 
e s’ e r 


G s’ G r 
's- 'r 
'ts- 't r 

t 

v s- v r 


= Cross-sectional area of stringer and ring, respectively 
= Distance between stringers and rings, respectively 
= Young’s modulus of stringer and ring, respectively 

= Distance between centroid of stiffener cross-section and middle surface of the shell 
(positive inward) 

= Shear modulus of stringer and ring, respectively 

= Moment of inertia of stiffener cross-section about the middle surface of the shell 
= Torsion constant of stiffener cross-section 
= Thickness of shell 

= Poisson’s ratio of stringer and ring, respectively 


then one obtains the following additional "smeared" stiffener terms that must be added to the previously 
derived shell wall stiffness coefficients for layered composite shells 


s E c A o 

Al 1 = A i i + -J— 


A A s E r A r 

A22= A22+ 


B 11 = B ii + e 


‘r 

11 T ° s -d 


s _ E s A s 


B 22= B 22 + e r 


s . _ E rAr 
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where 


d 22= D 2 S 2 + 


Es's 

"^T 

E r 'r 

d r 


D 66“ D 66 + 


1 , G s't s 

4 d<. 


G r 't r 

“d 


G s 


E S 

2(1 +v s ) 


G r 


E r 

2(1 + v r ) 
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APPENDIX B: The Periodicity Condition 

If the solution is to satisfy periodicity in the circumferential direction then the following expression 
must hold 

2tiR 

/ v .y d y = 0 (B.1) 

o 

where from Eq. (128) 


v,y A 12 F,yy+ AopF.xx A 2fi^’XV B P1 W fXX - BppW, 


2i v, ’xx D 22 vv -yy~ ^°26 vv ’xy 


2B * W + ^ 

<:D 9fi VV 'XV + -p- 


W,y(W, y+ 2W,y) 


Substituting for W, W and F from Eqs. (8)-(10) and regrouping yields 


— * " 


■y = ^- v, 0 |"( A i2^ +A 22 p e " A 26 t) +c ( W v +w p +w t) +A 22 f o ~-£-j^ B 21 w o +cw o 


”-2.An 2 [(wi +2A^)w-j +(W 2 +2A2)W2]| 
+ ^{v 1 (x)cosn0+V2(x)cos2n0+V3(x)sinn9+V4(x)sin2n9| 


(B2) 


where the functions V 1 (x),...,V 4 (x) are given by Eqs. (130a) - (130d). 

Substituting this expression into Eq. (B.1) and carrying out the y-integration one obtains 

|(cW v -A * 2 X) +(cW p -A 22 p e ) + (cW t +A 26 t) + A 22 f 0 + cw 0 

--^-l-n 2 [(wi +2A-|)wi +(W2+2A2)W2]|(2 jiR)=0 


Recalling that from Eq. (20) 

A 22 f o"^®2*1 w o +cw o--|-^ n2 [( w 1 + 2A 1 )w 1 +(w 2 +2A 2 )w 2 ]=C 1 x + C2 (20) 

Substituting this expression into Eq. (B.3) yields after some simplification 


{(cWy A.j 2 X.) +(cWp -A 22 p e ) +(cWj + A 2 gX) +C-|X +C2} = 0 


(B.4) 


Notice that if one lets 


W v"-3- A 12 X 


(B.5a) 



(B.5b) 


W P = -5- A 22 p e 

W t = -lA 2 + 6 x (B.5c) 


C-| -C 2 =0 


(B.6) 


then the periodicity condition specified by Eq. (B.1) is identically satisfied. 
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APPENDIX C: The Coefficients Used in the Governing Equations 
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where 


A-|= D 11 A 22 +B 21 B 21 

A * - ?A * + A * 

M xxyy dfl 12 ft 66 

D * _ OD * R * 

D xxxy~ "26~°61 

B xxxy = 2 ( B 16 + *W = 4B 16 

B xxyy = B 11 +B 22 _2B 66 

B xxyy = ^12 +B 21 +4B 66 = 2 ^12 +2B 66^ 

B xyyy" 2B 16~ B 62 

B xyyy = 2 ( D 26 +B 62^ 4B 26 

and 


q= 4c - 5 . q 

t 2 
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APPENDIX D: Derivation of the Reduced Boundary Conditions 

It is necessary to express the different boundary conditions in terms of the shell variables 
W, W and F. Considering successively 


1. Boundary condition: u = u 0 

From the strain-displacement relations 


u,y= Yxy- V, x - W, x W,y- W, x W,y- W, x W,y 


“■y" V ^ w 7- ™'V Wi y 


(D.1) 

(D.2) 


Differentiating the first equation with respect to y and the second with respect to x and then eliminating the 
term v >X y between the resulting expressions yields 


u >yy = Yxy.y - 8 y,x -p W, x W, x W,yy W, x W,yy- W, x W,yy 


Using the semi-inverted constitutive equations 


(D.3) 


Yxy = A i6 F -yy + A 26 F ’ XX A 66 F,x y B 61^’ xx B 62^’yy 2B 66 WlX y 


£y- A i2 F, yy + A 22 F ’ xx A 26 F>x y B 2i^ ,xx B 22 w -yy 2B 66^’ x y 


(D.4) 

(D.5) 


Recalling further the fact that if a function <(>(x,y) in an orthogonal reference frame satisfies the condition 


<|>(x,y) = C at x = x c 


(D.6) 


where both C and x 0 are constants then 


dy r 


<|>(x,y) =0 at x = x 0 for r = 1,2,... 


(D.7) 


hence since u = u 0 at x = 0,L, therefore u,yy= 0 at the shell edges. Thus specializing Eq. (D.3) to the 
shell edges, one obtains upon substitution for Yxy and e y and some regrouping 


"j — * — * — * “ * — * 

— { - A22p.XXX + 2A 2 0F, XX y- (A^ 4 ^ Agg)F >X yy + A-|gF,yyy } 


t ,5* 


+ 2c ^ B 2i^’ xxx+ ^ 2B 26 B 6i )w>xx y + ^ B 22 _ 2B 66^’ x yy _ B 62^’yyy ^ 


(D.8) 


R 


W, X - (W, x + W, x )W,yy- W, x W,yy= 0 


,yy- 
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Substituting now for W, W and F from Eqs. (8)-(10) and regrouping yields 


* A 22 f o J -ft B 21 w o + cw o~ J -ft n 2 C W -| (w-f + A-|) + W-, A 1 + w 2 (w 2 + A 2 ) + w 2 A 2 ] } 

+ { A 22 f^ - (A^ 2 + Agg)n^f.|- 2A 2 gn^fg+ A^gn^fg+ cw^- c _L n 2 [w Q (w-(+ A-|) + w-)A Q ] 


1 t * 


- 2 t B 21 w 1 - ( B 22- 2B eV n2w 1 + ( 2B 26" B 61 )nw 2 + B 62 n3w 2] > cos 

+ { A 22 f 2 - (A 12 + Agg)4n f 2 - 4A 2g nf 4 + 8A 16 n 3 f 4 

C t p ' ' ' ' 

~ "2 "R n t w -|( w 1 + A 1> + w 1 A i~ w 2 (w 2+ A 2 ) - w 2 A 2 ] } cos 2n0 


n0 


— ^ _ 


* . 7* \„2/ . „7 * 


+ { A 22 f 3 ' (A 12 + A 66 )n % + 2A 26 nf 1" A 16 n3f 1 + cw 2” c n2 t w 0 ( w 2 + A 2> + w 2 A Q ] 


1 t r ~ * 


j -ft [B 21 w 2 ' (B 22" 2B 66 )n2w 2' (2B 26“ B eV nw 1 " B 62 n3w lJ > sin 


n0 


. rrr 


( ^22^4 ^12 + ^26^2” ®^16 n ^2 

c t 2 ' * ' f 

~ ~2 -p n t w i (w 2 + A 2 ) + w-)A 2 + w 2 (w-|+ A-)) + w 2 A^] } sin 2n0 = 0 


(D.9) 


This expression must hold for all values of 9 = y/R, which implies that the expressions in the curly 
brackets must be equal to zero. This yields then 5 conditions that u = u 0 at x = 0,L/R implies. Notice that 
the first of these equations does actually not represent a new (boundary) condition, since it can be 
obtained from Eq. (20) by a single differentiation with respect to x = x/R. 


2. Boundary condition: v = v 0 

From the strain-displacement relations and using the semi-inverted constitutive equations one gets 


v -y = £ y + - j w, y (w, y+ 2 W, y ) 


(D.10) 


A 12 F >yy + a ?9 F >xx A 9« f >xv b o-|W, 


♦ 

■ , 22 r, xx" 


26 F >xy B 21 W ’XX B 22 W, yy“ 2B 26 W, *y + ~ 4 W,y(W,y+ 2W -y) 


Specializing to the shell edges at x = 0,L/R implies v, y = 0. 

Hence upon substituting for W, W and F from Eqs. (8)-(10) one obtains after some regrouping 
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{ ( c Wy A-| 2 ^) + (cWp e A 22 Pe) + (cWj+ A 2 gi) 

+ *22* o~ J -ft ®21 w o + cw o- n 2 [w-,(w 1+ 2A^ + w 2 (w 2 + 2A 2 )] } 

— £ 7/ — * Q * -j | ^ /7 — +p — + f 

+ { Ag 2 ^ - A-j 2 n ' 1 ~ ^26 n ^3 + cw 1 _ "p ^21 w 1” ® 22 n w-j+ 2 B 2 gnw 2 ) } cos n0 
+ { A 2 2 f 2 ~ 4A 1 * 2 n 2 f2“ 2A 2 * 6 nf 4 + -£ -1. n 2 [w-)(w-(+ 2A-|) - w 2 (w 2 + 2A 2 )] } cos 2nG 

+Q * 1 t "™ + ^ “+0 “ j|c 7 

+ { A^n 13+ A 2 gnf^+ cw 2 - _ _ (B 2 ^W 2 " B 22 n^w 2 - 2B 26 nw 1 ) } sin n0 


(D.1 1) 


+ { A 22 f 4 - 


4A 


12 


n 2 f 4 H 


2A, 


26 


nf 2 + -£ -1. n 2 [w 2 (w-|+ 2A-|) + w-|(w 2 + 2A 2 )] } sin 2n0 = 0 


This expression must hold for all values of 0 = y/R, which implies that the expressions in the curly 
brackets must be equal to zero. This yields then 5 conditions that v = v 0 at x = 0,L/R implies. Notice that 
the first of these equations does actually not represent a new (boundary) condition, since the underscored 
part is just Eq. (20). Thus it vanishes identically. The remaining terms are also identically equal to zero if 
one introduces for W v , Wp e> and Wf their values given by Eqs. (B.5). 


3. Boundary condition: W = 0 

Substituting for W from Eq. (9) yields 


W v + Wp e + W t + w 0 + w-jcos n0 + w 2 sin n0 = 0 (D.1 2) 

For this expression to be true for all values of 9 = y/R, the following boundary conditions must hold at 
x = O.L/R: 

w G = - (W v + Wp e + W t ) (D.1 3a) 

w 1= w 2=0 (D.1 3b) 

4. Boundary condition: W, x - 0 

Substituting for W, x from Eq. (9) yields 

/ / / 

w Q + w.|COS n0 + w 2 sin n0 = 0 (D.1 4) 

For this expression to be true for all values of 0 = y/R, the following conditions must hold at x = 0,L/R 
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w 0 = w rw 2 =o (D15) 

5. Boundary condition N x = - N 0 

Substituting for F from Eq. (10) one obtains 

2 

N x = F,yy= __ (- X - n 2 f-|Cos n0 - 4n 2 f 2 cos 2n0 - n 2 f 3 sin n0 - 4n 2 f 4 sin 2n0) = - N 0 (D.16) 

For this expression to be true for all values of 0 = y/R, the following boundary conditions must hold at 
x = 0,L/R 


f 1 = f 2 = f 3 = f 4~ 0 

Notice that there are only 4 new boundary conditions, since the condition given by Eq. (D.17) is just the 
definition of the nondimensional axial load parameter X. 

6. Boundary condition: N xy = S 0 

Substituting for F from Eq. (10) one obtains 

Et^ _ ' / / , 

N xy = “ F -xy = ~ (" x - nf 1 sin n0 - 2nf 2 sin 2n0 + nfgcos n0 + 2nf 4 cos 2n0) = S 0 (D.19) 

For this expression to be true for all values of 0 = y/R, the following boundary conditions must hold at 
x = O.L/R 


(D.17) 

(D.18) 


cR 


Et 2 

(D.20) 

f 3 = f 4 = 0 

(D.21) 


Notice that there are, once again, only 4 new boundary conditions, since the condition given by Eq. (D.20) 
is just the definition of the nondimensional torque parameter x. 


7. Boundary condition: H = 0 

From the derivation of the nonlinear equilibrium equations using the stationary potential energy 
criterion^ 30 ], by definition 


H = M x,x + (M X y + M yx)-y + N x(W.x + W >x ) + N xy (W, y+ W, y ) = 0 


(D.22) 
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Using the appropriate semi-inverted constitutive equations one obtains 

H B ll F, xyy _ B 21 F, xxx + B 61 F, xxy~ B n^’xxx _ B i2^’xyy _ 2B 16^’Xxy 

+ 2 (~ B i6 F, yyy” B 26 F, xxy + B 66 F, xyy“ B i6 w, xxy~ D 26 w, yyy” 2D 66 w, xyy) (d.23) 

+ F,yy(W, x + W,^) - F,xy(W,y+ W,y) - 0 

Substituting for W, W and F from Eqs. (8)-(10) and regrouping one obtains the following residual 

— * 1 t — * ' ' o ’ ’ ' ' ' ' 

e H“ ( - B 2l f o" 2 -R D 11 w o“ 2a < w o + A o ) " cn t f l( w i + + f 3 (w 2 + A 2 ) + f^w-i-e A-j)+ fg(w 2 + A 2 )] } 

+ { — — [- D 11 w 1 + (D 12 + 4D 66 )n^w 1 - 4D 16 nw 2 + 2D 26 n J w 2 ] - 2cA,(w 1 + A^ + 2cnx(w 2 + A 2 ) 

- < B lV 2B 66> n2, 1* (®6*r 2 ®26> n V 2B 1*6" 3| 3 

- 2cn^[f-|(w 0 + A 0 ) + 2f 2 (w^ + A.,) + 2f 4 (w 2 + A 2 ) + f 2 (w-|+ A-)) + f 4 (w 2 + A 2 )] } cos n0 


— if ,,f — if ~~ if O f — * — A- — g q 

{ - B 2 ^f 2 + 4(B^- 2Bgg)n f 2 + 2(Bg^- 2B 2 g)nf^+ 16B^gn f 4 


if / / / 


cn [8f 2 (w Q + A Q ) + f"| (w^ + A^- f 3 (w 2 + A 2 ) - f^w-) + A-|) + fg(w 2 + A 2 )] } cos 2n0 


/ / / 


{ - cn [4f 2 (w.|+ A^) - 4f 4 (w 2 + A 2 ) - 2f 2 (w-(+ A-|) + 2f 4 (w 2 + A 2 )] } cos 3n0 


(D.24) 


It ~ — £ — £0 /— '* n “ * Q / * — ' 

+ { g- -p [- D n w 2 + (D 12+ 4Dgg)n^w 2 + 4D 16 nw 1 - 2D 26 n’®w 1 ] - 2cX(vj 2 + A 2 ) - 2cnx(w-( + A-|) 


. w _ 


* ~ * 


- B 21<3 * < B lV 2B 66>"V < B 61- 2B 26) n, r 2B 16 nJ| 1 


2cn 2 [f 3 (w 0 + A Q ) + 2f 4 (w 1 + - 2f 2 (w 2 + A 2 ) + f 4 (w 1 + A-|) - fg(w 2 + A 2 )] } sin n0 


{ - B 21 f 4 + 4(B; r 2Bgg)n 2 f 4 - 2(Bg* r 2B 2 g)nf 2 - ISB^n 3 ^ 


t / / 


cn [8f 4 (w Q + A Q ) + f 3 (w 1 + A.,) + f-|(w 2 + A 2 ) - f 3 (w-|+ A-|) - f 1 (w 2 + A 2 )] } sin 2n0 


/ / / 


{-cn [4f 2 (w 1 + A^ - 4f 4 (w 2 + A 2 ) - 2f 2 (w-|+ A-|) + 2f 4 (w 2 + A 2 )]}cos3n0 


/ / / 


{ - cn [4f 4 (w 1 + A.,) + 4f 2 (w 2 + A 2 ) - 2f 4 (w-|+ A-|) - 2f 2 (w 2 + A 2 )] } sin 3n0 = 0 


Here the 0 = y/R dependence is eliminated by using the following 3 Galerkin integrals 


J end0 = 0 


(D.25) 
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2k 

J ercos n0 d0 = 0 
o 

2n 

J CRSin n0 d0 = 0 
o 


(D.26) 


(D.27) 


Notice that these are the same Galerkin integrals that have been used earlier to eliminate the y- 
dependence from the out-of-plane equilibrium equation (3). Substituting for e H and carrying out the 
integrals involved one obtains the following boundary conditions at x = 0,L/R 


— * i t — * ’ ' 

B 21 f o + 2 -R D 11 w o + 2c Mw 0 + A 0 ) 


/ / 


ft i 


cn [f 1 (w 1 + A.,) ♦ f 3 ((w 2 + A 2 ) + f^w-j + A-j) + f 3 (w 2 + A 2 >] = 0 


(D.28) 


IB* OD * \r,2i /D * 


/ / 


B 21 f 1 ( B 11“ 2B 66) n<if 1" ( B 61 - 2B 26^ nf 3“ 2B 16 n3f 3 + 2c M w i + A.,) - 2cnx(w 2 + A 2 ) 

J ^ f°1>1 - (°12 + 4 D 6 * 6 ) n2w i + 4D* 6 nw 2 '- 25 26 n 3 w 2 ] 


t r 


f r 


r f t 


(D.29) 


+ 2cn £ [f 1 (w 0+ A q ) + 2f 2 (w 1+ A^ + 2f 4 (w 2 + A 2 ) + f 2 (w 1 + A-|) + f 4 (w 2 + A 2 )] - 0 


* * 


B 21 f 3 “ ( B lV 2B 6>V (B 61" 2B 26 )nf 1 + 2B 1*6 n3, 1 + 2c Mw 2 + A 2 ) + 2cnr(w-j + A-|) 


/ / 


1 t 


* \ 2... jn * 


+ 2 -R- [D 11 W 2 “ < D 12 + 4D 66^ n ^ w 2~ 4D 16 nw l + 2 D 26 n 3 w-,] 


(D.30) 


/ / 


ft t 


+ 2cn [f 3 (w 0 + A 0 ) + 2f 4 (w 1+ A.,) - 2f 2 (w 2 + A 2 ) + f 4 (w 1 + A-|) - f 2 (w 2 + A 2 )] = 0 


But now, differentiating Eq. (20) once with respect to x = x/R one obtains 

i t B 21 c ' c t n 2 ' ' ' ' 

Vj r -=T w o~ -=r^~ w o + j r- — t w i( w 1 + A l) + WfA 1+ w 2 (w 2+ A 2 ) + w 2 A 2 ] 

A 22 A 22 A 22 

Thus upon substituting and regrouping the first condition, Eq. (D.28), becomes 


(D.31) 


1 < m.y B Jf?i 

2 R 11 7* 

A 22 


B 


21 


t t 


) w o ~ c ~=~— w o + 2c Mw 0+ A 0 ) 


'22 


B. 


/ / 


c t 21 2 ' 

+ "2 "R ~ = T' n f w 1^ w 1 + A 1^ + w 1 A i+ w 2 (w 2+ A 2 ) + w 2 A 2 ] 

A 22 


(D.32) 


f / 


/ / f 


+ cn [fi(w 1+ A.,) + f 3 (w 2 + A 2 ) + f^(w-| + A-,) + fg(w 2 + A 2 )] = 0 
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8. Boundary condition: M x = - N 0 q 

Using the appropriate semi-inverted constitutive equation one gets 

^x = “ B 1 1 F> yy _ B 21 F, xx + B 6i F -xy- i'^'XX _ *- ) i2^’yy _ 2D 16 W, xy = “ N 0 q (D.33) 

Substituting for W and F from Eqs. (9)-(10) and regrouping one obtains the following residual 


e M~ { B ii^ B 2iPe + B 61 t+ B 21 f o + \ D 11 w o~ J ^ ^ ^ 


— x rr — 


d + — -if ,r ” +0 — * r 

+ { B 2 -\U~ B 11 n f 1" B 61 nf 3 + _ _ (D^w^ D 12 n w-|+ 2D 16 nw 2 ) } cos n0 


— s|t f/ — ^ p ' 

+ { B 21 f 2 - 48^0^2- 2B g1 nf 4 } cos 2n0 


(D.34) 


+ { B 21 f 3 - B^n^ Bg* 1 nf 1 + .1 1. (D 1 * 1 w 2 - D^* 2 n 2 w 2 - 2D 1 * 6 nw 1 ) } sin n0 


+ { B 21 f 4 - 48^0^4+ 2Bg.jnf 2 } sin 2n0 = 0 


where q = 4c(R/t' ; )q and q is the load eccentricity. 

Here the 0 = y/R dependence is eliminated by using the previously defined (see Eqs. (D.25)-(D.27)) 3 
Galerkin integrals. Substituting for and carrying out the integrals involved yields the following boundary 
conditions at x = 0,L/R 


* — — * — 


B 21^o + ~2 Tf ^H W o ( B 1 1 + B 2lPe + B 61 t = 0 


(D.35) 


* O — * ' It— * 

^ d . ' L /r\ 


* „ 2 l( 


B 21 f 1~ B 11 n f 1" B 61 nf 3 + 2 -R (D 11 W 1" D 12 n w 1 + 2D 16 nw 2> = 0 


(D.36) 


* 2x , d * ^ t /r\ * 


B 21 f 3~ B 11 n ^ f 3 + B 61 nf 1 + ~2 ~k (D 11 W 2" D l2 n2w 2“ 2D 16 nw 1 ) = 0 


(D.37) 


Substituting for f Q from Eq. (20) and regrouping the first condition, Eq. (D.35), becomes 


— * — * 

It,:;* b 21 b 21 


B. 


_ _ (D 
2 R 


* — ~ * — 


11 


) w o" c w o" ( B iV \ <0*- ~ B 2iPe + B 61 t 


22 


22 


g 

+ fl _L _£1 n 2 [w-|(w-|+ 2A-|) + w 2 (w 2 + 2A 2 )] = 0 

a 22 


(D.38) 
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The 8 different edge conditions discussed previously can be combined to yield besides the 4 different 

simply supported and the 4 different clamped boundary conditions, the free edge conditions and, if the 

initial imperfections are chosen accordingly the symmetry conditions at x = _1 Ji 

2 R 

SS-1 Boundary Condition 


N x = - N 0 

— > 

O 

II 

II 

CO 

II 

II 

T“ 

(D.39) 

Nxy = S 0 


f 1 = f 2 = f 3~ f 4“ 0 

(D.40) 

:> 

it 

o 


w 0 = - (W v+ W pe+ W t ) 

(D.41a) 



II 

$ 

ro 

ii 

o 

(D.41b) 

M x = - NoQ 


// 

W Q = C 41 X + C 2 20^ 

(D.42a) 



rr rr / 

w 1 = - B-|f 1 - B 2 w 2 

(D.42b) 



rr rr / 

w 2 = " B 1 f 3 + B 2 w i 

(D.42c) 

where 




Bl<2 
* 6,*, 

7 

n 1 r 2R ,7 * n * 7 * n * s 

c 41- — { — ( a 22 b i 1 A 12 B 21 ) 

+ A 22 q } 

(D.43) 

B 2 = ^ n 
D 1*1 

1 

~ 1 2R-T+Q* 7 * “ * 

c 22°- — ( A 26 B 21 A 22 B 6l) 

and Ai= ^5^+ B 21 B 21 



(D.44a) 

(D.44b) 

(D.44c) 


(D.44d) 
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Nxy - S 0 


f 1- f 2 = f 3 - f 4 = 0 


W = 0 


w 0 = - (W v+ W pe+ W t ) 


W 1 “ w 2 = 0 


M x = - NqCI 


w = C 41 A- + C 220 x 


where 


w^= - B-|f.| + B3f-|- B4W2 

n // / 

w 2 = * B 1 f 3 + B 3 f 3 + b 4 w i 


B, 


Bg^Illn 2 

d t r * 


D 


11 


□ 1 t B 21 

B 8- T 7 -FT 


2 R 7 * 
A 22 


2D 1r 

64= n = Bp 

— * 

d h 


Bq. 1 1 " 

y on 


2 R 7* 
A 22 


B 5 - t 2A2* 6 - <2i 2 * 6 - B 6 *,) ] 


22 
.3 _ 


D 


11 


B, 


n v “ * 1 — * — * 

B 6 = ~~r~ l A 16 - "T~ ( 2B 26- B 61> 1 


22 


D 


11 


1 1 t 9 “ * “ * “■ * — * 

b 7- ( c * 4 ^ n 2 [ B 22 - 2B 26 - (2B 26 - B 61 ) ] } 

A 22 °11 


B 10= 2n _ 


26 


'22 

* 


B H= n 


3 A 16 


22 

* 


Bi 2 . n 2 

°11 


SS-3 Boundary Condition 


N x = - N 0 

— > 

f 1= f 2 = f 3 = f 4~ 0 

O 

II 

> 

-> 

tr t r 

f 1 = B 15 f 3 + b 16 w 2 


f 2 = 2Bi 7 f 4 


(D.45) 

(D.46a) 

(D.46b) 

(D.47a) 

(D.47b) 

(D.47c) 


(D.48) 


(D.49) 

(D.50a) 


(D.50b) 
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rf , r 




f 3 = " B 15 f i " b 16 w i 

(D.50c) 



f 4 = - 2B 17 f 2 

(D.50d) 

:> 

II 

o 


w 0 = - (W v+ W pe+ W t ) 

(D.51a) 



w-j - W2= 0 

(D.51b) 

M x = - N 0 q 


// 

w 0 = C 41^ + c 220 t 

(D.52a) 



" / / 

w 1 = B 18 f 3 - B 1 9 w 2 

(D.52b) 



// / / 

w 2 = " B 1 8 f i + b 19 w i 

(D.52c) 


where 


b 15" ( A 26 D 11 + B 21 B 61^ 


Bl8 2 T Xj" ^ A 22 B 61 A 26 B 2l) 


Bi6 ~ TT (B 26 D 11” B 21°16) 


B ig- 2 -L (A 22 D 16+ B 21 B 26 ) 


b 17= n 


'26 


'22 


A 1~ A 22 D 11 + B 21 B 21 


SS-4 Boundary Condition 


U - U r 


V = V r 




rrr r 


U = 


B 28 f i ~ B 29 f 3 + b 8 w i + b 30 w i + 2B 9 a 1W ( 


^2 _ ^ b 31^2 _ 8B 32^4 + B 9 A 1 w -j - BgA2W 2 

rfr r frt/f 

f 3 = B 28 f 3 + B 29 f 1 + b 8 w 2 + B 30 w 2 + 2B 9 A 2 W 0 

nr , , , 

f 4 = 4 ^31 f 4 + 8B 32 f 2 + b 9 a 2 w i + B 9 A 1 W 2 

” t r 

f 1 = B 20 f 1 + B 15 f 3 + b 16 w 2 
f 2 = 4B 2 4f2+ 2B 17 f 4 


(D-53) 


(D.54a) 

(D.54b) 

(D.54c) 

(D.54d) 

(D.55a) 


(D.55b) 
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rr r r 

f 3 = B 2 o^3“ B 15 f 1 - B i6 w i 

(D.55C) 



f 4 = 4B 24 f 4" 2B 17 f 2 

(D.55d) 

1 

II 

o 

— > 

w 0 ~ - (W v+ W pe+ W t ) 

(D.56a) 



w-j - W2= 0 

(D.56b) 

M x = - N 0 q 

— > 

// 

w 0 = C41X, + C 2 20 t 

(D.57a) 



rr / / 

w 1= B 25 f 1 + B-ieV B 19 w 2 

(D.57b) 



rr j / 

w 2 = B 25 f 3 _ B i gf-j + B igw 1 

(D.57c) 


where 

b 2 o= (A-| 2 D 11 + B^B^) 


B 23 = 4c — — App 
" t A-| 22 


b 24= n 


A * 

2 12 


22 


B 25- 2 y -jj- ( A 22 B 11 A 12 B 21^ 


2 - 

n - t * t * 


B 28- — ( A 12 + A 66 — 
A 22 


2A 26 ,7 * ^ * . S * n * x . (2B 26" B 6l) ,7 * o 


( a 26 d ii+ B 21 B 61 ) 


Al 


( a 26 B 21 A 22 B 61^ ^ 


— * — * — * — * — * (^®pc ^f>i) — * — * — * — * 

B 29~ -= j ~ I A 16" (A 12 D 11 + B 11 B 21 ) " (A 22 B 11“ A 12 B 21 ) } 

A 22 


(D.58) 


b 30 - ^ { ^ u 2 [ - 2 b; 6 ) + 
A 22 


2A, 


A 1 


26 < B 21 5 lV B 26 5 lV 


(2B 


26 


B 


Al 


«i) — * — * 

— ( a 22 D 16 + 


B 21 B 26^ 1 


c } 
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4 : + 

2 pA A 

o n ,7 * T * ^ M 26 26 , 

B 3 1_ ^7- ( A 12 + A 66 — — ) 

A 22 A 22 


b 32 = jl i ( A ; 6 - _ i 2 * 26 ) 


'22 


'22 


A 1 A 22 D 11 + B 21 B 21 


C-1 Boundary Condition 


N x = - N 0 


N X y= S Q 


f T f 2 = f 3 = f 4 = 0 


f 1 = f 2 = f 3~ V 0 


W = 0 


w 0 = - (W v + W Pq+ W t ) 


W, x =0 


w-| = w 2 = 0 


r / / 

w Q = w 1 = w 2 = 0 


C-2 Boundary Condition 


f/t rt 


U = Ur 


f i = B 10 f 3“ B 11 f 3 + b 8 w i + b 13 w 2 

f/t rr 

f 2 = 2B 10 f 4~ 8B 11 f 4 


^xy = B o 


f 3 - " B 10 f i + B 11 f 1 + b 8 w 2 " B 13 w i 

rrr ,/ 

f 4 = - 2 Biq f 2 + 8B-| -jfp 


r f r t 


-> f r f 2 =f 3 =f 4 =0 


W = 0 


w 0 = - (W v+ W pe+ W t ) 


w-| = w 2 = 0 


(D.59) 

(D.60) 

(D.61a) 

(D.61b) 

(D.62) 


(D.63a) 

(D.63b) 

(D.63c) 

(D.63d) 

(D.64) 

(D.65a) 

(D.65b) 



73 


r r r 


W, x =0 


o 

II 

CM 

II 

<T 

ii 

o 

£ 

(D. 66 ) 

C-3 Boundary Condition 



N x = - N 0 

— > 

fl= f 2 = f 3 = f4= 0 

(D.67) 

o 

> 

II 

> 

-4 

// r rr 

f 1 = B 17 f 3 + b 8 w i 

(D. 68 a) 



f 2 ~ 2B 17 f 4 

(D. 68 b) 



f 3 = " B 17 f i + b 8 w 2 

(D. 68 c) 



Ur - 2 B 17 { 2 

(D. 68 d) 

w = o 

— > 

w 0 ~ - (W v + Wp e + W t ) 

(D.69a) 



w-| = w 2 = 0 

(D.69b) 

W,x=0 

— > 

r r r 

w 0 = W 2 ~ 0 

(D.70) 

C-4 Boundary Condition 



C 

II 

c 

o 

— > 

/// / ;// // 

f 1 = B 31 f.,- B 32 f 3 + B 8 w 1 + B 14 w 2 

(D.71a) 



w / 

f 2 = 4B 31 f 2' 8B 32 f 4 

(D.71b) 



/// r irr rr 

f 3 = B 31 f 3 + B 32 f 1 + B 8 w 2 ' B 1 4 W ! 

(D.71c) 



rrr r 

f 4 = 4 B 31 f 4 + 8B 32 f 2 

(D.71d) 

o 

> 

II 

> 

— > 

rr r rr 

f 1 = B 24 .fi + B 17 f 3 + b 8 w i 

(D.72a) 



f 2 = 4B 24 f 2 + 2B-| 7 f 4 

(D.72b) 



rr r rr 

f 3 = B 24 f 3 - B 17 f 1 + B 8 w 2 

(D.72c) 



V 4B 24 f 4“ 2B 17 f 2 


W = 0 


w, 


o= - (W v + Wp e+ W t ) 


W, x =0 


w-| - W 2 = 0 


r r t 


-* w Q = w 1 = w 2 = 0 


(D.72d) 

(D.73a) 

(D.73b) 

(D.74) 


Symmetry Condition at x = — — 

2 R 


u = 0 


Nxy= 0 


H = 0 


W, x =0 


rrr tr 


f 1 = 


B 33 f 3 + B 34 f 3 + b 35 w 2 + ( B 36 + B 37 t ) w 2 + b 37 a 2 t 


fg - 2B-|of4~ ® B 11^4 

rn n „ 

f 3 = " B 33 f i - B 34 f 1“ b 35 w i" ( b 36 + b 37^) w 1 “ B 37 A-,t 

w n 

f 4 = - 2B-)of 2 + BB 1lf2 


/ / / / 


f 1 = f 2~ *3 = f 4~ 0 


w o =0 


W 1 “ “ B 38 f 3 + B 39 f 3" b 40 w 2 + ( B 41 + B 2 3x)w 2 + B 23 A 2 i 

rr „ 

w 2 = B 38 f i “ B 39 f 1 + b 4 qW 1 - (B 4 i + B 23 7)w-| - B 23 A 1 x 


r t / 

-> w Q = w 1 = w 2 = 0 


(D.75a) 

(D.75b) 

(D.75c) 

(D.75d) 

(D.76) 

(D.77a) 

(D.77b) 

(D.77c) 

(D.78) 


where 
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B35 ' J Tf XT ‘ (S5 2V 5 61> 5 1*1- 45 21 5 1*6 1 


Q 

1 t n — * — * — *— * 

B 36 = 2 R ( B 62 D 11 + 2 B 21 D 26 ) 


b 37 = 2cn 


B 38 - 2 y 177- t 2B 26 B 61 + f 2 A 26 D 11 B 21^ 2B 26 B 61 ^ J ^ 
1 D* A 1 


(D. 79 ) 


B 39- 2§ 16 - ^ ( 2B 12 B lV A 16 5 1 + 1> 1 

D 11 


b 40 . JL ( 46^* [ (2i 2 V b 6 >,V 4B 2 *,5; 6 ] 1 

D 11 


B 41 - ^ [ 2D 2 * 6 - (B 2 * 6 D,V 2B 2 > 2 * 6 ) ] 

n a 1 


Free Edge Condition 


N x = - N 0 


fl= f 2 = f 3 = f 4 = 0 


(D.80) 


f r r r 


Nxy = B o 


V f 2“ f 3 = f 4 = 0 


(D-81) 


H = 0 -> w Q = (B42- B43 ^)w q - B43 A q X.- B44 [ w^w-] + A-|)+ w-jA^ w 2 (W2+ A2) + v^Ag 1 ( D - 82a ) 

rrr rrr rr r r ^ 

W-j = - B-| f ^ + 645^+ (B40- B47 ?i)w^ - B47A-JX + (643+ B4gi)W2+ B49A2T (D.82b) 

/// rrr rr r r 

w 2 = - B-jfg - 645^ + (643- B47 ^)w 2 - E^A^ - (B43+ B4gt)w-|- B4gA-|X (D.82c) 


^X = - NqQ -* 


w = B 42 w 0 + B 50 X + B 51 p e - B 52 x - B53 [ W!(W 1+ 2A-|) w 2 (w 2 + 2A 2 ) ] 


(D.83a) 


/^= B-(2 w 1 _ B2W 2 ~ B-|f 1 


v 26 = B 12 w 2 + B 2 w r B^g 


(D.83c) 
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where 


B 42 = 2c 1 

t A-| 


b 49= 4c v R 


1 B 1*1 


b 43 = 4c -5. _L 

t A-, 


B 5 q= 2 — — 22 ( B^ + 1 _L q ) 
SU t Ai 1 1 2 R M ’ 


B 44 = cn 


2^1 

Ai 


B 5 1=2^ 


R A 22 - * 


t A! 


B. 


21 


(D.84) 


B 45=2 JL ( Bg r 2B 2 * 6+ 

D 


11 


$ ♦ 
4B 21 D 16 


'11 


) B 52 = 2 " 


R A 22 “ 


t Al 


B ; 


61 


n 2 — * — * 

b 46 - A- ( d 12 * 4D 66 - — IL15 
D 11 D 11 


n ** 

R c n 2 B 21 

B 53= T n -t— 
2 A-| 


b 47= 4c .5. 1 


' 5 1 + 1 


R _ 2n 3 , 2D 12 D 16 x 

B 48 — ( D 26 - — — ) 

D 11 D 11 


Symmetry or Anti-symmetry condition at x = ii 

2 R 

It has been shown in Ref. [26] that if a layered composite shell is layed-up in such a fashion that in 
the stiffness matrix the coupling terms between the bending and shear strain of the middle surface are not 
zero, then the buckling and the postbuckling modes are skewed with respect to the shell axis. To 
represent this skewed pattern it is necessary to include in the circumferential Fourier representation both 
cosnQ and sinn0 terms (see also Eqs. (8) - (10)). It has been shown by Booton f 27 ^ that if the axial 
dependence of the cosnO terms is symmetric with respect to the mid-section (at x = L/2) of the shell, then 
the sinnG terms are anti-symmetric with respect to the same location or vice-versa. This implies that, if 
one uses the appropriate initial imperfection shapes (say, A 0 and A-) are symmetric and A 2 is anti- 
symmetric with respect to x = L/2) then one can work also for the general collapse analysis with half of the 
shell length only by enforcing appropriate symmetry and/or anti-symmetry conditions at x = L/2 ( or 
x = (1/2) L/R). 

In order to derive these conditions one must recall that symmetry at x = L/2 implies that 
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u = N X y= H = w, x = 0 

whereas anti-symmetry at x = L/2 requires that 


N x~ - N o- V = W = 0, M x = -N 0 q 
at the shell mid-section. 

Thus if A 0 and A-| are symmetric at x = L/2 then, considering Eqs. (8) - (10) and Eqs. (34) - (44), 
one obtains 


u-| = u 2 = 0 



H-| = H 2 = 0 


f / 



whereas the fact that A 2 is anti-symmetric at x = L/2 implies 


Nx 3 = N X4 = 0 


v 3 = v 4 = 0 


W2= 0 
M X2 =0 


Next one must express these conditions in terms of the shell variables W, W and F . 

First boundary condition 


, rrr 


^ 22*1 ^ 12 + ^ 66^ n ^1 2 ^ 26 n * 3 + ^ 16 n ^ 3 + cw i 


It * % o' * 4; ff O 

~2R l B 21 W 1 " (B 22' 2B 66 )n ^1 + (2B 26 _ B 61> nw 2 + B 62 n w 2 


t „2 


c -r n 1 w o (w 1 + A l) + W 1 A 0 3 = 0 


(D.85a) 

] 


U-| = 0 -> 
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u 2 = 0 


A 22 f 2 " (A 12 + A 66 )4n f 2" 4A 26 nf 4 + 8A 1* 6 " 3f 4 

" j n2 [ w i ( W 1 + Ai) + w-|A 1 - w 2 (w 2 + A 2 ) - w 2 A 2 ] = 0 


Nx 3 =° 


fq- 0 


Nx 4 =0 


U= o 


Second boundary condition 


Nxy^O 


— > 


f r° 


N xy 2 = 0 


f 2 = 0 


v 3 - 0 -> A 22 f 3 - A 12 n 2 f 3 + A 26 nf 1 + cw 2 - 1 JL (B 21 w 2 ~ B 22 n 2 w 2 - 2B 26 nw 1 ) = 0 

"* " "* «2,.. f \ c t ? 


v 4 0 — > ^22^4 *^26^2 + t ’d' ^ [W 2(w-| + 2A*|) + w-j(W2 + 2Ap)] = 0 

4 '*“ 4 R *” ' 

Third boundary condition 


H 1 =0 


1 t ~ * - * - * 


fff / 


-2 -p ( A 22 D 11 + B 21 B 21 ^ w o cB 21 w o + 2CA 22 Mw 0 + A q ) 

+ -| -1 B 21 n 2 [w 1 (wi+ A i) + w 1 A 1 + w 2 (w 2 + A 2 ) + w 2 A 2 ] 


* „2 


/ / 


+ cA 22 n Cf 1 (w 1 + A,) + f 3 (w 2 + a 2 ) + A-j) + fg(w 2 + A 2 )] 


h 2 = 0 


f r 


B 21 f 1 - ( B 11" 2B 66> n f 1' ( B 61~ 2B 26^ nf 3~ 2B 16 n f 3 + 2cX( Wl + A 1 ) -2cnt(w 2 


1 t r r * 


-L [D 1 * 1 w 1 - (D* 2+ 4Dg 6 )n 2 w 1 + 4D 1 + gnw 2 - 2D 26 n 3 w 2 ] 


f f 


/ / 


+ 2cn [f-,(w Q + A q ) + 2f 2 (w 1+ A^ + 2f 4 (w 2+ A 2 ) + f 2 (w 1+ A-,) + f 4 (w 2 + 


(D.85b) 

(D.85c) 

(D.85d) 

(D.86a) 

(D.86b) 

(D.87a) 

(D.87b) 

(D.88a) 

= 0 
+ a 2 ) 

(D.88b) 
A 2 )] = 0 


W 2 = 0 


(D.88c) 
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Fourth boundary condition 

r 

w Q = 0 (D.89a) 

r 

w 1 = 0 (D.89b) 

m x 2 = 0 -> B 21 f 2 - B-|* 1 n 2 f 3 + Bg 1 nf 1 + 1 _L (D^Wg- D 12 n 2 w 2 - 2D 1 * 6 nw 1 ) = 0 (D.89c) 

/ -r 

Recalling further that A 0 , A-j symmetric implies A Q = A 1 = 0 at x = 172, whereas A 2 antisymmetric implies 


A2= 0 at x = L72. 


— * — * " it — * — * — * " 

A 22 f 1 ' 2A 26 nf 3~ 2" "R ^ B 21 W 1 + ( 2B 26~ B 61^ nw 2^ = 0 

(D.85a) 

— * — * " 
A 22^2 ~~ ^ A 26 n ^4~ ^ 

(D.85b) 

— * " i t — * ” 

A 22 f 3~ J R - B 21 w 2 = 0 

(D.87a) 

A 22^4~ ® 

(D.87b) 

1 t i — * — * — * — * 

2 R -* ( A 22 D h + B 21 B 21^ w o” 0 

A 22 

— * — * — * "If-*"' “ * " 

B 2i f i ~ ( B 61 _ 2B 26^ nf 3 + "2 r - ( D 11 w 1 + 4D 16 nw 2^ = 0 

(D.88a) 

(D.88b) 

- * " i t - * " 

B 21 f 3 + J "R D 1 1 W 2 = 0 

(D.89b) 


From Eqs. (D.87b), (D.88a) and (D.85b) it follows that f^= w Q = f 2 = 0 . 

Considering Eqs. (D.87a) and (D.89b) it can be seen by inspection that the solution is w 2 = fg= 0 . 
Taking all these results into consideration the remaining 2 equations become 


- * i t - 

A 22 f 1 ~ J "r B 21 w 1 = 0 

- * 1 t “ * 

B 21 f 1 + 2 R D 11 W 1 = ° 


(D.85a) 

(D.88b) 


Once again, by inspection, their solution is f 1 = w 1 = 0 . 

Summarizing, if A 0 , A-| are symmetric and A 2 is anti-symmetric with respect to x = L/2, then the 
appropriate symmetry and/or anti-symmetry conditions at x=(1/2)L/R are 


/ / / t rr rr // /// rrr nr rrt 

f 3= f 4” w 2 = f i = f 2 = w 0 = W 1 = f 3 = f 4 = w 2 = f 1 = f 2 = w o ’ W 1 = 0 


(D.90) 
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APPENDIX E Derivation of Cohen’s Ring Equations^] 

When ring stiffeners are attached to the shell edges, boundary conditions for the shell equations must 
be generated which represent the ring behavior. Ideally the ring reactions enter the shell at a single 
meridional station, the ring boundary. A set of suitable ring equations based on moderate rotations are 
derived in this appendix by using the stationary potential energy criterion. These equations are based on 
the following assumptions 

all geometrical and mechanical properties are axisymmetric; 

the ring material is homogeneous and isotropic; 

the effects of non-uniform warping of ring sections, transverse shear strains and shear center 

eccentricity related to the section centroid are neglected. 

The origin of the ring cross-sectional xyz-axis is assumed to be at the centroid of the ring cross- 
section. The sign convention used is indicated in Fig. 2. 


For a one-dimensional theory of rings the centroidal hoop strain eq is the only extensional strain of 
consequence. The strain-displacement relations have been derived in Ref. [21] as 


ee= 1 (v'- w) + .1 (P* + p2) 


Px= 1 (v + vf) ; P z = - -1 


K x= - P 


x 72 


(v'+ w") 


Kz= i p z + i p y s ¥ ( " ir + 


and the twist per unit length of the curved bar is given by 


(E.1) 


=i P ; + iu' 
a p y a 2 


(E.2) 


where () = _(). 
de 

If one includes the initial stress-free imperfection w(0) of the ring then e 0 and p x must be modified as 
follows 


se- -1 <v'* W) * .1 <£* D=) - i. 


1 A/ 1 


P x = P x + — w=_(v+w+w) 


(E.3) 


Neglecting the effect of nonuniform torsion for a homogeneous, one dimensional ring theory the 
constitutive equations are 
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N e = EAee 

^X = ^lx K X" ^*XZ K Z 


^z _ ^I xz k x + EI z k z 


M|= GJcp 


(E.4) 


By definition the potential energy of the ring is 

n r = U r - W r 


where 


271 

U r = .1 J (Nq£ 0 + M x k x + M z k z + M^(p)ad0 
o 


(E.5) 


(E.6a) 


271 _ „ _ _ _ _ 

W r = J (N x u + Sv + Pw + M x p x + MyPy+ M z p z )ad0 (E.6b) 

o 

and 

N x , S, P, M x , My, M z are the external line loads applied to the ring. See Fig. 2 for the notation and 
sign convention used. 

The variational statement of equilibrium for a homogeneous ring under the specified external line 
loads is 

2 7t 

8n r = 5 J F(u, u , u , v, v , w, w , w , Py, p y )ad0 = 0 (E.7) 

o 

where 

F = 1 (EA£q+ EI x k^+ EI z k^- 2EI xz k x k z + GJ(p 2 ) 

— _ — — — — (t.oj 

- (N X U + Sv + Pw + M X P X + MyPy+ M z p z ) 


thus in this case the equilibrium condition 


,r an' 


an r 


an r 


8ir = _ 5u + _ 5v + 

3u 3v 3w 


5w + iE_ 8 p v = 0 

ap y y 


(E.9) 


implies the 4 conditions 
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an r = an r = an r = an r = 

3u 3v 3w 9Py (E. 1 0) 

Carrying out the indicated variations, integrating by parts whenever necessary and using the fact that the 
variations 5u, 5v, 5w and 5Py must be periodic one obtains the expressions 


an r 

3u 


2n 


m r I 3F _ d_ 5F + _d^ 
J Hu dG 3u' --9 


de 


3F 


2 3u" 


} 5u adG = 0 


an 


271 


f r 9F d 3F ls . A _ 
5v adG = 0 

3v J dv dG 3v 
o 


2k 


(E-1 1) 


an 1 _ f . 3F d aF d" 3F , s .Q 
ST ■ j { 7ST • Hff * — ■=» > 5w ade ■ 0 


3w dG aw 7 d0 2 3w" 


2rc 


an 1 _ r 3F _ d 3F . SR 

a^ " J 1 a^ dG ap7 } 5|i y 


adG * 0 


For this to be true for arbitrary nonzero variations 5u, 5v, 5w, and 5p y one must satisfy the following Euler 


equations 



3F 

d 

3F 

3u 

dG 


3F 

_ d 

3F 

3v 

~ dG 

dv' 

3F 

_ d 

3F 

15w 

dG 

aw' 

3F 

d 

3F 

3py 

dG 

aPy 


+ iiiL = 0 


2 3u" 


dG 


+ = 0 
dG 2 


(E.12) 


Carrying out the indicated differentiations, substituting and regrouping yields the following ring equations 

El z u lv - GJu + El xz v + El xz w lv - a(El z + GJ)P y = a 4 N x + a 3 (M z *)' (E.13a) 

E| XZ U + (El x + a^EA)v + El x w - a^EAw - aEl xz P y = - a 4 S *- a 3 M x (E.13b) 

El xz u lv + El x v - a 2 EAv + El x w lv + a 2 EAw - aEl xz p = a 4 P*- a 3 (M*)’ 

y * 


(E.13c) 
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(El z + GJ)u + El xz v + El xz w + aGJ(5y- aEl z (3 y = - a 3 M y (E.13d) 

Notice that in these equations the nonlinear terms have been considered as effective additional loads and 
moments applied to be linearized ring equations. These effective additional loads and moments are 


p*=p + $ -1(^-) 2 ] 

a d c. a 

s*-s+M[ ’ $ x + p z ) - 1 (i.) 2 ]' 

a 2 x z 2 a 


M x = M x - 


M*= M z - N 0 P Z 


(E.14) 
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APPENDIX F: The Boundary Stiffness and Boundary Flexibility Matrices 

Substituting in the separated form of Cohen's ring equations, Eqs. (28a)-(28d), for the ring 
displacements and the ring forces their equivalents in terms of the shell variables from Eq. (41) and Eqs. 
(39) or (40), respectively one obtains the following equations at the lower edge (at x = 0) 

El z (u s - e z w ®), 00 0 e- GJ(u s - e z w ®) I0e 

+ El xz{(-^- vS_ e x u 'y” e z w >y)>000 + ( w S+ e X W ’x)>0000} (F.la) 

- a(El z + GJ)(- w,®), 00 = a 3 { (R - q) N 0+ RN®} 


El xz( u S " e z w - x ).000 + ( EI X + a 2 EA)(^. v s - e x u,®- e z w ,®), 00 

+ El x (w s + e x w , x ), 900 - a 2 EA(w s + e x w ,®), 0 

- aEl xz (- w,®), 0 = - a 3 { - (R - q t )S 0 + RN® y ) 


(F -1b) 


E| xz (u s - e z w, x ), eee0 + El x (i. v s - e x u,®- e z w ,®), 


-yA 00 0 

r S , 


a2EA ( -p V S e X U -y e Z W 'y)«0 + El x( W S+ e X W 'x)>0000 

+ a 2 EA(w s + e x w,®) - aEl xz (- w,®), 00 = a 3 {RH s } 


(Fl.c) 


(El z + GJ)(u s - e z w,®), 00 + EI XZ (JL v s - e x u ,®- e z w ,®), 0 


+ El xz (w s + e x w, x ), 00 + aGJ(- w,®) >00 - aEl z (- w,®) = 

-a 2 {- e z RN®+ (q - e z )(R - q)N 0 + e x RH s + RM S } 


(F.ld) 


Substituting for the shell variables their Fourier representations from Eqs. (43) and ( 44 ), respectively, 
yields after some regrouping the following expressions 

1111 o p 

( a n u l + a i 2 v 1 + a l 3 w i + a 14 w i,x " N x-|) cos ne + ( a i i u 2 + a i 2 v 2 " N X 2 )cos2n0 

,3 3 3 3 44 (F.2a) 

+ (a 1103 + a 12 v 3 + a 13 w 2 + a 14 w 2 ~ - N X3 )sin n9 + (a 11 u 4 + a 12 v 4 - N X4 )sin 2n0 = 0 

( a 2 i u i + a 22 v-|+ a 23 w-| + a 24 w 1x - N xyi )sin n 0 + (a|.,U 2 + a^^-N^Jsin 2 n 0 

, 3 3 3 3 4 4 (F.2b) 

+ (a 21 u 3 + a 22 v 3 + a 23 w 2 + a 24 w 2i x ~ N xy 3 ) cos ne + ( a 21 u 4 + a 22 v 4 - N X y 4 )cos 2n0 = 0 
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e H“ { a i°i(W v + W pe+ W t + w 0 ) + a 1 ° 2 W Q - - H 0 } 

+ (a^u-)* a32 v 1 + a 33 w 1 + a 34 w -( x ” H l) cos n0 + ( -)cos 2n0 + (,..)cos3n0 
3 3 3 3 

+ (83^3+ ag2 v 3+ a 33 w 2 + a 34 w 2x ” *"*2) s ' n n0 + (•••) s ' n 2n0 + (,..)cos3n9 


e M~ ( a 2i( w v + W pe+ W t + w 0 ) + a 2 °w o - - [M X()+ Xq(1 - * 


)] } 


1 


1 


1 


1 


( a 41 U 1 + a 42 v 1 + a 43 w 1 + a 44 w i,x " M x-,) cos n0 + (-)cos2n0 
( a 41 u 3 + a 42 v 3 + a 43 w 2 + a 44 w 2,x “ M x 2 ) sin n0 + <-) sin 2n0 


where 


a 1 ° 1 =<Y )2 Dll 


a 21 ~ 


- <T |2 5l2 


a 1°2' <”> 2 d 12 


0 

a 22 = 


- r p e x - 

- [d 22 + {II) 2 JL d 12 ] 


a ii= — (-) 2 n 2 D 33 
11 4c a 


a lV 4~ (~) 2 n4 °35 
IJ 4c a 


a io = - — (— )* n° — D34 
12 4c v a 7 d 44 


1 / 1 \2 „3 a p'. 
R 


a iV i 4> 2 " 2 688 


a lV 7 4> 2 n 2 °77 


4 - - ! 4 )2 " 3 w 578 


a?i= — (-) 2 n 2 D33 
11 4c a 


a lV i 4> 2 " 4 °35 


3 _ ^ / t ^2 n 3 3 n 

12- 4J <a' n R°34 


a «- i 4 )2 n2 688 


a lV 7 4> 2 " 2 °77 


a 12" 1 4> 2 n3 W 5?s 


1 _ 1 /t,2„3n 


21 


^ W n D 43 


4c a 


a 23 = - 4s 4 )2 n ° 45 


1 1 , t ,2 „ 2 a n 

a 22 _ T - ( — ) n -pr D44 

44 4c a R 


a 24' 4? 4 )2 " 048 


(F -2c) 


(F -2d) 


(F 3) 


(F .4) 


(F.5) 


(F -6) 


(F.7) 


(F.8) 
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2 _ 2 A 2 3 n 

a 21 n D 87 


4=l(l) 2 n 2 a 6 88 

(F.9) 

a 3 _ 1 / 1,2 _3 p; 

a 21 4 c af n ° 43 


4= i 4 )2 n °45 


a 3 _ 1 #t »2 n 2 a = 

a 22 4 c "a n ^ ° 44 


4= - 4 4 2 - d 46 

(F . 10 ) 

a 2 4 1 = - | 4 ) 2 n3 5 87 


4=l(l ) 2n2 a 5 88 

(F.11) 

4 = n 2 £) 2 6 53 


a 33’ <-J )2 ° 55 


a 32 = " n & D 54 


4= 4 )2 ^ 5e 

(F.12) 

4 = n2 <Y )2 5 53 


4= A 2 555 


a 32~ n <4 ^54 

a 


4= A 2 656 

(F.13) 

_ 1 _ n 2 R /n . e z n ® x 

a 41" n Y ( ° 63 T ° 33 " T 

D 53 ) 

a 43“ ( n2 D65 + " 4 4 ° 33 ' 4 D 55 ) 

a a 3 


"1 R 3 "” 967 ^ 6y “ 

a 42 = n -j (" -r D 64“ n -g- °34 + -p- d 54) 

1 R 0 67 6v — 

a 44 = 7 T (° 66 + n — d 36” — d 56) 
c* a a 

(F .14) 

4 = n2 “ (°63 + £ ^ 33 - 

D 53 ) 

a 43" IT ^ n2 ^ 65+ p4 it ^ 35 " 4 


a 42 = n Y ^ ^ 64+ H 2 4 ^ 34 

6v — 

- t? °54> 

a 44 = ■j ^66 + H 2 4 ^ 36 - -£■ D 56 ) 

(F.15) 
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d 22 

^33 

D 34 

^35 

^36 

^43 

644 

D 45 

^46= 
°53= 
d 54= 
d 55 = 
d 56= 
d 63~ 
d 64= 
d 65 = 
^66 
D 75 

^76 


EI Z 

aD 


O EU Gy 

n d [ t - _ ( 

aD R aD 


x Z + GJ 
aD 


El 


xz 


aD 


El 


xz. 


(1 - — ) ( 

R aD 


!£ - (± + n 2 f£) (fli) -(-*—) (— ) 

R aD v R R ; aD ^ V R r' W 


El 


_ (. 

R aD 


x + a 2 EA 


El 


) -(-£) 
aD 


aD 


El 

■aD 


x + a EA 


aD 

e. 


El, 




2 e z, , a^EA 
R M aD 


(i. ♦ n 2 f£) ( El 
R R aD 


xz. 


[n 2 & + 1?!*] 
R aD aD 


.2 r El xz 


e x E'x, 


a 2 EA, 


n [___ - (_JL)] - _n ( 

aD R aD R aD 


n 2 (fli) . if!* 
aD aD 


El, 


a 2 EA, 


n (1 - _±) (^1) + (1 - n c —) ( 

R aD v R ’ aD 


6v A Ely 

— [n 4 ( _) + 

R v aD 


a 2 EA 


.2/ a 


] - n (_ + n‘ 
aD ' R 


El 


xz, 


R ) ( aD 


El- 


GJ 


El 


xz, 


i- + WU A / 

aD aD "R" aD 


El 


xz 


aD 


e 7 El 

(1 - — ) ( 

R aD 


xz, 


El 


xz 


— ( 

R v aD 


El, 


) - (“ + n* If) (Z±) - n 2 (1 + _±) (^i) 
R R aD R R aD 


, GJ, 


n 2 [Ili£ 

aD 


e X Ely, 


a 2 EA, 


16 -1 (_)] - 4 _ ( 

R aD' R aD 


4n 2 (f!i) + 

aD aD 


(F.16) 
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P El z 

D 77 -4n^_2 


e x El 

_ ( 

R aD 


“>j * ^ 

aD 


D 78= 


El 


xz 


aD 


6 79 =^ 

/y aD 


“ ©v Ek 
°87= t? <-£ 


R 

Ek 


aD 


a 2 EA 


aD 


D 88 = + £?!* 

88 aD aD 


El 


) - (-£) 

aD 


and the load parameters are 


^ = N 0 and q = 4c il q (F.17) 

Et 2 t 2 

Equations (F.2a) and (F.2b) must hold for all values of 0 = y/R, which implies that the expressions in the 
brackets must be equal to zero. This yields then the following 8 boundary conditions at x = 0 


a 11 u 1 + a 12 v 1 + a 13 w 1 + a 14 w i,x ' N x-| = 0 
a fl u 2 + a i2 v 2 - n x 2 = 0 
a 1 3 1 u 3 + a T2 V 3 + a f3 w 2 + a U w 2,x ' N ><3 = 0 
and a 1 4 1 u 4 + a lV4- n x 4 = 0 

a 21 u 1 + a 22 V 1 + a 23 W 1 + a 24 w 1 ,x “ N xyr 0 
a 21 u 2 + a 22 v 2“ N xy 2 = 0 
a 21 u 3 + a 22 V 3 + a 23 w 2 + a 24 w 2,x ' N *y 3 = 0 
a 2l u 4 + a 22 v 4~ ^xy4 = 0 


(F .18) 


(F.1 9) 


For equations (F.2c) and (F.2d) the 0 = y/R dependence is eliminated by using the following 3 Galerkin 
integrals 
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2n 

J EjdG = 0 
0 

271 

J £j cos nQdQ = 0 (F.20) 

o 

2n 

J £j sin n0d9 = 0 
o 

where for i = 1, £j= e-| = £(.) and for i = 2, £j= £2= e^. Notice that the above are the same Galerkin 
integrals that have been used earlier to eliminate the y-displacement from the out-of-plane equilibrium 
equation (3). Substituting, in turn, for eh and e^ and carrying out the integrals involved one obtains the 
following 6 boundary conditions at x = 0 

a 1°1<W v + w p e + W t+ w 0 ) + a° 2 w o - - H 0 = 0 

®31 U 1 + a 32 v 1 + a 33 w 1 + a 34 w l ,x ' H 1 = 0 (F.21) 

a 31 U 3 + a 32 v 3 + a 33 w 2 + a 34 W 2,x ‘ H 2= 0 

a 2 °|( w v+ W Pe+ W t + w 0 ) ♦ a 22 w o - - (M Xq + Xq) = 0 

a 41 u 1 + a 42 v 1 + a 43 w 1 + a 44 w 1 ,x " M Xl = 0 ( F - 22 ) 

a 41 u 3 + a 42 v 3 + a 43 w 2 + a 44 W 2,x ' M x 2 = 0 

Notice that in the first of equations (F.22) q/R is neglected as compared to one. This same approximation 
has been used earlier in equations (46). 

The boundary conditions at the upper edge (at x = L) have the same form as the ones at the lower 
edge (at x = 0), except that all boundary stiffness coefficients a^ at x = L involve a minus sign in their 
definition as per Eqs. ((F.3) - (F.15)). Notice that these minus signs are due to the definition of the forces 
and moments acting at the ring centroid, which according to Eqs. (39) and (40) involve a minus sign, 
depending upon wether the ring is located at x = 0 or at x = L. 

Finally, it must be mentioned that the elastic boundary conditions given by Eqs. (F.18), (F.19), (F.21) 
and (F.22) are valid in the limit as E r -> 0 (free edges). 


Another set of elastic boundary conditions valid in the limit as E r -4 00 (fully clamped edges) can be 
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derived by inverting the boundary stiffness matrix to yield the boundary flexibility matrix. To facilitate this 
operation initially the elastic boundary condition obtained earlier are regrouped into the following matrix 
notation 


“ 

" 


"1 


1 

0 

a 11 

0 

a 12 

< 

W v + Wp e+ W t+ w 0 

1 = 1 

I 

► * 

o 

X 

a ° 

a 21 

a ° 

22 


2= 

O 

XI 


M + Xq 
x o J 


(F.23) 


1 1 1 1 


' 


r 

a 11 a 12 a 13 a 14 


U 1 


^x-| 

1111 





a 21 a 22 a 23 a 24 


V 1 


N xyi 

J J 4 J 



r - ’ 

> 

I I 1 1 

a 31 a 32 a 33 a 34 


W 1 


H 1 

1111 





a 41 a 42 a 43 a 44 


. W 1,x . 


, M xi . 


(F .24) 


2 

a 11 

1 

CVJ ^ 
cfl 


u 2 


n x 2 ' 

2 

a 21 

a 2 
22 j 


' v 2 - 


.NxyJ 


(F.25) 


3 

3 

3 

3 

l 11 

a 12 

a 13 

a 14 

3 

3 

3 

3 

l 21 

a 22 

a 23 

a 24 

3 

3 

3 

3 

l 31 

a 32 

a 33 

a 34 

3 

3 

3 

3 

l 41 

a 42 

a 43 

a 44 


u 3 


\ 

< 

w 

II 

%3. 

w 2 


h 2 

xl 


, m x 2 , 


(F.26) 


4 

a 11 

r 


u 4 


N x 4 ' 

4 

a 21 

4 

a 22 


i 

v 4 J 


^xy 4 . 


(F .27) 


Inversion of these (2x2) and (4x4) boundary stiffness matrices yields the following (alternate) elastic 
boundary conditions valid in the limit as E r -» °° (fully clamped edges) 



91 


b 1 ° 1 H 0 + b 1 ° 2 (M Xo+ Xq) = W v+ W pg+ W t+ w 0 
b 2l*"*o + b 22^x 0 + = w 0 x 

bliNxi* b l2 N xy-| + *- ) l3^1 + b i4M x ^= Ui 
b 2i^Xi + b 22^xy-| + b 23^1 + ^ 4 ^x 1 = V 1 
b 31 N x-| + b 32 N xy-| + b 33 H 1 + b 34 M x-) = W 1 
b 41 N x-( + b 42 N xy-| + b 43 H 1 + b 44 M x-| = w -| x 

2 2 

b 11 N X2 + b 12 fv| xy2 = u 2 
b 2l^X2 + b 22^xy 2 = v 2 

b 1 3 1 N x 3 + b l 3 2 N xy 3 + b i 3 3 H 2 + b l 3 4 M x 2 - u 3 
b 21 N x 3 + b 22 N xy 3 + b 23 H 2 + b £ 4 M X 2 - v 3 
b 31 N x 3 + b 32 N xy 3 + b 33 H 2 + b 34 M x 2 = w 2 
b 41^x 3 + b 42 N xy 3 + b 43 H 2 + b 4 3 4 M x 2 = w 2,x 


(F.28) 


(F .29) 


(F -30) 


(F.31) 


b 11 N X4 + b l2 N xy 4 - u 4 
4 4 

b 21 N X4 + b 22 N xy 4 = v 4 


(F -32) 


Notice that once again the boundary conditions at the upper edge (at x = L) have the same form as the 
ones at the lower edge (x = 0), except that all boundary flexibility coefficients b^ at x = L involve a minus 
sign in their definition as per Eqs. ((F.3) - (F.15)). 
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APPENDIX G: Components of Jacobians J (Eq. (180)) and J (Eq. (208)) 

Using the 28-dimensional vector variable Y defined on p. 31 the nonlinear governing differential 
equations (Eqs. (21)-(27)) can be put into the form given by Eq. (146), where the 28-dimensional vector 
function f has the following components 


*1= Y 8 

*8 = y 15 

*15= y 22 

*2= y 9 

*9 = y 16 

*16= y 23 

*3= Y 10 

*10= y 17 

OJ 

>■ 

II 

r- 

*♦— 

*4” Y 11 

*11= y 18 

LO 

II 

00 

»*— 

*5= Y 12 

*12= y 19 

*19= y 26 

*6= Y 13 

*13= y 20 

*20= y 27 

*7= Y 14 

*14= y 21 

*21= y 28 


(G.1) 


and 


*22“ C 1 y 15“ c 2 y 1“ c 3 Y 20" c 4^( y 20 + a -|) + C 5 Y 6 + C 201 Y 24 + C 202 Y 10 + C 203 Y 28 

+ C 204 y 14 + c 21lPe( Y 6 + A l) + C 212 T ( Y 14 + A 2 ) + C 6t Y 19( Y 6 + 2A l) + Y 6( Y 19 + 2A 0 )] 
“ CyYigOV A l) + C 8 Y 5 (Y 6+ A l) - 2 Cio Y l(Yl9 + A q) 

- C 9 [Y 6 (Y 6 + 2A-|) + Y7(Y7+ 2A 2 )](Y 6 + A-|) - C-| 0 {Y 16 (Y 6 + A-|) + 4Y 9 (Y 13 + A.,) 

+ 4Y 2 (Y 2 o+ A.,) + Y 18 (Y 7 + A 2 ) + 4Y 11 (Y 14 + A 2 )+ 4Y 4 (Y 21 + a 2 )} 

*23= C 11 y 16“ C 12 Y 2 + C 205 Y 25 + C 206 Y 11 + c 13( Y 20( y 6 + 2A l) ~ 2Y 13( Y 13 + 2A -|) 

+ Y 6( Y 20 + 2 A i) - t Y 2l( Y 7 + 2A 2) - 2Y 14( Y 14 + 2A 2 ) + Y 7< Y 21 + 2A 2 )1> 

*24= ^1 y 17“ ^2 y 3“ ^3 Y 21“ C 4 X(Y 2 i+ A 2 ) + C 5Y7- C 2 qiY 22 - C 2 q 2 Y 8 

f 

“ C 203 y 27~ c 204 y 13 + C 21lPe( Y 7 + A 2) “ C 212 T ( Y 13 + A i) 

+ c 6t Y 19( Y 7 + 2 A 2) + Y 7( Y 19 + 2A 0 )1 " C 7 y 19( y 7 + a 2> + C 8 Y 5( Y 7 + k 2) 

- 2Ciq y 8 ( y i 9+ A q ) - C 9 [Yg(Yg + 2A-|) + Y7(Y7+ 2A 2 )](Y7+ A 2 ) - C-| o( Y 1 8 ( Y 6 ■*" A l) 

+ 4Y 1l( Y 13 + A i) + 4Y 4 (Y 20+ a') - [Y 16 (Y 7 + A 2 ) + 4Y 9 (Y 14 + A 2 ) + 4Y 2 (Y 21 + A*)]} 

*25= c 11 y 18“ c 12 y 4“ c 205 y 23" c 206 y 9 + c 13( Y 2l( Y 6 + 2A l) " 2Y 14( Y 13 + 2A -() 

+ y 7( y 20 + 2A i) + Y 20( Y 7 + 2A 2) " 2Y 13( Y 14 + 2A 2 > + Y 6< Y 21 + 2A 2 )} 


(G.2) 


(G.3) 


(G.4) 


(G.5) 
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*26= C 14 Y 19 - C 15 Y 5 - C 19 MY 19+ A 0 ) + C 17 [Y 6 (Y 6+ 2A^ 4 Y 7 (Y 7+ 2A 2 )] 

- C 16 {Y 20 (Y 6 + 2A-|) ♦ 2Y 13 (Y 13+ 2a') + Y 6 (Y 20 + 2a') + Y 21 (Y 7+ 2A 2 ) 

+ 2Y 14 (Y 14 + 2A 2 ) + Y 7 (Y 21+ 2a")} - C 18 {Y 15 (Y 6+ A^ + 2Y 8 (Y 13+ a') 

+ Yi(Y 20+ a") + Y 17 (Y 7 + A 2 ) + 2Y 10 (Y 14 + a') + y 3 (y 21 + a 2 )} 

*27= C 20 Y 20~ C 21 Y 6 + c 22 Y 15~ c 24 Y 1~ C 23 X(Y 2 q+ A^ + C 2 ogY 28 + C 210 Yi 4 + C 207 Y 24 
+ c 208 y 10 + c 213Pe(Y6 + A l) + C 214^( Y 14 + A 2 ) “ c 25[Yig(Y 8 + 2A-|) 

+ Y 6 (Y 19+ 2a")] - C 26 Y 19 (Y 6+ A-|) + C 27 Y 5 (Y 6 + A-|) - C 28 [Y 6 (Y 6 + 2A-|) 

+ Y 7 (Y 7 + 2A 2 )](Yg+ A-|) - 2C 29 Y-|(Yi 9 + A")- C 29 {Y 16 (Y 6+ A^ + 4Y 9 (Y 13+ a') 

+ 4Y 2 (Y 20 + a") 4 Y 18 (Y 7 4 A 2 ) 4 4Y 11 (Y 14+ A 2 ) 4 4Y 4 (Y 21 + a")} 

*28= C 20 Y 21~ c 21 Y 7 + c 22 Y 17 _ C 24 y 3~ c 23^(Y 2 i + A 2 ) - C 20 gY 27 - C 210 Yi 3 

/ 

“ c 207 y 22 - c 208 y 8 + c 213Pe(Y7 + A 2) " C 214*(Y-| 3 + A^) - C 2 5[Y-| 9 (Y 7 4 2 A 2 ) 

4 Y 7 (Y 19 4 2A")] - C 26 Y i g(Y 7 + A 2 ) ♦ C 27 Y 5 (Y 7 4 A 2 ) - C 28 [Y 6 (Y 6+ 2A-,) 

4 Y 7 (Y 7 4 2A 2 )](Y 7 4 A 2 ) - 2C 29 Y 3 (Y 19 4 a") - C 29 {Y 18 {Y 6 + A^ + 4Y 11 (Y 13+ a') 

+ 4Y 4 (Y 20 4 a") - [Y 16 (Y 7 4 A 2 ) 4 4Y 9 (Y 14 + A') 4 4Y 2 (Y 21+ a"]} 

Next one must evaluate analytically the components of the Jacobian matrix 
Introducing the short-hand notation 

3 . 

'J SYj 

one obtains the following values 

/ / / f t r f i t t 

J 1,8 = J 2,9 + J 3,10 = J 4,11 = J 5,12 = J 6,13 = J 7,14 = J 8,15 = J 9,16 = J 10,17 = 1 

FFFrrtFFFFt 

J 1 1,18 = J 12,19 + J 13,20 = J 1 4,21 = J 15,22 = J 16,23 = J 17,24 = J 18,25 = J 19,26 = J 20,27 = J 21,28 
Further 

J 22 t i = - ^2" 2 C-|q(Yi 9 4 A q ) ^22, 4 = ” 4 C-| q(Y 2 -| + A 2 ) 


(G.6) 


(G.7) 


(G.8) 


(G.9) 


(G.10) 


(G-1 1) 
= 1 


(G.1 2) 
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^ 22 , 2 = - 4 Clo( Y 20 + A -|) 


J 22,5” ^8( y 6 + A l) 


J 22,6 _ C 5 + C 21lPe + 2C 6( Y 19 + A 0 ) " C 7 Y 19 + C 8 Y 5” C 9[Y6( Y 6 + 2A l) 

+ Y 7 (Y 7+ 2A 2 )] - 2C 9 (Y 6 + A-,)(Y 6 + A-)) - C-| 0 Y 16 

f 

^ 22 , 7 ~ ~ 2 C9(Y6 + A l)(Y 7 + A 2) - C-|C)Y-|8 

J 22,9 = " 4C 10( y 13 + A -|) J 22,18 = “ C 1C)( Y 7 + A 2) 

/ 

J 22,10 = C 202 

f f 

J 22,11 = ~ 4C 10( y 14 + A 2^ 

f 

J 22,13 = " 4C 10 y 9 
/ 

J 22,14 = C 204 + C 212^ ~ 4C 10 Y 11 J 22,24 = C 201 


J 22,19 = 2C 6( y 6 + A l) “ c 7( y 6 + A l) ~ 2C 10 Y 1 
^22,20 = “ C3- C4A, ~ 4 ^10 y 2 


J 22,21" " 4C 10 Y 4 


J 22,15 


= Ci 


J 22,28 _ C 203 


J 22,16 = “ C 10( Y 6 + A l) 


Further 


J 23,2~ " C 12 

J 23,6“ 2C 13( Y 20 + A t) 

J 23,7 = “ 2C 13( y 21 + A 2^ 

/ 

J 23,1 1 = C 206 

' / 
J 23,13 = - 4C 13( y 13 + A i) 


J 23,14 _ 4C 13( y 14 + A 2^ 

/ 

J 23,16 = C H 

J 23,20 = 2C 13( y 6 + A l) 

J 23,21 = ” 2C 13( y 7 + a 2) 
/ 

J 23,25~ C 205 


Further 


^24,2” 4 Cl0( Y 21 + A 2) 


J 24,16 = C 10( y 7 + a 2) 


(G.12) 


(G.13) 


(G.14) 
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J 24,3 = ' C 2~ 2C 10( y 19 + a 0 ) 
J 24,4 = ~ 4 C 10( Y 20 + A i) 

^ 24 , 5 = ^ 8 ( Y 7 + a 2 ) 


J 24,17~ C 1 
/ 

J 24,18 = " C 10( Y 6 + A l) 

^24,1 9 = 2 C 6( y 7 + a 2) - c 7( y 7 + a 2) “ 2C 10 Y 3 


J 24,6 = " 2 C 9( y 6 + A l)( Y 7 + a 2) " C 10 Y 18 J 24,20 = “ 4C 10 Y 4 

J 24,7 = C 5 + c 21lPe + 2C 6( Y 19 + a 0 ) “ C 7 Y 19 + C 8 Y 5 

_ Cg[Yg(Ye+ 2A-|) + Yy(Yy+ 2A2)] - 2Cg(Yj+ A2)(Yy+ A2) + C-|oY-|6 


J 24,8“ “ C 202 

/ t 

^24, 9 = 4< ^10( y 14 + A 2^ 

^ 24 , 1 1 = - 4 ^ 10 ( y 13 + A -|) 

/ 

J 24,13 = " C 204~ c 212 x " 4C 10 Y 11 
/ 

^24,14 = 4C 10 Y 9 


^24,21 ~ ^3 C4 X + 4C-|qY2 

r 

^24,22 = - ^201 

f 

J 24,27 = _ C 203 


Further 


J 25,4 = " C 12 

^ 25 , 6 = 2 < ^ 13 ( Y 21 + A 2 ^ 

^25,7” 2 ^13( Y 20 + A i) 

/ 

J 25,9 = “ C 206 

/ // 

^25,1 3 = - 4< -'13( y 14 + A 2^ 

Further 

/ n 

^ 26 , 1 = “ ^ 18 ( y 20 + a i) 


J 25,14 = ' 4C 13(Yl3 + a i) 

/ 

J 25,18 = C 11 

^25,20 = 2 ^13( Y 7 + A 2) 

J 25,21 = 2C 13( y 6 + A l) 

/ 

J 25,23 = _ C 205 

J 26,13 = " 4C 16(Yl3 + a -|) " 2C 18 Y 8 


(G.14) 


(G.15) 


(G.16) 
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J 26,3 = ~ C 18( y 21 + A 2 ) J 26,14 = “ 4C 16(Yl4 + A 2 ) “ 2C 18 Y 10 

' / 

J 26,5 = ' C 15 J 26,15 = “ C 18( Y 6 + A l) 

J 26,6 = 2C 17(Y6 + a 1 ) " 2C 16 ( y 20 + a -|)~ c 18 y 15 J 26,17 = “ C 18< Y 7 + A 2) 

J 26,7 = 2C 17( y 7 + a 2) “ 2C 16( y 21 + a 2 ) " C 18 Y 17 J 26,19 = C 14 _ c 19*- 

J 26,8 = " 2C 18( y 13 + A l> J 26,20~ " 2C 16( Y 6 + A l) “ C 18 Y 1 

J 26,10 = “ 2 ^18C Y 14 + a 2 ) J 26,21 = ~ 2C 16( Y 7 + A 2) " C 18 Y 3 

Further 


J 27,1 = " C 24” 2C 29( Y 19 + A 0 ) 

r ft 


J 27 2 = - 4C 2 g(Y 20 + a i) 
J 27,4 = “ 4C 29( y 21 + a 2 ) 


J 27,14 = C 210 + c 214 t " 4C 29 Y 11 

t 

J 27,15 = C 22 

f 

J 27,16 = _ C 29( Y 6 + A l) 


J 27,5 = C 27( y 6 + A l) J 27,18 = " C 29( Y 7 + A 2) 

J 27,6 = “ C 21 + C 213Pe _ 2C 25( Y 19 + A 0 ) ~ C 26 Y 19 + C 27 Y 5~ c 18t Y 6( Y 6 + 2A l) 

+ Y 7(Y7 + 2A 2 )] - 2C 2 s(Y 0+ A-))(Yg+ A-|) - C 2 gY-|e 


J 227 = - 2C 2 8(Y 7 + A 2 )(Y 6+ Ai) - C 2 gY-|8 
J 27,9 = " 4C 29( y 13 + a i) 

r 

J 27,10~ C 208 

J 27,11 = “ 4( ^29( y 14 + a 2 ) 

l 

J 27,13 = ' 4C 29 y 9 


J 27,19 = " 2C 25(( y 6 + a 1 ) 

-C 26 (Y 6+ Ai)-2C 29 Yi 

t 

^27,20 = ^20“ C 2 3^ ” 4C 29 Y 2 
/ 

J 27,21 = ~ 4C 29 y 4 
/ 

J 27,24" C 207 

r 

J 27,28 = C 209 


Further 


J 28,2 = 4C 29( y 21 + a i) 


J 28,14 _ 4C 29 y 9 


(G.16) 


(G.17) 
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J 28,3 = ” C 24“ 2C 29( y 19 + a 0 ) 

/ // 

J 28,4 = " 4C 29( y 20 + a i) 

J 28,5 = C 27( Y 7 + A 2) 

J 28,6 = -2C 28( Y 6 +A l)( Y 7 +A 2) ~ C 29 Y 18 


J 28,16 = C 29( Y 7 + a 2> 

J 28,17 = C 22 

J 28,18 = " C 29( y 6 + A l) (G.18) 

J 28,1 9 = " 2C 25( Y 7 + a 2) _C 26( Y 7 + a 2) _2C 29 Y 3 


J 28,7" " C 21 + C 213Pe“ 2C 25( Y 19 + A 0 ) " C 26 Y 19 + C 27 Y 5 

- C 2 8[Y6( Y 6 + 2A l) + Y 7(Y7 + 2A 2)] “ 2 C 28 ( Y 7 + A 2)(Y7 + A 2) + C 29 Y 16 


l 28,8 = " C 208 

^28,20 = 4 ^29 y 4 

/ / 
*28, 9 = 4( -'29( Y 14 + A 2 ) 

r 

J 28,21 = C 20” C 23^ + 4C 29 Y 2 

*28,1 1 = “ 4C ^29( y 13 + a i) 

/ 

J 28,22 = " C 207 

f 

*28,1 3 = ~ ^210“ < -'21 4 T " 4 ^29 y 11 

/ 

^28,27 = ” ^209 


All other components of the Jacobian matrix J are equal to zero. Notice that since the Jacobian J is a 
function of Y (or U or V) therefore the variational equations (176)-(179) depend step-by-step on the 
results of the associated initial value problems. Thus the variational equations depend on the initial 
guesses S v . Hence the variational equations must be integrated together with the corresponding 
associated initial value problems. 

In shell buckling analysis the use of load increments to locate the limit point of the prebuckling 
equilibrium states leads to a singular problem. It has been shown by several authors^ 28,16 ! that with the 
introduction of the appropriate generalized displacements as an auxiliary constraint this singular behavior 
can be removed. It has been explained earlier (see Eqs. (152)-(153) and Eqs. (162)-(163)) that this 
auxiliary constraint condition can be written in the following form 


^ = f A (x, Y ; X, p e> x) 
dx 


A(0) = 0 
^ ) = A o 


(G.19) 


(G.20) 


where, if the variable load A = X, then A = 6 nl and 
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f A = C 50 Y 19 + c 51 y 5” c 52t Y 6( Y 6 + 2A l) + Y 7( Y 7 + 2A 2)1 

+ C 53{ y 12( y 12 + 2A 0 ) + j t Y 13( Y 13 + 2A -|) + Y 14( Y 14 + 2A 2 )]} 

whereas if the variable load A = p e , then A = W ave and 
f A — y 5 

and finally, if the variable load A = x, then A = y nl and 


f A - c 550 y 19 _ c 551 y 5 + C 552t Y 6( Y 6 + 2A l) + Y 7( Y 7 + 2A 2)J 

" c 553( y 13( y 7 + a 2) + y 7 a i ~ I Y 14( Y 6 + A l) + Y 6 A 2 ^ 

Considering the extended (29x29) dimensional Jacobian of Eq. (208) 


A' 3 A — A — - 

J = ^ f (x, 9 ; X, p e , x) 


then it can easily be seen that for the 29th column 

A' 3fj 

s ■ 0 

whereas for the 29th row the values of 


A' 3f A 

J 29,j = Wi 


depend on the generalized displacement used. 


Thus for A = 5 


nl 


J 29,5' C 51 


A' 


J 29,6 = " 2C 52( y 6 + A l) 


a' 


J 29,7 = " 2C 52(Y7 + a 2) 


a' 


J 29,12 = 2C 53( y 12 + a 0 ) 


A' 


J 29,13 = C 53( y 13 + a i) 


A' 


J 29,14 = C 53( y 14 + A 2^ 


A' 


J 29,19 = C 50 


whereas for A = w. 


ave 


(G.21) 


(G.22) 


(G.23) 


(G.24) 


(G.25) 


(G.26) 
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(G.27) 


A' cR 
J 29,5- ~ 

and finally for A = y n ' 

A' a' 

J 29,5 = " C 551 J 29,13" ' C 553( y 7 + a 2) 

J 29,6 = 2C 552( y 6 + A l) + C 553( Y 14 + A 2 ) J 29,14 = C 553( Y 6 + A l) 

J 29,7 = 2C 552( Y 7 + a 2) " C 553( Y 13 + A -|) ^29,19" C 550 

All other J 29 s= 0. 


(G.28) 
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APPENDIX H: The Forcing Function jL f and the Generalized Unit Vectors f 

dA ~ ~A 


The form of the forcing function df/dA of the inhomogeneous variational equations (204)-(207) 
depends on the choice of the variable load parameter A. Thus if A = A., then from Eqs. (G.2)-(G.8) and 
from Eq. (G.21) 

a? af 

dA dA, 


= {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, ?22> °- <24- °> <26- <27- <28- °} T 
where 


(H.1) 


o dioo " 

<22= = - C 4 (Y 20+ 


df 2 4 

IT 


<24= ■ - C 4(Y21 + a 2 ) 


<26=-^ = -C 19 (Y 19+ A 0 ) 


o dio-7 " 

<27= = -C 23 (Y 2 0 + a i> 

o dioo " 

<28= -gjj- = - c 23(Y21 + a 2 ) 


(H.2) 


_ A 


whereas if A = p e , then from Eqs. (G.2)-(G.8) and from Eq. (G.22), assuming that A. = p e R 


df df 


3a 


dp ( 


= {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, f 22 , o, f 24 , 0, f 26 , ? 27 , f 28> 0} 
where 


T 


(H.3) 


df. 


° ° ] 22 " & o " a 

<22= — — = “ c 4( y 20 + A -|)k + C 2 ii(Y 6+ A-)) ; f 27 = __ = - C 23 (Y 20 + A^R + C 213 (Y 6 + A-|) 


dp 

df. 


dp E 

3 28 


(H.4) 


Q UlOA "A o GOO " A 

<24= — - = ~ C 4 (Y 21+ A 2 )R + C 2 ii(Y 7 + A 2 ) ; f 28 = ___ = - C 23 (Y 2 -|+ A 2 )R + C 2 -| 3 (Y 7 + A 2 ) 

dp e dp e 

<26= -=^- = “ c 19( y 19 + a 0 )^ 

^Pe 


Finally, if A = x then from Eqs. (G.2)-(G.8) and from Eq. (G.23) 


df df 


dx (H.5) 

= {0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, f 22 , 0, f 24 , 0, 0, ? 27 , f 28 , 0} T 

where 


o dioo ' 

<22= — ^ = c 212( y 14 + a 2 ) 
dt 

o dfoA ' 

<24= — — = ~ c 212< y 13 + a i) 
dx 


o dio-7 ' 

<27= -4r- = C 214 (Y 14+ A 2 ) 
5t 

o dioa ' 

<28= = - C 2 1 4 (Y 13 + A.|) 

dx 


(H.6) 
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The form of the generalized unit vectors i , used as the initial condition for the inhomogeneous 
variational equations (204) and (207), depends on both the choice of the variable load parameter A and 
on the boundary condition one wants to enforce. Introducing the following notation 


I - generalized unit vector for the specified boundary condition (for instance, BC = SS-1 etc) 


i = JLf 

~X dk~BC 


I =i_f 

-P ap e ~BC 


t -£_t 

~BC 


(H.7) 


one has the following options. 


SS-1 Boundary Condition 

]_ ss _.|= (0, o, 0, 0, 0, 0, 0, 0, 0, 0, S-I , s 2 , s 3 , s 4 , s 5 , s 6 , s 7 , fy, fjj, fy, s 8 , s 9 , s 10 , s 1(l s 12 , s 13 , s 14 , 0} T 

? = (0, 0, 0, 0, -& v , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C 41 , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} T 

—A 


t 

~P 


{0, 0, 0, 0, -w c , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C 41 fi, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} T 


t = (0, 0, 0, 0, -w t , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C 22 o, o, 0, 0, 0, 0, 0, 0, 0, 0, 0} T 

~T 


where 


fu= - oV + WpPe + W t t) 


2 - 
f U = C 41^ + C 220 T 


f U = “ B 1 s 4- b 2 s 3 


4 

f|j = “ B 1 s 6 + b 2 s 2 


(H.8) 


and 


A 1 “ * 

w v =1a 12 


A 1 — * 

W P" -5- A 22 


A 1 — * 

W t = " 7: A 26 


A AAA 

W C =W V R+Wp (H.9) 


SS-2 Boundary Condition 


f 

-SS-2 


!x 


! P 

t 

~T 


{Si, S 2 , S 3 , S 4 , fy, 0, 0, 0, 0, 0, 0, S 5 , S 6 , S 7 , S 8 , S 9 , S 10 , S^, fy, fy, ty, fy, fy, f^, fy, S^, S^, S^, 0^ 

{0, 0, 0, 0, -w v , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C 41 . 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} T 

{0, 0, 0, 0, -W c , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C 41 A, O, 0, 0, 0, 0, 0, 0, 0, 0, 0) T 

{0, 0, 0, 0, -W t , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C 2 2Q, 0, 0, 0. 0, 0, 0, 0, 0, 0, 0) T 



102 


where 

,1 


A _ A _ 


f u= - (W V X + w p p e+ W t x) 


f u = B 5 s 10~ b 6 s 3~ b 7 s 6 + b 8 s 13 + 2B 9 a 1 s 5 


f u = C 41^ + c 220 x 


f U = “ b 1 s 8 + b 3 s 1 ” b 4 s 7 


f U= 2B 10 S 11~ 8B11S4+ BgAi S q - BgA 2 S7 


f U" " B 5 s 8 + B 6 s 1 ~ b 7 s 7 + b 8 s 14 + 2B 9 a 2 s 5 


(H.10) 


fu= - B 1 s io + b 3 s 3 + b 4 s 6 


f U = ” 2B 10 s 9 + 8B1-1S2+ B 9 A 2 s 6 + B 9A-1S7 


SS-3 Boundary Condition 


IsS~3 = {0 ' °’ °' °’ °' °' 81 ’ S2 ' S3 ' S4, S5 ' S6 ’ s? ' 4 f u- C *U' f u' f u' i s 8' s 9' s 10. Sn, s 12 , s 13 , s 14 , 0} T 


U = {°. 0. 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C 41 , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} T 


| = {0, 0, 0, 0, -W c , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, C 41 ft, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) T 


t 

-p 

A 


I = ( 0 . 0 , 0 , 0, -W t , 0, 0, 0 , 0 , 0, 0, 0, 0, 0, 0, 0 , 0, 0, C220. 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} T 

where 


,1 


A. — A _ 


f u = - (W V X + W pPe + W t T) 


f U = ■ 2B 17 s 2 


f u = B 15 s 3 + b 16 s 7 


f U = 2B 17 s 4 


6 

f u = c 41^ + c 220 x 


f U = B 18 s 3" b 19 s 7 


(H.11) 


f U = ~ B 15 s 1 _ b 16 s 6 


“ B 18 s 1 + b 19 s 6 


SS-4 Boundary Condition 


~SS-4~ {S1 - S 2- S 3' s 4- f u' °. s 5. s 6’ *. SB. SB. s 10 . Sl1 , f 2 l 3 f 4 f®. f®. I* I®. I® f™. fj 1 . f^ 2 s 12 , s 13 , s 14 , 0) 


t. 


{ 0 , 0, 0, 0, -w v , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 , 0, C 41 , 0, 0, 0, 0. 0, 0, 0, 0, 0, 0} T 


ip = (°* 0. -W c , 0, 0 ( 0, 0 r 0, 0, 0, 0, 0, 0, 0, 0, 0, C 41 A, 0, 0, 0, 0, 0, 0, 0, 0, 0. 0} T 

^ = { 0 , 0 , 0 , 0 , -W t , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , C 220 » o, 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 } T 

~T 



103 


where 

f J,- ' + WpPe + W t t) 

2 

f u = B 2 ()Si+ B 2 is 7 - B 2 2S-)i 

f U = 4B 24 s 2 + 2B 17 s 8 
4 

f U = B 20 s 3~ b 21 s 5 + b 22 s 10 

f u = 4B 24 s 4” 2B 17 s 6 

6 

f u = C 41 A, + C 2 20 t 


f U = B 25 s 1 + b 26 s 7“ b 27 s 1 1 

f U" B 25 s 3" B 26 s 5 + B 27 s 1 0 (H I 2) 

9 

f lT B 28 s 5~ b 29 s 3 + b 8 s 13 + b 30 s 10 + 2B 9 a 1 s 9 

fy°= 4B 3 -|S6" 8B 32 s 4 + b 9 a 1 s 10~ B 9 A 2 s 1 1 

f Jj 1= b 28 s 7 + b 29 s 1 + b 8 s 14 + b 30 s 11 + 2B 9 a 2 s 9 
12 

f U = 4B 31 s 8 + 8B 32 s 2 + b 9 a 2 s 10 + b 9 a 1 s 11 


C-1 Boundary Condition 


?c-V 


{0, 0, 0, 0, f^, 0, 0, 0, 0, 0, 0, 0, 0, 0, s-|, S2, S3, S4, Sg, sg, S7, Sg, S9, Sig. s-) -j , S12. S13, S14, 0} 


A 

I 

~P 

f 

~T 


where 


{0, 0, 0, 0, -W v , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0}' 

(0, 0, 0, 0, -W c , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} T 

{0, 0, 0, 0, -W t , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} 7 


f Jr - oV + WpPe + 'V) 


(H.13) 


C-2 Boundary Condition 

A 1 

[n_n m {S 1’ s 2' s 3’ s 4’ f u . 0, 0, 

A A 

I = {0, 0, 0, 0, -W v , 0, 0, 

t = ( 0 , 0 , 0 , 0 , -w c , 0 , 0 , 

~p 

t = { 0 , 0 , 0 , 0 , -w t , 0 , 0 , 

~T 


0, 0, 0, 0, 0, 0, 0, S5, sg, S7 Sg, S9, s-| q 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 


,2 ,3 A .5 _,T 

s 11' f u' f u- f u- f u' s 12' s 13’ s 14' °> 

0, 0, 0, 0, 0, 0, 0, 0, 0) T 

0, 0, 0, 0, 0, 0, 0, 0, 0} T 

0, 0, 0, 0, 0, 0, 0, 0, 0} T 



104 


where 


fi= - (W v ^ + W p p e+ W t x) 

2 

f u = B 10 s 7 _ b 11 s 3 + b 8 s 13 + b 13 s 11 
3 

f u = 2B 10 S 8~ ® b 1 1 s 4 


4 

f U= - B 10 s 5 + B 11S1+ B 8 s 14 - Bi 3 s 10 

5 

f u = " 2 B 10 s 6 + 8B 11 s 2 (H.14) 


C-3 Boundary Condition 


C-3 


-P 


{0, 0, 0, 0, f u , 0, 0, Sf, s 2 , s 3 , s 4 , 0, 0, 0, ijj, fy, fy f®, S5, s 6 , s 7 , s 8 , Sg, s 10 , s t1 , S12, s 13 , s 14 , 0} T 

{0, 0, 0, 0, -W v , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} T 

{0, 0, 0, 0, -W c , 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} T 

{0, 0, 0, 0, -Wj, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 f 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} T 


where 


fj* - (W V X + WpPe+ Wft) 
2 

f u = B 17 s 3 + b 8 s 6 
f u = 2B 17 s 4 


f u = ■ e 17 s 1 + b 8 s 7 
f u = " 2B 17 s 2 


(H.1 5) 


C-4 Boundary Condition 


Ic-4- 

(si. 

s 2 , 

s 3- 

s 4 , 

f U’ 

o, 

0, 

s 5’ 

s 6> 

s 7» 

s 8> 

0, 

o, 

0, 

< u 2 

'u 3 

C 

<u 5 

s 9. 

s 10- s 11* 

«!■ 

£ 

'u 8 

£ 

s 12> s 13» 

s 14> 

0} T 


(0, 




A 























lx 

0, 

0, 

0, 

-w v , 

o, 

o, 

0, 

0, 

0, 

0, 

o, 

0, 

0, 

0, 

0, 

o, 

0, 

0, 

0, 0, 

0, 

o, 

0, 

o, 

0, 0, 

o, 

0) T 

A 

(0, 




A 























1 

0, 

0, 

o, 

-W c , 

0, 

o, 

o, 

0, 

0, 

0, 

o, 

0, 

o, 

0, 

0, 

0, 

o, 

0, 

0, 0, 

0, 

0, 

0, 

0, 

0, 0, 

0, 

0) T 

~P 

























f = 

(0, 

0, 

0, 

o, 

A 

-W t , 

0, 

o P 

o, 

0, 

0, 

0, 

o, 

0, 

o, 

0, 

0, 

0, 

o, 

o, 

0, 0, 

0, 

0, 

0, 

0, 

0, 0, 

0, 

0) T 


where 

f U = " (VV + w pPe + w t*) f u = b 31 s 5" b 32 s 3 + b 8 s 13 + b 14 s 11 
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2 

, u = B24S1+ B 17 s 7 + B 8 s 10 

3 

f u = 4B 24 s 2 + 2B 17 s 8 

4 

f U = B 24 s 3" b 17 s 5 + b 8 s 1 1 
f U = 4B 24 s 4~ 2B 17 s 6 


f U = 4B 31 s 6" 8B 32 s 4 

Q 

f U = B 31 s 7 + b 32 s 1 + b 8 s 14“ b 14 s 10 (H.16) 

f U = 4B 31 s 8 + 8B 32 s 2 


Symmetry Condition at x = .1 Jl 


t 

~sym 

k - 


t 


f 

-X 


{Si, s 2 , S3, s 4 , s 5 , s 6 , s 7 , 0, 0, 0, 0, 0, 0, 0, Sg, sg, s 10 Sn, s 12 , s 13 , s 14 , fy, ijj, fy, «y, 0, f®, I®, 0) T 

{ 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 } T 

{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0} T 

{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0. i 1t 0, B 2 , 0, 0, 63, § 4 , 0} T 


where 


f u = B 33 s 10 + b 34 s 3 + b 35 s 14 + B 36 s 7 + b 37 a 2 t 
f U = 2B 10 S 1 1 “ 8B 1 1 s 4 

3 — 
f u = " B 33 s 8" B 34 s 1 “ b 35 s 13" b 36 s 6“ b 37 a 1 t 

4 

f U = “ 2B 10 s 9 + 8B 11 s 2 

5 — _ 

f u = ' B 38 s 10 + B 39 s 3 “ b 40 s 14 + ( b 41 + B 23 ' c ) s 7 + B 23 A 2 X 

0 _ 

f U = B 38 s 8“ b 39 s 1 + b 40 s 13" ( b 41 + B 23 T ) S 6“ B 23 A 1 T 


(H.17) 


and 


^1 =b 37 a 2 
^2= “ b 37 a 1 


^3 = b 23( s 7 + a 2) 
^4 = “ b 23( s 6 + A l) 


(H-18) 
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Free Edge Condition 


{0. 0, 0, 0, S1 , s 2 . S 3 , 0, 0, 0, 0, S 4 , s 5 , s 6 , s 7 , s 8 , s 9 s 10 , fj. f* f* Sl1 , s 12 , s 13 , s 14 , f« f u 5 , f ® 0} T 

{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, B50, 0, 0, 0, 0, 0, 0, B^, § 2 , § 3 , 0) T 

(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, § 4 , 0, 0, 0, 0, 0, 0, §5, § 6 , 67, 0} T 

{0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, -b 52 , 0, 0, 0, 0, 0, 0, 0, § 8 , 69, 0} T 


where 


f u = B 42 s 1 + b 50^ + B 5lPe“ b 52 t " B 53l s 2( s 2 + 2A l) + S 3(S3+ 2A 2 )] 
f u = B 12 s 2 “ b 2 s 6 “ b 1 s 7 


V B 12 s 3 + b 2 s 5 " b 1 s 9 

^ f * f 

f u = ( B 42 “ b 43 ^) s 4 “ b 43 A 0 ^ “ B 44[ s 5( s 2 + A l) + s 2 A i + s 6( s 3 + A 2) + s 3 A 2 ] 

5 ' _ _ 

f u = ■ B 1 s 11 + B 45 s 9 + ( b 46 “ b 47 ^) s 5 “ b 47 a i^ + ( B 48 + B 49 T ) S 3 + B 49 A 2 T 

6 ' _ - _ 
f u = " B 1 s 13" B 45 s 7 + ( b 46 “ b 47 ^) s 6 ~ b 47 A 2^ " ( B 48 + B 49 T ) S 2“ B 49 A 1^ 


(H.19) 


^1 = _ b 43 ( s 4 + a 0 ) 
6 2 = - B 47 (s 5 + 

^3 = " b 47 ( s 6 + a ?) 


^ 4 ~ b 50 ^ + b 51 
^5 = ~ b 43 ( s 4 + a 0 )^ 
^6 = ” b 47 ( s 5 + a i)^ 


^7 = - b 47(S6 + A 2 )A 

^8 = b 49 ( s 3 + a 2) 

^ 9 = - B 4 g(s 2 + A-|) 


(H.20) 


Mixed Symmetry or Anti-symmetry Condition at x = _ _ 

2 R 

f • H = 

- mixed 


{si. ®2> S3> s 4> 0, S5, Sg, 0, 0, S7, S0, Sg, 0, 0 , s-jq, s-j-j, 0, 0, 0, s 4 2 . Sl3> 0, 0, s-| 4 , 0}”^" 


A A A T 

I p = Ix = ( 2 J 
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APPENDIX I: The Extended Jacobian j(S v ) - i_ & (S v ) 

' a| ‘ ' 

Using elastic end rings to impose the boudary conditions at the shell edges and employing Newton’s 
Method with increments in deformation to obtain the location of the limit point involves the use of the 
following extended Jacobian J(S V ) given by Eq. (196) 



9gi 

9gi 





9gi 

9s-| 

9s 2 8 

0 


0 


9A~ 

9gi4 

dgi4 





9gi4 

9si 

9s 2 8 





9A 

9<P1 

9<pi 

9(Pl 

9<Pi 



9cpi 

9s i 

9s 2 8 

9pi 

9P28 

0 


IT 

9(p 2 8 

9(p 28 

9cp 2 8 

9cp28 



9(p 2 8 

9s-j 

9(p 2 8 

9pi 

9P28 



dA 



9(p 2 9 

9cp 29 

9q>29 

9(p 2 g 

9(p29 



9pi 

dP28 

3ti 

9*28 

9a 

0 


9<P56 

9<P56 

3<P56 

9<P56 

9<P56 



9pi 

9P28 

9*1 

9*28 

9A 





9h-, 

3h -] 

9hi 





9H 

9*28 

9A 

0 


0 


9h i4 

9hi4 

9hi4 





9*1 

9*28 

3A 

9A 

9A 

9A 

9A 

9A 

9A 

9A 

9si 

9s 2 8 

9pi 

9P28 

9*1 

9*28 

9A 


( 1 . 1 ) 


A careful look at the structure of this block-diagonal matrix reveals that it is assembled out of the results of 
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different computations. The components involving the derivatives of the matching and constraint 
conditions, the 3cpj/5Sj terms, are obtained by solving the appropriate 29-dimensional variational equations 
(see Eqs. (200)-(203), whereby the 29th components are used to fill the last horizontal row. The 
components of the last column are assembled out of the solutions of the corresponding inhomogeneous 
variational equations (see Eqs. (204)-(207)). Finally those components of this Jacobian which involve 
derivatives of the specified boundary vectors g and h can be calculated analytical, as will be described 
in the following. 

Limiting Case as E r 0 

in this case one must satisfy the 14-boundary conditions specified by Eq. (143), the first of which 
reads 


9l = b-j i + * 3 -|2 ^x 0 + ^9) ~ < w o + W v* + W pPe + Wfc) = 0 


Replacing the shell variables involved H * and M s by their equivalent expressions in terms of the 

u 

variables used in the anisotropic shell analysis W, W and F derived earlier (see Eqs. (99) and (107), 
respectively), one obtains after some regrouping the following nonlinear boundary condition 

91= - 2c ^ n 2 b 1 ° 1 [(Y 13+ A')Y 1+ (Y 14+ A 2 )Y 3 ] + (2c b° 2 JlL -1)Y 5 

A 22 


.2 21 rh o A v 1 . o 


D f 

[b-| U |A 1 Y 6 + i. b 1 ° 2 (Y 6+ 2A 1 )Y 6 ] - cn 2 _£1 [b°A^Y 7+ .1 b° 2 (Y 7+ 2A 2 )Y 7 ] 


B 

- 2cJ! n 2 b°[Y 6+ A-|)Y 8 + (Y 7 + A 2 )Y 10 - -L {J± - 2A.)Y 12 ] (| 2b) 

n A 22 

B + b + B * 

- cn2b l°l^- K Y 6 + A l) Y 13 + (Y 7 + A 2 )Y 14 ] - (6; 1+ JL1 1 ) (b° 2 Yi 9+ b°Y 26 ) 

A 22 A 22 

- 4cf. b 1 °,XA 0+ 2.5. b 12 (B 1 + 1 X + B 21 p e - B^x) * b^q - (W V X + W p p e+ Wt) = 0 
Next one can proceed to calculate the partial derivatives 3g-|/3Yj analytically yielding the following nonzero 


- -2c f n 2 b°(Y 13 . a') 


- - cn 2 b ,°[2 ft Y 3 . (Y 7 + A 2 )] 

A 22 


- - 2c T "VllYu- V 


991 h 0 * B 21 B 21 , 
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391 ^ R h 0 B 21 i 

= 2c b-, n -1 


301 


3Y 5 t 


12 ~ 


3Y 


22 


26 


- - bii(D. 


- * - * 
o -* B 21 B 21 


11^11 


'22 


^ = - 2c J! r|2 b 1 ° Y 8 - cn 2 [b 1 °,(Y 13+ a') + b 1 ° 2 (Y 6+ A-,)] 

A 22 


= - 2c ^ n^Yio- ^ 2 -il ^(Y^* A 2 ) + b° 2 (Y 7+ A 2 )] 
/ A 22 


3 9l 

3Y^ 


= - 2c f n 2 b 1 ° 1 (Y 6+ A-|) 


391 

3A” 


R u° 


0 R — * — A 

u ir\ n n ^ \ A 


= " 40 11 b 11 (Y- )2 + A 0 ) + b 12 (2 11 B 11+ q) - 


W„ 


- - 20 " „Vl(Y 7 * A 2 ) 


3g i _ , 9g i> 6 o R K°5* \h 
^ = ( ir )B T b 12 B 21" W P 


= 2c J! b°( ^£1 - 2X) 
3Yi2 t 11 T* 

m 22 


39 ^ = - 9 R h 0 R * - W t 

__ ^ — d 12 b 61 w t 


391 

9Y 13 


- cn 2 b 1 °[2 Y,. Ill (Y 6 . A,)] 

*22 


Proceeding similarly with the other boundary conditions specified by Eq. (143) one obtains for 


0 1 1 s ■. o ... s 


g 2 = b 21 H o+ b 22 (M; o+ Xq)-w o _-=0 


(1.4) 


the following nonzero partial derivatives 


^ * - 2c 4 n2b 2°< Y i3* *;> 


3g2 

9Yl4 


- cn 2 b 2 °[2 J! Y 3 h 


B 


21 


'22 


(Y 7 + A 2 )] 


^ - - 2c T n2b 2°l< Y 14* A 2> 


992 

9Y i9 


— * — * 

. 0 “* B 21 B 21 , 
b 22< D 1 1 + — — — ) 

A 22 


(1.5) 


3g 2 R K o B 21 

9Y 5 t 22 -* 

m 22 


9g2 

3Y 2 6 


- » 2 > 11 ‘ - 41^1 ) 


'22 
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302 

3 Y 6 

302 

^7 

ag 2 

avr 


2c _ n 2 b 21 Y 8 cn 2 —ill [b 2 ° 1 (Yi 3 + + b 22 (Y 6 + A-j)] 

Ar 


B . 


.2 u _21_ ri ^o 

i * 

'22 


2c -5. n 2 b 2 °,Y 10 - cn 2 _£1 [b 2 °(Y 14+ A 2 ) + b 22 (Y 7+ A 2 )] 

Ar 


Br 


,2 *21 r .° 

i * 

'22 


_R „ 2 k o 

t 


2c 11 n‘b 2 '' 1 (Y 6+ A-)) 


^jr - - *= T b 2° (Yl2i A o> * b 2°2< 2 T B 1*1* ® < l 5 > 


^-- 2 c T " 2 b 2 > 7 * A2) 


^ 2 " b 2 “ 

3p e t 22 21 




R uO , B 21 
A 22 


392 - - 2 R b°B* 

^T“ T b 22 B 61 


B 


21 


= - cn 2 b 21 [2 R Y r 

° y 13 1 t T * 

m 22 


(Y 6 + A-,)] 


the third condition 


93- b ’l N xV b 12 N x S Vl * b 13 H 1 S * b i>xV "I*' 0 


( 1 . 6 ) 


yields the following nonzero partial derivatives 
^ - - n 2 b - 4c B n 2 b 1 ' l 3 (Y, 2 . A^ . 2 B n^B,*, 


903 

aTJ 


= - 8c .5. n 2 b 1 1 3 (Y-| 3 + a') 


503 

3y 15 


2 y b > 2 , 


383 -4 B n 3 bli*- H. A* 


3 Y . 


f " 13 16 Y "16 


303 

3Y 17 


2 B nb^V B 6 *,) . A. a 2 * 6 


303 

3 Yj 


8c f. n 2 b 1 1 3 (Y 14+ A 2 ) 


303 

3 Y ^ 


20 


b 14°11 


^ ’ - 40 T n2b '3 Y 9* AVlV * <Yl2* V ^7 





Ill 


= b 1 1 3 (2n 3 D 2 * 6 - 4c — n 2 Y 11+ 4c — nx) + 1 _L - B 62 
3Y 7 13 26 t 11 t 2 R c 62 


-JL - nb 12 + 2 n 2 b 13 (B 11 2B 66 ) + 1 (A 12 + A 66 ) 


= - 2 R b'l*- _L A* 

^22 t 13 21 2 22 


393 O R nh 1 5* 
__=2_nb 14 B 61 


. 1 =;* 11 1 H* 

13 11 p r p 21 
^ n cn* 


^ ■ - 4c T n2b lV7* *2) 


4c — b 1 g(Y-|3+ A^) 




.2*1 ro„ R 


R X, 1 1 t ,= * 


n b-j 3 [8c 1L Y 2 - (D 12+ 4D 66 ) ♦ 4c J2JL] - - _L JL (B^- 2B 66 ) 


n n 


- - 8c " n ^4- 2" hVl'e ^ ■ 4c T " b '3 (Y7 " Aa) 


Further, the fourth condition 


. 1 s .1 s .1 ,,s .1 s s . 
94= b 21 N x/ b 22 N xvi + b 23 H 1 + b 24 M xr V 1 = 0 


yields the following nonzero derivatives 

= - n2b 21 ‘ 4c y n2b 23 (Y 12 + A o^ + 2 y p2b 24^11 + y A 12 


^94 r p 1 ' 

^ - 80 " n b 23< Y 13* A,) 


iSL • - 2 Ji bl B 2 *,--La* 

3Yi5 t 24 2 1 cn 22 
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394 * R n 3 k 1 5* 

W^ T b 23 B 16 


9Y^7 ” 2 T nb 23 (2B 26 B 61 ] 


^94 


R r 2^ 1 

t 


8c _ n t b 2 ' 3 (Yi 4 + A 2 ) 


^94 

9y 20 


- b' 6*. 114. SI 

24 11 2 R cn 21 


^ = - 4c B n 2 b 2 1 3 Y 9+ n 2 b 2 1 4 D 12 - 1 - 


1 \ n R* 

2Rc 22 


^ = b 23^ 2n 3 °26 - 4c j n2 Y 11+ 4cf nt) 


3g 4 . .1 -* 

_ - -4nb 23 D 16 


^ * n 4r 2 j n2b 23< B i‘r 2B 66 


- - 2 « b 16 * 
3Y 22 t 23 21 


^94 R 9 1 

- - 4c " n 2 b 23 (Y 6 + A,) 


394 

3Y 2 7 


u 1 Pi * 

b 23°1 1 


3g 4 o R nh 1 n * 1 7 * 

__ - 2 _ nb 24 B 61+ _ A 26 


394 R 1 ' 

_ - -4c_b 23 (Y 13+ A.,) 


3g4 R 9 1 

■3777 - - 4 = - n 2 b 23 < Y 7 * A 2> 


394 = , 394. ^ 
3Pe » 


3g4 p pi 

= - 4c _ n b 23 Y 1 


394 R 1 

! 4c -1- nb 23 (Y 7 + A 2 ) 


3t 


394 

3Y-I3 


,2l, 1 


= - n-b 23 [8c 11 Y 2 - (D 12+ 4Dg 6 ) + 4c B A] 


1 n 2 


8c — n 2 b 23 Y 4 2nb 24 D 16 + 1 .1 B 26 
Further, the fifth condition 


b >xV b i N x S Vl * » 33 H 1 S * b 34 M xV w 1 S - 0 


(1.10) 
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yields the following nonzero partial derivatives 


^ - - "*1*31- 40 T n2b 33< Y 12* V * 2 J " 2b 34 g 1*1 


^.-ec^^Vu. a;> 


T " 2b 33 Y 4- 2 " b 34 D 16 


5 95 _ . R 3. 1 r * 
3Y 3 t 33 16 


i?L - - 2 n. b ' 5 * 

3Y 15 t 34 21 


^ - 8 b x " 2 b 33 < Y 14 * * 2 > 


^■- 2 7 " b 33< 2i 26- S 61> 


^ ‘ - “b T n2b 33 Y 9* " 2b 34 D 1*2- 1 


995 

9 y 20 


1 - * 
»34 D 11 


^ • b ^(2n 3 D 2 * 6 - 4c £ n 2 V ir 4c f. m) 


995 . . 1 - * 

= - 4nb Q oD 1c 

3Y 21 33 16 


^ - nbg 2 + 2 " n 2 b^ 3 (B,V 2l^) 


*!=-.- 2 f b '6 * 

3Y 22 t 33 21 


^■- 4 c £ n^fYe.A,) 


995 


3 Y 


27 


b 33°1 1 


395 O R nh 1 R* 

3Y^ = T 34 61 


3g 5 R U 1 /v a\ 

_ = 4c _ b 33 (Y 13 + A.,) 


395 = - 4c — n 2 b 33 (Y 7 + A 2 ) 


3 Yn 


t 


995 , 995 


3p- 


a* 


995 

9 Yl 2 


- 4 b " b 2 b 33 Y 1 


’ 4b x " b 33 (Y 7’ A 2> 


995 

9^13 


»2|1 ro_ R 


R X 


= - n b 33^ c T Y 2" (D 12 + 4D 66> + 4C T ^ 


(1.11) 
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Further, the sixth condition 


g 6 = b 41 N Xi + b 42 N ^ 1 + b 43 H 1 S+ b 44 M x S -| " w 1^= 0 
yields the following nonzero partial derivatives 

wf . - n 2 b^,- 4c JjL n 2 b) 3 (Y 12 * A Q ) ♦ 2 £ rA^B,* 


dgfi R o 1 ' 

--8c^n 2 b 43 (Y 13 * A,) 


^96 


9Y 


14 


R .2.1 

t 


8c 11 n"b 43 Y 4 - 2nb^ 4 D; 6 


3g 6 _ . Ft 3. 1 r* 
- 4 T n b 43 B 16 


a 96 


2 T b « 5 21 


- - 8c ~ r. 2 b^(Y 14 . A 2 ) 


396 


3Y 


17 


2 « nb 4 1 3 (2B 2 g- B 6l) 


= - 4c y n 2 b^ 3 Y 9+ 


^ = b 43^ 2n 3 °2*6 - 40 y n 2 Y 11+ 4c J± nt) 


396 

aY 20 

396 

9YiT 


b it 5 1*1 


4nb 43 D 16 


3g l ^K 1 .o R .2.1 to- 
re 


_ - nb 42 + 2 11 n^b 43 (B ir 2B 66 ) 


396 

9Y 22 


- 2 R h ^ R * 
2 — b 43 B 21 


= - 4c y n 2 b 4 1 3 (Y 6 + 


396 

3Y 27 


b 43°11 


396 9 R nh 1 5* 

3Y^q = T b 44 B 61 


- - * ? -jw *;> 


396 .-4c^n 2 bJ 3 (Y 7 .A 2 ) 


3Y 


11 


396 , 3g6 


e 


3p- dX 


3g 6 , R 2. 1 v 

^5" 4c T n b "3 Y l 


*p- - ” nb 43 (Y 7 * A 2 ) 


dt 


( 1 . 12 ) 


(1.13) 
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- - „ 2 b>c " y 2 - (6,V 45 6 * 6 , * 40 R JL] - 1 

Further, the seventh condition 


.2 . . s .2 .. s s „ 
97 =b 11 N x 2 +b 12 N xy 2 - u 2 = 0 


yields the following nonzero partial derivatives 


3 97 . 2 h 2 

= - 4n 

5Y 2 11 


3g 7 _ 1 t , v . . 

3Y 13 8 R 6+ ^ 


a 97 = _ 2n T * 


3Y. 


16 


3 97 1 t r\j » \ 

8"R 7+ 2) 


3g 7 _ 1 t „ ' 

3Y 6 8 R (Yl3+ ^ 


3 g 7 
3Y -j q cn 


1 A* 
a 26 


897 1 t „ A \ 

s;"ir (Y ' 4 * a 2> 


997 

3Yii 


4cn 


L_ A,‘ 
2 22 


3Qt p 1 — * — * 

__ = 2nb 12 + _ (A 12 + A 66 ) 

Further, the eighth condition 


.2 .. s .2., s s . 
g 8 = b 21 N x 2 + b 22 N xy 2 v 2 ' 0 


yields the following nonzero partial derivatives 


398 .- 4n 2 b, 2 .,2ilA* 


9Y 


21 


12 


398 Onh 2 

Wg = 22 


3g 8 1 t n(V A , 

a?i ’ ' 4 R n<Ye * Al) 


998 = 1 A + 


9Yn 


26 


= 1 _L n(Y 7 + A 2 ) 
3Y 7 4 R W 2; 


^98 

3Yl6 


1 A* 
2cn A22 


(1.14) 


(1.15) 


(1.16) 


(1.17) 
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Further, the nineth condition 


_ u,3 s u 3 k ,s . 3 u s .3 . . s s - 

g 9 b 11 N x 3 + b 12 N xy 3 + b 13 H 2 + b 14 M x 2 u 3 = 0 


yields the following nonzero partial derivatives 


3gg 

W7 


4 — n 3 b 1 3 oB 1 * fi + - A* 
l 13 16 q 16 


^99 ry R . 3 ,An * D + \ 2 A * 

_ = 2 _ nb 13 (2B 26 - B 61 ) - _ A 


26 


n2b 1 3 3< Y H* A 2 ) 


Sg 9 - . 2 R h 3 B * 
3Y^‘ 2 T b 1“ B 21 


^ - - n 2 b 3 - 4c R n 2 b, 3 3 (Y 12 » A 0 ) * 2 R 


3gg 

3Y^ 


R 

t 


8c_n-bf 3 (Y 13+ 


1^20 ' 4nbl33D l* 3 ' 1 Ti ^ {2B ^~ B ‘* 


26 61 


^99 

3Y£ 


b, 3 3 ( 2n3 D 2 * 6 . 4c R n 2 Y,,. 4c R nt) - 4- 4 ^ B f * 


2 W c «2 


^ ' 4c x n2b f3 Y 9* " 2b 1 3 4 5 lV ^l 2 * A o> - b 1 3 4 5 1*1 


399 o R nh 3 o* 

" ‘ T 14 si 


3gg 

3Y 24 


2 — b? 0 Bo«- — L A * 


13 21 


cn 


22 


3g 9 R _2 k 3 , v . . 

__ - 4c _ n b 13 (Y 7 + A 2 ) 


iEL - - b?,5;,. 1 1 _L 5,*, 

3Y 28 13 11 2 R cn 2 21 


-^2- nb 12 + 2 .5. n 2 b 13 (B 11 2B 66 ) + 1 (A 12 + A 66 ) 


399 R 2u 3 /v/ . v 

w _ = -4c_n 4 b 13 (Y 6+ A-j) 


399 _ . R . 3 , v ' 

- 4c _ b 13 (Y 14 + A,) 


- - 4° T n2b 1 3 3 Y3 ‘ Tf <Y? * A2) 


3gg 3gg a 

4^- - r 

3p e dX 


(1.18) 


) 


(1.19) 
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9gg 


R „2.3 


.3 - * 


- 8c J n ^13 Y 4 + 2nb 14 D 16 


S- = - 4c f. nb* 3 (Y 6 + A-,) 
7 )t ^ 


^99 


.2.3 


R 


R X 


n b i 3 [Sc 11 Y 2 + (D 12+ 4D 66 ) - 4c 11 11] - ^ 


1 

2c 


R 


(B 22 2B rr ) 


66 ' 


Further, the tenth condition 


.3 .. s 

9lO = b 21 N x 3 + 


,3 ..s . 3 u s . 3 s s _ 

b 22 N xy 3 + b 23 H 2 + b 24 M x 2 " 3 “ 0 


(1.20) 


yields the following nonzero partial deviatives 


3 9lO 
" 9 Y T 



3^3 „ * 
b 23 B 16 


9910 


,2.3 


R X 


— = n bl(8c _ Y 2 + (D 19+ 4D rr ) - 4c _ 11] 

9Y-|4 23 t ^ 12 66 t 2 J 


9910 

~ 3 Y^T 


8c J1 n 2 b 23 (Y 14 + A 2 ) 


a 9 lO 

9 V?i 



3 

23 




) 


a 9 lO 

~W 


n ^ b 21 ~ n ^ b 23 ^ 12 + A 0 ) 


* 2 t 



■ - 80 T n2(Y ' 3 * A '> 


9910 „ R .3 1 t* 

3Y 17 t 24 21 cn 22 


9910 

W 


= - b 23 (2n 3 D 26 + 4c n 2 Yii+ 4c 11 nx) 


_R 

t 


99-10 _ . .3 -* 

■ ■ - — AnDooU^o 

9 Y 20 23 16 


9910 _ . R 2.3 v „2.3 p; * 1 1 t n p * 

■ — 4c — n brj q Y q + n b^^D-iQ* — + — — — — — B^n 

3Y 7 t 23 y 24 12 n 2 R c 22 


(1.21) 


9910 

W 



3 R* . 

24 D 61 



9910 

9 Yi? 


.3 n* 
b 24°1 1 


1 

2 


i _L i * 

R cn 21 


9 gio 


4c J1 n 2 b| 3 (Y 7+ A 2 ) 


9910 

9^24 


? R h 3 R * 
2 Y b 23 B 21 


9910 

9 Yl 0 


- nb 22 + 2 5. 
22 t 


n ^ b 23^ B 11 2 ®kr) 


J 66 ; 


9910 


3Y 


28 


.3 n* 
' b 23°1 1 
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a 9lO 


3Y 


= - 4c _ n 


R r. 2.3 


11 


b 23 (Y 6+ A!) 


4p u . 4 


0g 1O . R 2.3 v 
-=« 7 — = - 4c _ n^booYo 
3Y-(2 t 23 4 


a 9l0 a 9l0 


9Pe 


(_i!£) A 
ax ' 


9910 - - 8c 5. n 2 b 2 3 3 Y 4 . 2nb 2 3 4 5,V 1 ' B 


3Y 


13 


t 


R 26 


8910 = - 4c R nb| 3 (Y 6+ A-j) 


ax 


Further, the eleventh condition 


_ .3 Kl s l, 3 s ,3..s .3 ,, s s „ 

g 11 = b 31 N x 3 + b 32 N xy 3 + b 33 H 2 + b 34 M x 2 “ w 2 " 0 


{1.22} 


yields the following nonzero partial derivatives 

a 9l1 _ A R 3. 3 r * 

__ - -4 - n b 33 B 16 


"" - n^tsc $ Y 2+ <5,V «&> - 4c » A, 

n * 


a y i4 


t 


^7 - 80 T n2b 3 3 3< Y 14* * 2 ) 


a 9l 1 

9Y i5 


* 2 y bb 3 3 (2i 2 * 6 - B 6 V 




(1.23) 


= - sc 4 „ 2 4<Y 13t A ;, 


8911 - - 2 R h 3 B* 

a 7^7 ' 2 T 34 21 


a 9l1 

3Y 6 


b 33( 2n 3 °26 + 4c y n 2 Y i1+ 4c R nx) 


a 9n _ , ,3-* 
3Y 20 33 16 


“37“ = 4c -y n 2b 33 Y 9 + n2b 3 3 4 D 12- 1 


a 9ll 

^17 


h 3 -* 
b 34 D 1 1 


3g 11 R 3 -* 

17^ - - 2 _ nb 34 B 61 


3g 1 1 
3Y 24 


2 “f b 33 B 21 


W = 4C T p2b 33( Y 7 + A 2) 


3g 1 1 _ _ .3 ~ * 
3Yii = " ^ 33 ' 1 
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^11 - - nb 3 V 2 ^ " 2b M< 6 1*r 2B 66> 


a 9l1 _ A R .3 v ' 
— — - - 4c _ b 33 (Y-j 4 + A 2 ) 


3911 = - 4c *. n 2 bf 3 (Y 6+ A!) 


9Yn 


t 


3 9l1 = 3911 & 

* r 


3911 --4cf n 2 b| 3 Y 3 


9Yi2 


t 


3911 =-4c^nb| 3 (Y 6+ A-,) 


3t 


^!i - - 8c B n 2 b 3 3 3 Y4* 2nb 3 3 4 5,* 6 

Further, the twelveth condition 


„ k 3 m s k 3 ki s l3 u s .3.. s s _ 
9 12 = b 41 N Xq + b 42 N XVq + b 43 H 2 + b 44 M Xp~ W 2 ,U = 0 


yields the following nonzero partial derivatives 


a 9l2 . R 3. 3 ~ * 
__ = - 4 _ n b 43 B 16 


5g 12 „2 k 3 


9Y U 


- n‘b 4 “ 3 t 8c T V 2* < D 12* 4D 66> - 4c y 


■ 8b T n2b 43( Y 14* A 2 ) 


3912 -? R nh 3 f?R* R* 1 
- 2 _ nb 43 (2B 26 - B 61 ) 


- - n 2 b 4) - 4c B n 2 b 4 3 3 (Y 12 » * 2 * n 2 b£,i;, 


3 9l2 

Hy7 


= - 8c _ n 


R „2 k 3 


^43(^1 3 + 


3912 _ R 3~* 

97^7 " 2 T 44 B 21 


3g 12 _ ,3 R „2 n 


3 Ye 


= - b 43 (2n D 26 + 4 c _ n Yi -| + 4c _ nx) 


R 

t 


3 g 12 _ . . 3 -* 
9Y 20 43 16 


^7 ’ 40 T " 2b 43 Y 9 + n2b 4 3 4 5 1*2 


3912 


dY 


21 


b £* 6 i'i 


3g 12 _ _ R .3 -* 

~3Y^ T 44 B 61 


3gi2 

3^24 


p R h 2 R * 

2 — b 43 B 21 


(1.23) 


(1.24) 



(1.25) 
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3 9l2 

"377 


R „2 k 3 
t 


55 4c _ n t b4 3 (Y7+ A2) 


^!£--nb| 2 *2 Rn 2 b 4 3 3 (5,V2i 6 * 6 ) 


d 9l2 . R 2u 3 , v . . 

wtt ' T n b43< ^ 6 * Al) 


3 9l2 . R 2.3 v 

- "c _ n b 43 y 3 


^li --80 Rn 2 b 4 3 3 Y 4 »2nb 4 3 4 5; 6 

Further, the thirteenth condition 


913= b 1 4 1 N xV b 1 4 2 N 4- U 4 S = 0 


3 9l2 

3y 28 


h 3 -* 
b 43°1 1 


a 9l2 

ax 


- " 4c y b 43^ Y 14 + A 2 ) 


9 9l2 _ ^912 ^ £ 

3p e ^ 

3Qi p p q 

~--^T nb 43^ Y 6 + A 1> 
dx 1 


yields the following nonzero partial derivatives 


3 913 2n 7* 

3 9l3 

3Y 2 c Al6 

3Y 13 

3 9l3 , 2 h 4 

- - 4n b 4 4 

3 9l3 


w " 377 


i n ^ A *> 

i i <*♦ A ’> 


3913 

w 


i ff <Y 1 4 * A 2> 


3 9l3 

177 


i * < Y ' 3 * V 


3913 1 7* 

3Y-I6 crT 26 


3 9l3 

3Y 25 



a 9l3 

3Y11 



+ 


1 , 7 * T * 

_ (a 12 + a 66 


(1.25) 


(1.26) 


(1.27) 
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Finally, the fourteenth condition 


, 4 .. s .4 . . s s 
914= b 21 N X4 +b 22 N xy 4 - v 4 = 0 


(1.28) 


yields the following nonzero partial derivatives 


5 9l4 . 2l 4 2n 7* 

= - - — A 


dY 4 




12 


3 9l4 

W 


1 A * 
c 26 


3914 1 t , v A V 

_ - 7 _ n(V 7 . A 2 ) 


3914 - _ 9nh 4 

3V„ 22 


(1.29) 


3 9l4 _ 1 t , v . v 

dY 7 4 R ( 6+ 1 


3gi4 

3Y 18 


1 A* 

2tS A 22 


Limiting Case as E r 

In this case one must satisfy the 14-boundary conditions specified by Eq. (141), the first of which reads 


V W p . W t > - 3,°^ - = 0 


(1.30a) 


Replacing the shell variable H Q by its equivalent expression in terms of the variables used in the 
anisotropic shell analysis W, W and F derived earlier (see Eq. (99)), one obtains after some regrouping 
the following nonlinear boundary condition 


o * 

R b 21 


R r,2r 


h 1= a n Y 5 + (a 12 ~ 2c 11 )Y 12 + 2c 11 n [(Y-| 3 + A^Y 1+ (Y 14+ A 2 )Y 3 

t « * 1 


(1.30b) 


22 


,2 B 21 


+ (Y 6 + A-|)Y 8 + (Y 7 + A 2 )Y 10 ] + cn [(Y 13 + A 1 )Y 6 + (Y 14+ A 2 )Y 7+ A^* A 2 Y 14 ] 


22 


(D 


* B 21 B 21 
11 + ^^ 


)Y 26 + 4c (Y 12 + A q ) * 3°^ + Wp Pe + W t t) = 0 


'22 


Next one can proceed to calculate the partial derivatives ah-|/9Yj analytically yielding the following nonzero 
terms 


ah-| 

W7 


= 2cJln 2 (Yl3 + a') 


3h-| 


a° + 2c — 
12 t 


(2X - ^31) 

a 22 
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ahi r o ' 

- 2c J± n 2 (Y 14 + A 2 ) 


3Y 3 

9hi o 


t 


3hi r o 

^ ■ 2C T n Y 8 * cn‘ _il (Y 13 * A,) 


,2 B 21 


3h*i r o 

■gy- = 2c T n 2 Y 10+ cn- _ti. (Y 14 + A 2 ) 

' A 


-v 22 

,2 B 21 


22 


1 = 2c y n 2 Y-|+ cn 2 _il (Yg+ A-|) 


3Y i3 


22 

* 


3h -i r 0 

_ = 2c _ n Y 3+ cn- _11 (Y 7+ A 2 ) 


.2 B 21 


5h 1 FI* B 21 B 21 
= U 1 1 + — =T — 


22 


3Y 2 6 


'22 


3hi o a r 7 

ar ■ a n w v* * t (V i 2 4 a o> 


■ 20 T n2 < Y 6* Al ) 


3hi 3h-i a o A 

~=~ » (-grl) R + a ilWp 
3p e ^ 11 P 


- 2o nV 7 - a 2 ) 


3hi _o,a 
3t 


a 11 W t 


Proceeding similarly with the other boundary conditions specified by Eq. (141) one obtains for 


h 2 = a 2 °l( w o + w v + W p+ W t ) + a 2 ° 2 w 0 S _- (M x S o+ Xq) = 0 


the following nonzero partial derivatives 


ah 2 

a 2°r 

2 c . 

R B 21 

1 A 2 2 

9h 2 

3Yi9 

-* ' 
= D H + . 

IT * “ * 
B 21 B 21 

a 22 

3h 2 

3Y 6 

cn 2 

B 21 
A * 

(Y 6 + A-|) 

ah 2 
~5T = 

a °Wv- 

2 T 5 i‘,-o 



a 22 





3h 2 
3Y7 = 

cn 2 

B 21 

a 22 

(Y 7 + A 2 ) 

9h 2 

3Pe 

( Sh2 )^ 
dX ' 

* a 2 °lV 2 j B 2 *, 

8 h 2 

0 



3h 2 

0 * 

R “ * 

3y 12 

= a 22 



3x 

a 2 iWt+ 

2 T B 61 


(1.31) 


(1.32) 


(1.33) 
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The third condition 


h 3 = a^u®* a^v®* a i3 w -( S+ a i4 w ®- - N x ®- 0 


yields the following nonzero partial derivatives 
5h 3 - - n a 1 A* + n 2 

5Y7 - c a 12 A 12 + n 


3h 3 _ 1 t 1 = * 

9Y 14 c R 12 26 


3hg 

9Y 3 


n a 1 A* 
7 a 1 1 A 16 


3h 3 _ 1 It* 
W ^5 = cn a 12 A 22 


5h 3 t _ 1 /v a' \ a 1 / 1 1 t n ~ * \ 1 

- - 5 - a ii(Yl2 + A 0 ) + a 12 (_ + _ _ _ B 22 ) + a 13 


avi" r 


ah 3 _ _ 1 t nip* 

377 " 2 R c 11 62 


3h 3 


2 1 T* 

~ 3-i -i A 


9Y-|7 cn 


11 "26 


3h 3 

3Y 8 


1 a 11 (A 12 + A 66 ) 


3h 3 1 t 1 a 1 R * 

" ' 2 R ct I 2 21 


3h 3 


9Y 


10 


1 It* 
7 a 12 A 26 


3h 3 

3Y 2 i 


- — — a i 1 l(2Bo R - Boh) 

2 R cn 11V 26 61' 


9hq * 1 


Sh 3 

9Y 2 2 


1 a 1 A* 
— 2 a 1 1 A 22 


cn 


9hq -i * — * 1 


3h 6 1 t 1 Ip* 

77^'~2RZ5 " *1 


cn 


The fourth condition 


h 4= a 2i u -| S+ a 2 2 v -| S+ a^W®* a 24 W®_- N x S yi 


yields the following nonzero partial derivatives 


0h 4 _ n 1 t * 


9Yi 


22 1 2 


3h 4 

Mu 


1 1 a 1 R * 

c R 22 26 


(1.34) 


(1.35) 


(1.36) 
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n a 1 a * 

T 21 a 16 


1 1 

m 322 22 


4 4<y,2* a 0 > * 4(i * 1^1 4> ♦ 4 


9h 4 _ 1 t nig* 

9 \> " ' 2 R T 21 B 62 


9h 4 2 It* 

9 Y ^7 = ' cn a 2 1 26 


3h 4 1 1 * t * , 

„ - a 21 < A 12 + A 66 ^ n 


4 = _ 1 t 1 _ 1 5 * 
9 Y 2 Q 2 R cn 22 21 


dh 4 1 1 ~* 

_ = - _ a 22 A 26 


9 ^4 1 t 1 1 - * - * 

W— = - ± 4- — a ?l(2B PR - B R1 ] 

9 Y 2 i 2 R cn 21 26 61 


4^1 (Y 6 + AO 


1 1 A* 

— 2 a 21 22 
cn 


ll It *" 3 !' — * ■! 

gy— = a 2 ll — 2 + 2t H (B 22 " 2 E W 1 + a 24 


11 1 a 1 R * 

2 R — 2 a 2 1 21 
cn 


The fifth condition 


h 5 = a 31 u 1 S+ a 32 v 1 S+ a 33 w 1 S+ a 34 w Tx" 0 


yields the following nonzero partial derivatives 


dhc n 1 — * R o ' 

W 7 = "^ a 32 A 12 + 4 c T n ( Y 12 + V 


9h S _ 1 a 1 A* 
cn 32 A 22 


dhc R rs ' 


|f 44* 2 T " (2B 26“ B 61> 


^ - H 4 a * - 4 B n 3 i * 
3Y 3 c 31 16 t 16 


11 1 a 1 i* 

"2 "R cn 32 B 21 


dnc R o ' 

4 = 8 c T n (y 14 + a 2 ) 
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-Jl = - _L al(Y 12 + A ) + a L (1 + 1 _L - B* ) + aL+ 4c — n 2 Y 9 
3Yg R 31 ld 0 32 n 2 R c 22 33 t y 


-S- = " - i - a 31 B 62 + 4c — n2y 11- 2 " 3 D 2 V 4c — nT 

77 j 2 R C 31 62 t 11 26 t 


ah 5 11 


R r"i 2/0 + 


3h 5 1 t 1 1 /0 5 


"c a 31^ A 12 + A 66^ 2 — n(B 11 2B 66> 77^ ~ 2 "R cn a 31 (2B 26 B 61> + 4nD 16 




3h 5 

^'^ C,31 ^ 22 ^~ LJ21 


1 alA* + 2 *B,* 


ahs _ _ 1 1 a* 

3 Y^~ 7 a 32 A 26 


sy 27 2 R cn 2 31 21 11 


(1.39) 


3hc p o 

b = 4c " n 2 (Y 7+ A 2 ) 


3Yn 


t 


3 hn r ' 

b = 4c _ (Y-| 3 + A.) 

dx t 14 1 


3h5 = - JL aon(Y 6 + Ai) + 4c ii n 2 Y-| 

3 Y -12 R 31 b 1 t 1 


3 hs 3 hs 

3 p e » 


W13 ' a 31 [J 2 * i Tf ,5 22' 25 66>1 * 4.* 4c y ,2n2Y 2* W ‘ " 2(5 lV 4D K> 


3h 5 _ 1 1 J R* + fir R n 2 Y„ 

3Y^ ' c R a 32 B 26 + 8c T n Y 4 


^ = - 4c " n <V 7 * A 2 ) 
3x * 


The sixth condition 


h 6 = a^ Ul S + a 4 2 v -| S+ a i 3 w 1 S+ a 44 w lV M xV ° 


(1.40) 


yields the following nonzero partial derivatives 


fl - - JL a>* - 2 £ n 2 i* 

3Y-j c 42 12 t 11 


— — — = " — — a 4 ? B ?*fi + 2 nD i*fi 

3Y 14 c R 42 26 16 


3h 6 n 1 T * 
* c "116 


- -1 a’ A* * 2 £. 5* 

3Y 15 cn 42 22 t 21 


(1.41) 
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ah 6 

3Y f 


- - * 4< Y 12* \) * 44 * J TT 7 g 22> * 4- " 25 12 


3h 6 1 t nip* 

a ?7 2 -R c 41 62 


<% 

3Y 17 


2 1 7* 

cn 841 A 26 


ah 6 

3 Y 8 


■3- a 41^ A 12 + A 66^ 


3h 6 - - 1 1 1 J i* + R* 

SY-pq 2-R-5F 42 B 21 +D 11 


3he 


3Y 


10 


7 a 42 A 26 2 y n B 61 


ah P 


5 Y 2 i 


- — — aL(2B* - Bfii) 
2 R cn 41 26 6V 


dh 6 t 1 , v . > 

3Y 12 ‘ "R 4l( 6+ ^ 


5h 6 _ 1 a 1 A* 

^ ^2 41 22 


3hg 11 — * 1 

4 = ^ i * <B 22- 2B 6sH * 4 


ah 6 It 1 1r* 

9^7 = " 2 R ^2 a 41 B 21 


cn 


The seventh condition 


h 7 = a 2 u 2 S + ai 2 2 v 2 s - N x S 2 = 0 


yields the following nonzero partial derivatives 




dh 7 1 t 2 /v . . 

9Yn 8 R 11 6+ 


r 13 


3h? = 2 n a 2 A* 
3Y^ c 1H6 


ah 7 1 t 2 . . 

■8 -r a 11 (Y 7 + a 2) 


3Y 


14 


3h7 1 t 2 7 1 t 2 

_ = --^ 11 ^ 3 ^ 1 ) + l^nai 2 (Y 6+Al ) 


3h 7 _ 1 2~* 

3Y^ ~ 2HR 12 A 22 


3h7 1 t 2 7 1 t 2 

^ = 8 r a 11 (Yl4+ A 2 } - 4 -R na l2< Y 7 + A 2) 


ah 7 1 2 - * 

3Y^J = "-3H a 11 A 26 


3h 7 

3Yq 


_ ai1 (A 12 + A 66 ) 


3hy 

w^ 11 " 22 


1 2 7* 

Sh ^ A/* 


(1.41) 


(1.42) 


(1.43) 
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3h 7 

avTT 


1 2 - * 
_ a-,2 A 26 


(1.43) 


The eighth condition 


h 8~ 


2 S 
*| u 2 


a 22 v 2 


N 


xy 2 


(1.44) 


yields the following nonzero partial derivatives 


9hg 

3Y 2 



2 A * 
a 22 A 12 


8h 8 




Al) 


8h 0 

9Y 4 



2 7 * 
a 21 A 16 


3h 8 

3Yl4 




A 2 ) 


0h 8 

0Y 6 


* 4( y 13 + A i) + ~ na 2 2 2 (Y 6+ Al ) 


0h 8 

aY 16 


1 2 - * 
2^ a 22 A 22 


0hQ 

0Y 7 


1 t 
8 R" 


afi(V 14 H 


Ao) 


1 t 
4 "R 


na 22 (Y 7 H 


A 2 ) 


3h 8 

9Yis 


1 

cn 


2 A * 
a 21 A 26 


3h 8 

3Y 9 


1 

c 


a 2 i( A i 2 + a 66 ) 2 n 


3h 8 
BY 23 


1 

4cn 2 


2 a * 
a 21 A 22 


(1.45) 


3h 8 

3Yn 


1 2 ~ * 
_ a 22 A 26 


The nineth condition 


3 s 
Hl u 3 


3 s 
*12 v 3 


3 s 
*13 w 2 


3 s 
*14 w 2,x- 


v 


(1.46) 


yields the following nonzero partial derivatives 


3h 9 _ n 3 - * 
977 - ' c 11 16 


3hg 

3Y U 


, 3 r 1 

Hi 1 -? 


+ _L _L (b L- 

2c R 22 


2B 6 6 )] 


+ a 


14 


9h 9 

9Y 3 


n 

"c 


3 t * 
a 1 2 A 1 2 


+ n 


3h 9 _ 2 3~* 

3Y^ H 11 26 


(1.47) 
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ah 9 

av 6 


i_ _t_ 
2 'R 


n - 3 r * 
7 a 11 B 62 


dhg 

3Yi 7 


1 3 7* 

cn 9 12 A 22 


3hg 

3Y 7 


i 4< Y 12* A 0 ) 


- 44 



t n o * , 

_ _ B 22 ) + a 


3 

13 


dhg 

5Yi 


1 3 7 * 

_ a 12 A 26 


3hg 

3Y 20 


It 1 3 ~ * — * 

= i i ^ a lV 2 B 26 - B 61 ) 


2 R cn 


(1.47) 


dhg 

dYio 




+ A 66> 


9h 9 _ 1 t 1 3"* 

37^7 ’ 2 R cn 12 B 21 


dhg 


- 4 a^tYy* a 2 ) 


dhg 

3Y 24 


1 3 7 * 

— 2 3 11 A 22 
cn 


dhg _ 1 t 3-* 

‘ ' c R a 12°26 


dhg _ it 1 a 3-« 

3Y^ - " 2 R ^2 a H B 21 


The tenth condition 


. 3 s 

h 10" a 21 u 3 + 


a 22 v 3 + 


3 s 
a 23 w 2 + 


a 2 3 4 w 2^" 0 


(1.48) 


yields the following nonzero partial derivatives 


dh-io _ n 3 -* dh 10 

W = " c 3 21 A 16 3Y^ 


dh-10 n 3 ~ * dh 10 

"dYg" = c a 22 A 12 3Y^ 


dhio _ 1 t n 3 -* dh-,o 

9Y 6 2 R c 21 62 3Y77 


“ a 21 [ ~2 + Tf ^ B 22 2B 66^ + a 24 
n 

2 3t* 

= cn a 21 A 26 

= ' a 22 A 22 ( U9 > 


ah 10 

9Y 7 


Tf a 2>12* Aq) - 4(1 


1 _L n 

2 R 




dhl° 1 3 -* 

1 _ = - _ a 22 A 26 


3h 


10 


1 t 1 _ 3 - * - * 


- 4 4. _ apH(2B 9fi - B fi1 ) 
3Y 2 q 2 R cn 21 26 61 
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3h 10 1 3 “ * T * , 

= -_a 2 i( A i 2 + a 66 ) +n 


9h 10 1 t 1 3 - * 

aY^j" = 2 R cn a 22 B 21 


3h 10 t q 3 , v 

= ' R a 21 (Y7+ A2) 


dn 10 _ 1 a 3 7* 

’ 7? 21 22 


3h 10 _ 1 t 3 - * 

%"c R 3 22 B 26 


9h 10 _1t 1 3 “* 

9 YS ■ ‘ 2 R 72 a 2' B 21 


The eleventh condition 


. 3s3s3s3s..s„ 

h 11 = a 31 u 3 + a 32 v 3 + a 33 w 2 + a 34 w 2,7 _ H 2 = 0 


yields the following nonzero partial derivatives 


ah 11 - n o 3 a * . / 1 R n 3 S* 
9Yi c 31 16 1 16 


3 hi i p 3 — * R ~ * — * 

11 = 3oi A or - 2 _ n(2B OR - B Ci ) 

ay-15 cn 31 26 t v 26 61' 


Q 4 " 2 < Y i 4 * A 2> 


3h 11 - 1 n 3 a * 

WT7 ~ ~ cn a 32 A 22 


- £ 4*12* 40 T " 2 (Yi 2 * A o> 


3 h 1 *1 1 t 1 3 — * — * ~~ * 

-^ 7 — - - -FT — a qi(2Bp R - B fi1 ) - 4nD 1R 
3Y 20 2 R cn 31 26 61 16 


± 11.80 "nW M 


1 1 1 a 3 R * 

2 R cn 32 21 


_11 = 1 - -1 — a 3-|B02 + 4c — n 2 Y-| 1 + 2n 3 Dp 6 + 4c .5. n*c 

aY 6 2 R c 31 62 t 11 26 t 


^ 7 - - 4 4< Y i2* *„> - 44 * * * £ 4> * 4- 4c t " 2v9 


1 a 3 A * 
1 a 32 A 26 


±11 - _L a?,A* . 2 f. B * 

3Y 2 4 on 2 31 22 ' 21 


4c 1 n 2 (Y 7 . A 2 ) 


3h 11 _ It 1 3-* - * 

aVii" 2^72 a 31 B 21* D 11 
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^1 1 1 3 - * T * \ R — * Shi 1 R ' 

= ' S a 31 (A 12 + A 66 J " 2 — n ( B 11" 2B 66 ) -gj[— = 4c — ( Y 14 + A 2 ) 


£li - 4c * „ 2 <Y 6 * A,) 


9h 11 = , 5h 11 

9p e 


)A 


ah 


ii 


avi 2 


^ 4lW 7 + A 2) + 4c y n 2 Y 3 


a hi i o 

— =— = 4c 1_ n (Y 6 + A-,) 
3t 1 


(1.51) 


- - 1 JL a 3 B* + 8c — n 2 Y 4 
3Y 13 c R 32 26 t 4 


■\u 

3Y^7 S *31 ^ + lc i ( ® 22 " + a 34 _ 4c y (2n 2 Y 2 - \) - n 2 (D 12 + 4Dg 6 ) 


The twelfth condition 


h 1 2 = 


3 S 
a 41 u 3 + 


a 42 v 3 + 


a 43 W 2 + 


3 s .. s 
a 44 w 2,x" V 


0 


(1.52) 


yields the following nonzero partial derivatives 


ah 12 _ n - 3 I * 
"37T = " c 4 1 16 


ah 12 

3Yl4 


a 4l(— - 
41 1 o 


* 2? Tf <B 22- 2i 66)I * a 


3 

44 


ahi 2 

3Y 3 


n a 3 A* - 2 R n 2 R* 

_ a 42 A 12 2 y n 


ah 12 2 3 

3Y^5 = "crT a 41 A 26 


ah 12 1 t n 3 ~ * 
^Yg - " 2 R ~E 41 B 62 


ah 12 

3Y 17 


1 

cn 


a 3 a * + 2 R n * 

a 42 A 22 2 y b 21 


(1.53) 


ah 12 

SYt 


.3 ,1 


1 t n “ * 


■r a 41 (Y 12 + A 0 ) ' a 42 ( y + J y B 22 J + a 43" n2[) 12 


ah i2 

W 



3 T * 

42 A 26 




3h 12 _ 1 t 1 3 -* 

' 2 R cn a 41 (2B2 6“ B 61> 


a h 12 

^0 




A 66^ 


ah 12 

9Y 2 i 


111. 

2 R cn 
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a 42 B 21 


+ D 


* 
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ah 12 t 3 A 

* ‘ Tf a 41 (Y? * A 2> 


ah 12 _ 1 a 3-* 

Wil - ^2 a 41 A 22 


iHH. - - 1 1 a 3 S * - 2n6 • 
5Y 13 c R 42 26 16 


ah 12 

avii 


1 1 1 q 3 r * 

2 R ^2 a 41 B 21 


The thirteenth condition 


h,3- -,V *- <- 0 


yields the following nonzero partial derivatives 


3h13 - - 2 £ a. 4 .A.* 


3Y 2 


1116 


_ 13 = - 1 J. a 1 4 1 (Y 7 + A 2 ) 
3Y 13 8 R 11 7 4 


ah 13 0 n 4 t* . 2 
= 2 _ a i0 A H0 + 4n 

3Y 4 c 12 12 


^ - - 1 1 af.(YQ* A,) 
3Y 14 8 R 11V b 1 


ah 13 

3Y 6 


1 t .4 


1 t 


0h 13 1 4 t* 


1 j. a 1l( Vi 4 + A 2 ) - i ^ na 12 (V 7 . A 2 ) ^ a„A 26 


3hi3 

3Y 7 


^ ^ *i> - -1 4- "4< Y 6* A 1> 3,113 


4 R 


9Y 18 


1 4 -* 

2 w a 1 2 22 


3h 13 - 1 a 4 A* 

3Yi" c *12*26 


3hi 3 

3^25 


1 4 T * 

— 2 a 11 A 22 


4cn 


ah 13 1 4 ,7* 7 * , 

__ - - _ a 11 (A 12+ A 66 ) 


Finally, the fourteenth condition 


h 14 = a 4 ^ a 4 2 v 4 s - N^= 0 


yields the following nonzero partial derivatives 


ah 14 _ _ ? n 4 -* 

- 2 _ a 21 A 16 


3h-i 4 i ♦ 4 

= a oi( Y 7 + A?) 

3Y-| 3 8 R 21 7 4 


(1.53) 


(1.54) 


(1.55) 


(1.56) 
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3Y4~ 



4 A* 
22 A 12 


dhi4 

3YiT 


2 _t. 
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a 21^ Y 6 + A l) 


3h 14 

9Y 6 


a 21 (Y 14 + k 2> ~ \ na 22^ Y 7 + A 2) 


5h 14 1 4~* 

= cn a 2 i 26 


3h 14 

3 y 7 


~4< Y 13»A,'> 


7 i "4< Y 6» Al) 


3h 14 
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2cn 


4 - * 
a 22 A 22 


ah 14 4 - * 

- - - a 22 A 26 


8h 14 

9Y 25 
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4 7* 
a 21 A 22 
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a 21 (Al2 + A 66^ + 2n 
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Table 1 Comparison of calculated buckling loads - Khot’s glass-epoxy shell 

(-40, 40, 0) - N c ^ = -965.9814 lb/in,R = 6.0 in, L = 12.5 in, t = 0.036 in 



AN 1 LISA 

COLLAPSE 


ANILISA 

COLLAPSE | 

SS-1 

-282.324 (n=2) 

-282.323 (n=2) 

C-1 

-519.656 (n=12) 

-519.656 (n=1 2) | 

SS-2 

-282.702 (n=2) 

-283. (n=2) 

C-2 

-521.455 (n=12) 

-521.454 (n=1 2) | 

SS-3 

-518.116 (n=12) 

-518.115 (n=1 2) 

C-3 

-520.333 (n=12) 

-520.332 (n=12) 

SS-4 

-520.727 (n=12) 

-520.663 (n=12) 

C-4 

-521.597 (n=1 2) 

-521.597 (n=12) 


Note: Buckling loads are given in Ib/in 


SS-1 

: N x = -N 0 


Z 
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II 
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II 

o 
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II 

c 

o 
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II 
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II 
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SS-3 
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II 
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II 
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II 
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II 

o 

SS-4 
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II 

c 
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o 

II 

£ 

ll 

> 


M x =0 

C-1 

: N x =-N 0 


N x y = W = 0 


w, x = 0 

C-2 

c 

II 

c 

o 


z 

II 

II 

o 


jS 

X 

II 

O 

C-3 

: N x = -N 0 


o 

II 

=i 

II 

> 


w, x = 0 

C-4 

: u = u 0 


o 

II 

£ 

II 

> 


w, x = 0 


Table 2 Summary of Imperfection Sensitivity Calculations using ANILISA^ 
Khot's glass-epoxy shell (-40, 40, 0) - N c ^ = -965.981 41b /in 



onl 
K c 

n 

b 

a 

P 


% 

$2=0-1 

$2 =0.5 

$2=1-0 

SS-1 

0.292266 

2 

-0.0436 

1.0042 

1.0179 

44.97 

-59.08 

0.869 

0.677 

0.553 

SS-2 

0.292659 

2 

+0.0328 

1.0015 

1.0139 

44.97 

+17.67 

- 

- 

- 

SS-3 

0.536362 

12 

-0.3893 

0.8975 

0.6644 

44.12 

-122.26 

0.757 

0.436 

0.251 

SS-4 

0.539065 

12 

-0.1726 

0.9444 

0.8083 

44.09 

-93.68 

0.807 

0.542 

0.385 

C-1 

0.537958 

12 

-0.2654 

0.9468 

0.8267 

44.56 

-110.73 

0.781 

0.498 

0.340 

C-2 

0.539820 

12 

-0.1115 

0.9689 

0.8939 

44.55 

-73.98 

0.829 

0.591 

0.447 

C-3 

0.538658 

12 

-0.2176 

0.9525 

0.8403 

44.56 

-104.68 

0.793 

0.520 

0.364 

C-4 

0.539966 

12 

-0.1101 

0.9690 

0.8934 

44.55 

-74.24 

0.830 

0.592 

0.448 


Note: 0 C and 0 C are given in degrees. 









Fig. 2 
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Fig. 3 Determination of forces and moments at the ring centroid - at x = 0 (lower edge) 



Fig. 4 


Determination of forces and moments at the ring centroid - at x = L (upper edge) 
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Stringer stiffened shell AS — 2 -xibQr2«0.05,mai2,n« 1 4- ss-4 B.C 


di»ploc 0 m*nt W 1 — max vs axial load (both normalized) 



Fig. 6 Variation of IW, max l with axial load 
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Khot'a glose — epoxy shell ( — 4-0. 4-0,0) xibor2*0. 1 D-OS.nn* 1 ,n»=2 


W1— mox vs axial lood (both values normalized) SS — 1 0.C. 



a - W 1 max vs normalized axial load p s = Xg/X^ 


Khot'a glass — epoxy ahell ( — <40. 4-0,0) x!bor2 «0. 1 D — 05, t ,n*2 


end-* shortening (total) vs axfal load (both normalized) — SS— 1 S.C 



0.00 0.50 1.00 1.50 2.00 2. SO 

end — shortening (total) / 0.422100 


b. end-shortening vs normalized axial load p s = X s /X c 


Load vs deformation plots for SS-1 boundary conditions 
(N x =-N 0 , N X y=w=M x =0 ; N c? =-965.9814 Ib/in) 
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Khot'* gloiB-apoxy shell ( — 4-0. 4-0.0) x!tyor2“0. 1 D — 05,m» 1 .n— 1 2 
W 1 — max vs axial load (both normalized) — SS — 3 B.C. 



a. w 1,max vs normalized axial load p s = 


Khot'» gloss — epoxy shell (—40.40,0) xlbarS "0. 1 O — 05, m * 1 ,n» 1 2 . 



b. end-shortening vs normalized axial load p s = X s /X c 


12 


Load vs deformation plots for SS-3 boundary conditions 
(N x =-N 0 , v=w=M x =0 ; N c{ = -965.9814 Ib/in) 
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KHot* ■ glaaa-apoxy ■ h«Ti ( — -4-0, *0.0) xlbar2*0.1D“05,mw1,n»12 


■ nd-»hort*nFng (total) v * axfot load (both normalized) — SS — 4 B.C. 



and — ehortening (total) / 0.77644-1 


Fig. 14 Enlarged view of the limit point for SS-4 boundary conditions 
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Khot’s glass-apoxy shell ( — 40.40.0) xfbar2-0.1 D-05.m- 1 ,r> — 1 2 
W 1 — max vs axial load (both normalized) — C — 3 B.C. 



a. wi max vs normalized axial load p s = Xg/Xc 


Khot'a gloss-epoxy shell (-4-0, 4-0.0) xiba r2 “O. 1 O — 05. m * 1 .n • 1 2 


end — shortening (total) vs axial load (both normalized) — C — 3 B.C. 



b. end-shortening vs normalized axial load p s = Xg/Xc 


Fig. 17 


Load vs deformation plots for C-3 boundary conditions 
(N x =-N 0 , v=w=W, x =0 ; N cr -965.9814 Ib/in) 
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Khot'» gloss — epoxy shell ( — 4-0. -4-0.0) xiba r2 0. 1 O — 05. 1 ,n~ 1 2 
W 1 — max vs axiol load (both normalized) — C — 4 S.C. 



a. w-| max vs normalized axial load p s = Xg/Xe 


Khot’s gloss — epoxy shell (-4-0. 4-0,0) xibor2»“0. 1 D — 05.m» 1 ( n» 1 2 
end — shortening (total) vs axiol load (both normalized) — C — 4- S.C. 



end-shortening (total) / 0.779637 


b. end-shortening vs normalized axial load p s = X S A. C 

1 8 Load vs deformation plots for C-4 boundary conditions 
(u=u 0 , v=w=W, x =0 ; N c j=-965.9814 Ib/in) 
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b. W1(x). W2(x) 



d. SS-4 B.C. Xs= 0.538999 (n=1 2) 



o. Wv + W0(x) 



b. W1(x). W2(x) 



e. C-1 B.C. X s = 0.537957 (n =12) 




b. W1(x), W2(x) 



e. Response Mode 


f. C-2 B.C. X s = 0.53981 9 (n-1 2) 


Fig. 19 Response mode shapes at the limit points for various boundary conditions (continuation) 
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g. C-3 B.C. X s = 0.538657 (n=1 2) 


4 W +W 
| v 0 



K(mon)«i/2 


a. Wv + W0(x) 



c. Response Mode 


h. C-4 B.C. Xs= 0.539966 (n =12) 


Fig. 19 


Response mode shapes at the limit points for various boundary conditions (continuation) 
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c. Post-buckling Component Modes 


a. SS-1 B.C. X ^ - 0.292266 (n=2) 



c. Post-buckling Component Modes 


b. SS-2 B.C. C= 0.292659 (n =2) 


Fig. 20 Bifurcation buckling mode for various boundary conditions 



c. SS-3 B.C. X ' = 0.536362 (n =12) 
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* * * 
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i(mo>c)wL/2 


a. Prebuckling Shape 



c. Post-buckling Component Modes 


h. C-4 B.C. x" 1 - 0.539966 (n=1 2) 


Fig. 20 


Bifurcation buckling mode for various boundary conditions (continuation) 
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